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Synopsis 


Porous silicon was accidentally discovered in 1956 during the anodic etching of 
crystalline silicon in aqueous HF solution. It is a brittle, spongy material consisting of 
an intricate network of nanometer-size pores and silicon crystallites. The formation of 
a variety of pore morphologies in preference to uniform etching of silicon has evoked 
a great deal of interest in the scientific community. In 1990, high porosity (60-80%) 
samples of porous silicon were fabricated, which exhibited a surprisingly strong red 
luminescence on illumination with light of wavelength A = 514.5nm [1], Ever since 
this discovery, porous silicon has been the focus of attention of several workers in the 
field of semiconductor physics [2.3]. 

Crystalline silicon is an indirect gap material with a band gap of 1.1 eV, and 
consequently, the momentum conservation rule inhibits radiative transitions between 
the conduction band minimum and the valence band maximum. In porous silicon, the 
observed photoluminescence is in the range 1.3 - 2.2 eV, and has an external absolute 
quantum efficiency of 1-10%. This unexpected behavior has been explained in terms 
of the quantum confinement model [l] which posits an enhanced energy gap because of 
the increased kinetic energy due to carrier confinement in nanometer sized crystallites. 
However, this simple explanation is hard put to explain several vagaries observed in 
the photoluminescence behavior. Chapter 1 constitutes a brief introduction to the 
porous silicon problem, where we review the porous silicon formation process and 
major theories proposed to explain the photoluminescence behavior. 

Porous silicon formed under different anodization conditions exhibits a variety of 
rich and complex structures similar to those observed in diverse growth phenomena 
such as electrochemical deposition, viscous fingering and bacterial colony formation. 
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The presence of such varied behavior has led to the view that some interesting physical 
processes must be at work, besides the electrochemistry of silicon dissolution. 


Chapter 2 describes a drift diffusion 
model [4] which is successful in gener- 
ating the various morphological classes 
reported in experimental growth liter- 
ature. These are (i) fractal dendritic 
growth as seen in diffusion limited aggre- 
gation (DLA), (ii) homogeneous dense 
branching morphologies with a stable 
circular growth front, (iii) stringy mor- 
phologies. and (iv) inhomogeneous den- 
dritic morphologies. We propose that 



Fig.1: Phase diagram of various morphologies 
obtained by tuning two parameters (i) depletion 
width (ii) drift length (in lattice units) 


the four classes result from the compe- 
tition between two antithetical scales - a diffusion length which tends to keep parti- 
cles away from the growing cluster, and a drift length which drives particles towards 
surface inhomogeneities. A phase diagram depicting the transitions between morpho- 
logical classes with varying scales is shown in Fig. 1. 


In Chapter 3 we extend the drift diffusion scheme to silicon anodization in HF, 
by means of a diffusion induced nucleation model [5]. A competition between the 
semiconductor depletion layer and electric field enhancement due to surface irregu- 
larities induce random removal of silicon atoms (i.e., an inverse nucleation process). 
This scheme successfully reproduces experimentally observed phenomena such as (i) 
high porosity structures similar to samples exhibiting visible photoluminescence, (ii) 
a constant rate of growth, (iii) the dependence of the rate of growth on the anodiza- 
tion potential, and (iv) electropolishing of silicon in the high potential limit. Further, 
the effect of quantum confinement on porosity can also be demonstrated within this 
model. 


The photoluminescence (PL) spectrum of porous silicon is asymmetric on the 
energy scale, with a large full width at half maximum (FWHM) of 300 meV. 
Chapter 4 presents a model for the spectral line shape [6]. Porous silicon is 
considered to be a collection of disordered silicon crystallites having a Gaussian 
distribution of sizes in the range 2-4 nm. Based on effective mass theory, 
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analytical expressions for the photolu- 
minescence spectrum of such a distribu- 
tion is obtained. This can explain both 
the asymmetry in the PL spectrum and 
the large FWHM (see Fig. 2). Based on 
the peak position of the calculated spec- 
trum the exciton binding energy of the 
electron hole pair is estimated. This is 
found to be in better agreement with our 
understanding than those estimated by 
earlier electronic structure calculations. 

The low energy tail shown by some pho- 
toluminescence spectra is also explained using the Lifshitz argument. 

The simple quantum confinement model was hard put to explain several experi- 
mental observations related to the temperature and pressure dependence of lumines- 
cence from porous silicon such as: 

1. The PL intensity shows a maximum with temperature, the peak being in the 
range 50Iv - 150K 

2. The luminescence decay time falls by an order of magnitude or more as the 
temperature is increased from ~ 10K to room temperature. 


Fig.2. The porous silicon photoluminescence spectrum. 
The solid line is an experimental spectrum. The broken 
line represents a theoretical fit based on Ref. 6 


3. The PL intensity falls steeply with pressure. 

4. The luminescence decay time decreases systematically with emission energy in 
rhe range 1.4 - 2.5 eV 

5. The PL peak position exhibits an anomalous behavior (both blue and red shifts) > 
with temperature. 

6. The PL peak position exhibits an initial blue shift and a subsequent red shift; 
with pressure, when the pressure transmitting medium is alcohol. 

Chapter 5 describes a unified model to explain all the above observations within; 
an entirely analytical framework [7]. Based on conductivity studies, we propose; 
that carrier hopping in silicon nanocrystallites follows a Berthelot type temperature; 
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Fig. 3. The temperature dependence of the luminescence 
decay time in porous silicon. The model of Ref.7 yields 
excellent fits to the data reported by different groups. 


Fig.4. The temperature dependence of intensity in por- 
ous silicon luminescence. The solid line is a theoretical 
fit based on Ref. 7 to the data reported in Ref. 3. 


dependence, i.e. it increases exponentially with temperature. The competition be- 
tween the carrier hopping and an Arrhenius type radiative recombination is shown 
to explain the temperature and pressure dependence of the luminescence lifetime and 
intensity (see Figs. 3 and 4). We then extend this model within the quantum con- 
finement scheme and successfully explain the luminescence lifetime dependence on 
emission energy and the observed temperature and pressure dependencies of the PL 
peak position. 

Porous silicon has a large surface area. Several workers have contended that the 
luminescence originates from the molecular complexes found on the porous silicon 
surface [2]. Alternately, the disordered surface may give rise to several trap states 
which can play an intermediate role in the photoluminescence process. In this context, 
it is necessary to characterize the porous silicon surface. A common tool used for this 
purpose is infrared vibrational spectroscopy. The vibrational modes of the Si-H bond 
are known to be correlated to the electronegativity of the surrounding environment. 
We define a renormalized electronegativity based entirely on theoretical atomic scales. 
This new scale is successfully correlated to the vibrational frequencies of a given Si-H 
or Si- Cl bond in a particular environment [8]. This model is outlined in Chapter 6 
and used to systematize the assignment of vibrational frequencies observed in porous 
silicon. 

In conclusion, we have presented a unified model to explain the vast body of 
experimental data in porous silicon. A simulation model for the formation of diverse 
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pore morphologies in porous silicon has also been developed. We show that the 
luminescence properties are modified by the disordered nature of the system as well 
as carrier hopping on or between crystallites. The optical characteristics are thus 
seen to correlate with the mesoscopic nature of porous silicon. Chapter 7 presents 
a brief discussion of our results and suggests guidelines for future work. 
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Chapter 1 


Introduction 


Silicon is - arguably - the most widely studied element in history. The technologi- 
cal applications of this material have revolutionized modern civilization and it has 
emerged as the workhorse of the semiconductor industry. The latter half of this 
century may rightly be termed the “silicon age” - so ubiquitous is the influence of 
this element in our day-to-day existence. On the other hand, silicon has also helped 
in furthering our understanding of traditional solid state physics. It has been in- 
tensely studied by the band structure theorists, being used as a test case for every 
new scheme, and any variation in methodology whether major or minor. 

Nearly fifty years into silicon research, this thesis which purports to carry out fur- 
ther theoretical studies on silicon has been motivated by the possibility of extending 
the dominance of this material to an arena hitherto dominated by III-V semiconduc- 
tors - namely, optoelectronics. Crystalline silicon has an indirect band gap of 1.17 
eV, and the momentum conservation rule inhibits radiative transitions between the 
valence band maximum and the conduction band minimum. In 1990, photolumi- 
neseenee in the range 1.3 - 2.2 eV with an external absolute quantum efficiency of 
1 - 10% was observed from a novel material called Porous Silicon [1]. Since then, 
several review articles [2-5], conferences [6-9], and well over a thousand publications 
have been devoted to this interesting optical behavior. Theoretical studies of the 
formation of porous silicon and its surprising optical behavior constitute the subject 
matter of this thesis. 
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Introduction 


1.1 Porous Silicon 


Porous silicon (PS), as the name sug- 
gests, consists of a brittle silicon skele- 
ton permeated by an intricately woven 
network of nanometer sized pores. Dur- 
ing anodization of silicon in hydrofluo- 
ric acid (HF), material removal does not 
always occur in the form of a uniform 
etch, but through the formation of mi- 
croscopic trenches on the surface (see 
Fig. 1.1). Eventually these trenches de- 
velop branches and get interconnected, 
leading to the formation of a pore net- 
work. The resultant silicon structure is 
seen to consist of irregular columns and 
spheres with diameters in the range 1-10 



Figure 1.1: Trench formation during silicon 
anodization in HF acid. 


nm. 


The central problem in porous silicon research today is discovering the origin of 
light emission from silicon. In order to maximize the luminescence yield and tune 
the emission parameters such as frequency and decay time, a better understanding 
of the carrier dynamics associated with the light emission process is needed. The 
diverse pore morphologies found in PS have puzzled the semiconductor community 
for the last four decades, and continue to elude a concrete theory. The eventual 
goal of all research is to fabricate PS under controlled conditions leading to desired 
characteristics needed to make an optoelectronic device, which could ultimately be' 
integrated into a “chip.” Hence, this novel material opens up a new vista of truly 
interdisciplinary research involving semiconductor physics, non-equilibrium growth 
processes, electrochemistry etc. 


The present thesis aims to explore - at a theoretical level - the formation mech- 
anism as well as the photoluminescence behavior of porous silicon. The remaining 
sections of this chapter present a brief 'overview of existing theoretical models of PS 
photoluminescence. Microscopic attempts to validate these models are reviewed in 
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See. 1.4. The issues that can be resolved within these models are spelt out, and the 
problems that remain unresolved are discussed. 

Chapters 2 and 3 provide a short review of the formation problem, and develop a 
computational model that successfully simulates several characteristics of the porous 
silicon formation process. Chapters 4 and 5 study the optical properties of PS, 
and present models for the PL line shape, intensity and decay dynamics. Chapter 6 
proposes a phenomenological scheme for the assignment of vibrational modes observed 
in PS using infrared spectroscopy, which has gained ground as a powerful tool for the 
characterization of the PS surface. Chapter 7 summarizes our results and provides 
guidelines for future research. 


1.2 Formation of Porous Silicon 

Ever since the discovery of porous silicon in 1956 at the Bell laboratories by Ulilir 
[10], the causative mechanism associated with its formation process has been of some 
interest to the researchers. Unlike theories of electronic structure and photolumi- 
nesconce, theoretical studies on porous silicon formation span nearly four decades. 
In spite of this, formation is perhaps the least understood of the several enigmas 
posed by this material. Various models have been proposed to explain the formation 
of the observed pore patterns. These range from purely electrochemical theories to 
stochastic simulations representative of physical phenomena. The discovery of visible 
photolumincsccnce in 1990 has also served to enhance the interest in PS formation 
studies, as the synthesis of PS under controlled conditions is an essential pre-requisite 
for a better control over its optical properties. 

Porous silicon is normally synthesized in an electrolytic reaction where the silicon 
substrate constitutes the anode and an acid resistant metal such as platinum is used ' 
as the cathode. The electrolyte consists of ~ 49% HF, diluted in de-ionized water ; 
or ethanol at ratios ranging from 1:1 to 1:5. Currents in the range 10 - 30 lnA/cm 2 ! 
are passed through the electrolyte over a period of 5 - 60 minutes. A porous layer 
is formed on the surface of the silicon substrate, its thickness increasing uniformly j 
with time. PS layers of thickness upto 500 /xm have been fabricated by this method, j 
Note that in the absence of any current, the acid would have eaten away the silicon j 
uniformly at a rate of a few nm/hr. This process, distinct from the anodization ! 
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procedure, is termed “passive etching.” 

Transmission electron microscopy (TEM) studies of porous silicon have been re- 
ported by several workers [11,12]. They reveal that porous silicon consists of isolated 
• crystallites and/or chains of crystallites linked by narrow silicon walls. The crystal- 
lites are found to be in the size range 1-10 mil. This implies clusters of about 1() 2 — K)' 1 
atoms which retain the diamond like lattice structure of bulk silicon. The pore struc- 
ture exhibits a variety of morphologies depending on the experimental conditions 
governing the anodization process, namely, 

• doping type 

• doping level 

• anodization current 

• electrolyte concentration 

A p-type silicon substrate forms a uniform and highly interconnected network of 
random pores, resulting in a spongy structure. In n-type, the pores are formed only 
under illumination and are relatively linear. For degenerately doped cases (both p + 
and n + ), the geometries are similar to the interconnected p-type structure, but with 
shorter and wider pores. At higher current densities, the pores widen further. In 
n-type, the linear pores tend to become pipe like in this case. At very high currents, 
silicon is uniformly etched away, leading to electropolishing. This can also occur at 
very low HF concentration. 

Besides the electrochemical etching of crystalline silicon, several other techniques 
have been employed to fabricate Si nanocrystallite layers. Microwave plasma decom- 
position of Si-H 4 silane gas onto a quartz substrate gives rise to spherical nanocrys- 
tallites in the size range 2-7 nm. [13]. Laser breakdown of silane also results in silicon 
nanocrystallitcs with an oxide coating of ~1.6 nm. thickness [14]. From Si02, Si 
nanoparticles have been obtained by rf-magnetron co-sputtering techniques [4]. This 
method has also been employed to obtain Ge nanocrystallites. Another successful 
approach to fabricate porous silicon like layers is the high frequency spark erosion 
[15] of a silicon substrate using a tungsten anode. All these methods result in silicon 
nanocrystallites capable of luminescing in the visible range. Recently, light emitting 
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Si-Si0 2 superlattices have been grown by molecular beam epitaxy [16,17]. However, 
in this thesis, we will locus our attention on the anodization of silicon in HF. 

The most common characterization parameter employed to describe porous silicon 
morphology is the sample porosity. Assuming a silicon wafer is completely converted 
to porous silicon, the porosity p is calculated as 

mass of Si wafer - mass of porous film 

p — x 100 

mass ol Si wafer 

The preferential pore formation during anodization in comparison to a uniform 
etch raises the following questions, “Is it a chemistry related phenomena caused by 
competing reaction mechanisms, or is there more physics to it? Do the semiconduct- 
ing properties of the system play any role in the process?” The earlier phenomenolog- 
ical models have focused on the chemistry of dissolution as the determining factor for 
pore formation. Then the attention shifted to possible surface effects on the electric 
field distribution in the Si substrate. The role of the semiconducting properties of 
the system in effecting an electric field redistribution assumed importance as a major 
factor. With the better understanding of non-equilibrium growth phenomena which 
lms been achieved in the recent past, several statistical models based on computer 
simulations have also been formulated to explain porous silicon formation. 

A satisfactory theory of PS formation must be able to (i) understand the formation 
of pores in preference to a uniform etch, (ii) explain the observed transition between 
pore morphologies obtained under various experimental conditions, and (iii) predict 
features such as the rate of growth of the film and the influence of initial conditions, 
say, the substrate roughness, on the porous layer. The eventual aim, as mentioned 
earlier, is the fabrication of a PS device under controlled conditions. Hence, mod- 
els which can computationally simulate experimental configurations and predict the 
resultant device geometry and structure will be desirable. In chapters 2 and 3 we 
evolve a generalized computational model, premised on a possible underlying semicon- 
ductor phenomenology. This model successfully simulates the diverse morphological 
patterns described earlier and is successful in making experimentally verifiable pre- 
dictions such as the rate of growth of the porous film and its dependence on the 
anodization potential. 



6 


Introduction 



Figure 1.2: A typical photoluminescence spectrum. 

1.3 Luminescence from Silicon 

It is the visible photoluminescence (PL) of silicon and related optical properties that 
have invited maximum scientific attention. Porous silicon starts luminescing at a 
threshold porosity of around 45%, with the emission being in the range 1.2 - 2.8 eV, 
depending on the sample porosity. High porosity (> 60%) samples have an external 
absolute quantum efficiency [18] of 1 — 10%. A typical PL spectrum is depicted in 
figure 1.2. The spectrum is asymmetric about the peak, with a shoulder on the high 
energy side. The full width at half maximum (FWHM) is about 300 - 400 meV. 
Three bands have been identified in the emission spectrum of porous silicon. They 
are: 

• The dominant red band with energies in the range 1.2 - 2.3 eV. This accounts for 
roughly 97% of the emission and has a slow luminescence decay with lifetimes 
from 10-* to 10~ 3 s. 

• The weak blue band with energies in the range 2.5 - 2.8 eV. This constitutes 
roughly 3% of the emission intensity and has a decay time in the nanosecond 
range. 

• The infrared band with energies ranging from 0.8 to 1.3 eV. 

Though the photoluminescence from porous silicon presents an interesting puzzle 
to the scientific community, the technological importance of this material is based on 
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its electroluminescence behavior - i.e., the emission of light when a current is passed 
through it. Porous silicon can electroluminesce both in aqueous electrolytes (wet 
electroluminescence) and solid state structures such as Au/PS/Si layers (dry electro- 
luminescence). The luminescence efficiency for wet electroluminescence (~ 10 -3 %)is 
substantially higher compared to dry electroluminescence (~ 10~ 5 %). However, for 
practical device fabrication, liquid-solid junctions are of little use. Further, wet elec- 
troluminescence is also affected by chemical reactions taking place on the porous sil- 
icon surface. In aiiy case, the efficiency of electroluminescence is at best a thousand 
times less than that of PL. The spectra observed in photoluminescence and electro- 
luminescence are remarkably similar. This suggests that the radiative mechanisms in 
both processes are of similar origin. 

A variety of hypotheses have been advanced to explain the observed optical behav- 
ior of PS. The major models, namely, (i) the Quantum Confinement model, (ii) the 
Chemical Luminescence model, and (iii) the Hybrid model invoking radiative recom- 
bination via intermediate states are briefly reviewed below: 


1.3.1 The Quantum Confinement Model 


The carriers in high porosity samples of porous silicon are confined in microcrystallites 
whose typical size d is in the range 1-5 nm. This can be modeled by a simple picture 
- the particle in a box. From the uncertainty principle d A p ~ h, the energy of the 
carriers relative to the band edge is given by the relation 


(Ap) 2 = h 2 
2m 2 md 2 


( 1 . 1 ) 


This applies to both holes and electrons and a larger energy gap is plausible on these 
grounds. 

We can extend the particle in a box picture using the well known effective mass 
approximation. The upshift of the conduction band can then be expressed as 

where ?n.r(?U/,) is the transverse (longitudinal) effective mass, and d the confining 
length. For a 3nm wire, the conduction band upshift is estimated to be 252 meV 
[19]. Here we have assumed bulk values for the transverse (0.19m e ) and longitudinal 
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Figure 1.3: The effect of confinement on the band structure of silicon, (a) Crystalline 
silicon with an indirect gap of 1.17eV (b) On confinement, the conduction band (CB) 
and valence band (VB) edges are pushed apart (see arrows) and become pseudo- 
direct (c) The density of states n(E) in a realistic nanocrystallite, with a U shaped 
distribution of surface states S near the band edges. The peak near the mid-gap 
region is due to dangling bonds (DB). 


(0.92m e ) effective masses, m e being the free electron mass. Because of the valence 
band degeneracy, the variation in the hole kinetic energy depends on the direction 
of the wire growth. Assuming growth in the (001) direction, the effective mass ap- 
proximation yields the downshift of the valence band to be 371 meV [19] which is 
about one and a half times the conduction band upshift, the net enhancement in the 
gap being 0.G cV. Thus energy gaps in the range 1.3 - 2.5 eV are possible. This can 
quantitatively account for the red luminescence. 

Figure 1.3a depicts the bulk silicon band structure with a well defined indirect 
gap. The upshift due to confinement is depicted in Fig. 1.3b. A crystallite, however 
has discrete energy levels. The level spacing is ignored in the figure. Further, these 
calculations result in nearly flat bands at the valence band maximum and the con- 
duction band minimum. This too can be understood on the basis of the uncertainty 
principle. We see that A k ~ 1/d is large. The indirect gap is thus ‘smeared’ out and 
an almost direct gap may result. These arguments form the core of the basic quantum 
confinement model. Several microscopic electronic structure calculations have been 
carried out to verify this model. We shall discuss them in the next section (Sec. 1.4). 
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1.3.2 Chemical Luminescence Models 

Visible luminescence from porous silicon has been found to be extremely sensitive to 
surface treatments. This lias prompted some workers to contend that luminescence 
originates not from the silicon crystallites per se, but from chemical complexes found 
on the surface of porous silicon. Armed with an extensive knowledge of “ancient” 
literature, they came up with several molecules capable of luminescing in the observed 
spectral range. A particular favorite was ‘siloxene’- a complex of silicon, hydrogen and 
oxygen. This contention is supported by a comparison between the salient features of 
PL from both porous silicon and siloxene [20], A review of the properties of siloxene 
is found in the article by Stutzmann and coworkers [21] who notes that 

(i) annealed siloxene radiates in the same spectral range as porous silicon (600 - 
800 nm) with long tails extending into the deep infrared. 

(ii) The PL spectrum of annealed siloxene is broad, with a band width much 
larger than that of porous silicon. All the observed PL spectra fall within this range. 

(iii) The temperature dependence of PL intensity for T > 100K in both materials 
is well described by the relation 

7 ^ = 1 +exp(|) 

with a characteristic temperature To between 50 K and 90 K. 

(iv) The luminescence intensity is proportional to the excitation intensity over 
six orders of magnitude at both low and high temperatures for porous silicon and 
siloxene. 

(v) 'Phe luminescence decays are non-exponential and become faster with increas- 
ing energy. Both nanosecond and microsecond to millisecond decays are observed in 
siloxene, jus also with porous silicon. 

(vi) Both materials exhibit similar photoluminescence excitation spectra. 

'The siloxene model has been disputed by some other workers [22] who report PL 
in freshly prepared porous silicon samples which exhibit infrared vibrational spectra 
lacking all oxygen related modes. Such a sample cannot have a very high intensity 
siloxene based luminescence. Besides this, the luminescence induced by soft x-rays in 
porous silicon (peak ~ 750 nm) differs substantially from that produced in siloxene 
(peak ~ 520 nm) [23]. 
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Besides siloxene, some of the other molecules that have been hypothesized to 
l)e the luminescing agents in PS include polysilane chains [24,25] and non-bridging 
oX ygen-hole complexes [26]. The polysilane hypothesis has also derived support from 
theoretical calculations carried out by various groups, which will be discussed in the 
ne xt section (Sec. 1.4). 

1.3.3 Recombination via Intermediate States - The Hybrid 
Model 

The tug-of-war between the quantum confinement model and the chemical complexes 
model has led to to a “hybrid model” where both the interior and the surface are 
involved in PL. Consider the case where an electron is radiatively excited into the 
conduction band in the interior of a crystallite. For a crystallite of size 3 nm, over 
3(j% of the atoms will be on the surface, which can give rise to several surface states. 
Besides this, there could be dangling bond levels which fall in the mid-gap. The 
density of states in a nanocrystallite is represented in Fig. 1.3c. The surface states 
give rise to a U shaped distribution within the gap, and the dangling bond levels fall 
nC ar the mid-gap legion. The carrier has a high probability of reaching the surface 
and getting trapped there. It is then possible for surface related mechanisms to 
influence (.lie radiative process. 

Kanemitsu and coworkers [14] have proposed a simple phenomenological model 
based on their experimentally acquired data. The absorption edge for crystallites of 
pizes 9, 3.5 and 2 nm is observed at wavelengths 1200, 1000 and 850 nm respectively. 
The crystallite sizes were determined by TEM and analysis of the Raman spectra 
0 f the samples. The blue shift on size reduction may clearly be attributed to the 
quantum confinement effect. The PL peak energy is however invariant and located 
at 700-800 nm. They contend that the latter is thus an intrinsic effect common 
to all crystallites. The photogeneration of carriers is hypothesized to take place 
in the core of the crystallites and hence is size dependent. All crystallites have a 
large surface to volume ratio with a similar surface chemistry in terms of H and 
Oil complexes. The recombination of carriers take place on the surface and hence 
pL is size independent. In another work [27] it is proposed that the red band PL 
has a slow decay and is from surface localized states while the blue PL band has 
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a fast decay and originates from the crystallite core. An implication of the model 
is that PL efficiency for smaller crystallites is enhanced since both the surface to 
volume ratio and the carrier transfer rate from the core to the surface increase with 
a decrease in the crystallite size. The PL intensity has a non-monotonic temperature 
dependence. At low temperatures, carrier diffusion to the surface is inhibited, hence 
PL intensity is low. At high temperatures, nonradiative recombination in the near 
surface regime increases and the PL intensity is lowered. Hence this model can also 
explain the peaking of PL intensity at around 100 K, which is an experimentally 
observed phenomena. 

More detailed arguments for the surface states mechanism have been advanced 
by some workers [28,29] who also believe that radiative recombination occurs via the 
surface states on the silicon nanocrystallite. They propose that carrier recombination 
can occur in three ways (i) Carriers trapped in the bound states can recombine among 
themselves. These corresponds to the lower emission energies, (ii) Carriers from the 
band can recombine with carriers in the gap levels and (iii) there could be band to 
band recombination. They assign the infrared band to the capture of electrons by 
dangling bond states which fall in the mid-gap. The red band is mediated by the 
surface states in the gap which lie closer to the conduction band edge (see Fig. 1.3c). 
The fast blue emission is assigned to transitions from the conduction band to the 
valence band in the core of the crystallite. Thus the carriers can “choose” between 
several paths available for radiative recombination. Hence they term this the “smart 
quantum well” [30]. 

1.4 Microscopic theories 

The enigma of porous silicon is its visible luminescence. One would like to understand 
the photoluminescence (PL), its excitation and decay characteristics. In the present 
section, we review attempts to explain PL based on explicit electronic structure 
calculations. These attempts must be distinguished from a host of phenomenological 
and qualitative theories which pervade the field, some of which were discussed in the 
previous section. 

The calculation of the electronic properties of porous silicon is premised on an 
assumed geometrical structure and herein lies a major difficulty. Porous silicon is a 
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Figure 1.4: The levels of structural disorder in porous silicon may be staged thus. 
The left side represents the idealized structures normally used in electronic structure 
calculations. The right side is the more realistic view (a) Nanocrystallite at a micro- 
scopic level. The perfect crystal symmetry is disturbed by the presence of vacancies 
and multivacancies (not shown) accompanying dangling bonds, (b) The nanocrystal- 
lite is idealized in the form of a column (pillar) or dot (sphere). The surface of such 
a crystallite is highly irregular. The large surface area supports silicon hydride and 
oxide complexes and a large number of dangling bonds, (c) The columns or dots (not 
shown) are not of uniform diameter. They come in a distribution of sizes, (d) A cross 
sectional view of an idealized, regularly spaced and cylindrically or spherically shaped 
nanocrystallite. The more realistic view is one of irregularly shaped nanocrystallites 
forming a silicon skeleton against a backdrop of an intricate network of pores. 
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complex material which one would characterize by a hierarchy of disorders. Fig. 1.4 
illustrates this complexity. The left side of Fig. 1.4 represents the idealized view 
adopted in most electronic structure calculations. The right side depicts the grim 
complex reality. 

At the microscopic level one has nanocrystallites. The internal arrangement of 
silicon atoms is believed to be the same as the diamond like structure of crystalline 
silicon [31-33]. But one must accept a fair number of native defects: vacancies, diva- 
cancies, multivacancics and concomitant dangling bonds with bond rearrangements. 
The nanocrystallites occur in two forms: (i) columnar (also termed pillar or rod) or 
(ii) dots (spheres) or an admixture of both. However the surface is not all smooth. It 
is highly irregular and the surface area is large. The surface has dangling bonds and 
supports molecular complexes such as [Si n H m ] polysilanes, SiC >2 and siloxene deriva- 
tives. Further, the nanocrystallites may come in a variety of shapes and a range of 
sizes from 1-10 lira. Finally, these are not arranged in a regular fashion, but are ran- 
domly arranged with intervening spaces. One has a brittle, spongy silicon skeleton - 
hence the name porous silicon. 

A calculation of the electronic properties which encompasses all levels of disorder 
is not feasible. The following points can be used to guide the calculations: (i) The 
local structure of the nanocrystallite is akin to crystalline silicon with lattice constant 
ao = 5 . 43 A. (ii) The aspect ratio (height to radius) for the columnar structure is 
approximately 1000. (iii) The surface area of the nanocrystallite is large: 1-10 m 2 / 
cm 3 , (iv) The porosity of the structure is high with p « 80%. 

The calculations assume an idealized structure such as the ones depicted on the 
left side of Fig. 1.4. Further, one may use periodic boundary conditions. In this case, 
Bloch’s theorem holds and one can talk of the “band structure” of porous silicon. 
The other mode in which the calculation can be carried out is to use a cluster of 
TV silicon atoms with N b dangling bonds terminated by hydrogen. If one is using 
four atomic orbitals per silicon site, one obtains 47V + N b energy levels. The highest 
occupied (lowest, unoccupied) energy level corresponds to the valence band maximum 
(conduction band minimum). The corresponding wave functions are termed HOMO 
(highest occupied molecular orbital) or LUMO (lowest unoccupied molecular orbital). 

The methodology used can be classified as either semi-empirical or first principles 
based. Most of the semi empirical attempts are based on the tight binding scheme. 
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The advantages of this approach are: (i) it is easy to implement; (ii) large clusters 
of sizes comparable to the actual nanocrystallites can be studied (storage efficient); 
(iii) the algorithms run faster (time efficient). The disadvantage is the semi-empirical 
character itself, the reservation being that one can obtain a desired value by simply 
fine tuning parameters, or by adding a fresh one. The parameters of the Hamilto- 
nian are determined with reference to some property or properties of the bulk host 
silicon. Thus optimized, the parameters are employed unchanged for porous silicon 
calculations. 

The first principles methodology is based on the Kohn-Sham local density ap- 
proximation (LDA) [34,35] using pseudopotentials. No adjustable parameters are 
employed in this scheme. LDA predictions for the ground state properties of silicon 
have been verified by a variety of works and further, heralded as a major triumph 
of band structure methodology [36-38]. The lattice constant ao = 5. 43 A, the bulk 
modulus (0.99 x 10 12 dynes/cm 2 ), the total valence energy (-7.919 Ry), the cohesive 
energy and the frozen ( T -4 0) phonon frequency are all reproduced to an accuracy 
of « 1%. In spite of such an impressive list of achievements, there is a major problem 
with the LDA. The predicted band gap E g of crystalline silicon is 0.6 eV, almost 
half of the observed value of 1.14 eV at T — 0. This corresponds to luminescence 
(if at all) in the far infra-red (A = 2 x 10 4 A). Given the fact that a primary aim of 
porous silicon electronic structure calculations is to explain the visible PL, the band 
gap underestimation by LDA is unfortunate. In LDA based investigations of porous 
silicon, the bulk E 0 is engineered by upshifting the conduction band in various ways. 

Read and coworkers have performed LDA based calculations on columnar wire 
geometries with a uniform rectangular cross section [39] as well as undulating columns 
which can be visualized as a periodic array of overlapping spheres [40]. The surface 
is passivated with H atoms. Periodic boundary conditions were employed so as to 
obtain a band structure. For a rectangular wire of diagonal thickness d = 2.3 nm, a 
direct gap of 0.78 eV was obtained, which is larger than the reported LDA values of 
0.5 - 0.6 cV [37]. The band upshift A E u could be expressed ns 

A E u = C/d 2 

where C is a constant. Using this expression, they calculate the band gap for a 3 
nm crystallite to be close to 1.7 eV, which is larger than the observed PL peak of ~ 
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1.5 cV. Consequently, they have to invoke an unreasonably large excitonic binding 
energy of almost 200 meV. In chapter 4 we show that the inherent disorder in the 
PS system can lower the observed peak energy, and hence one need assume only 
physically reasonable exciton energies of about 60 meV. 

In a similar first principles study employing square columns, Buda and cowork- 
ers [41] also obtain a direct gap E g which is seen to vary with the crystallite size 
as 

E,~E ,( «) + ! + ! 

where rq , c 2 are constants. This enhancement of the bandgap with decreasing d is 
also observed by Olmo and colleagues [42] in their LDA based calculations. They also 
find that the oscillator strength for transitions between the valence band maximum 
and the conduction band minimum is enhanced by 2 to 3 orders of magnitude as a 
result of confinement, though it is found to be a non-monotonic function of the size. 

Semi-empirical calculations have also been carried out on free standing silicon 
wires [43] in the size range 7-27 A. A complementary geometry of circular and square 
pores with interpore space occupied by Si has also been employed by some workers 
[44] (see Fig. 1.4d). All of them yield direct or nearly direct band gaps in the range 2.0 
- 3.0 eV. The oscillator strength has also been found to increase with d [45], justifying 
light emission. Thus the quantum confinement model is theoretically validated for 
periodic arrays of Si crystallites. 

Calculations have also been carried out for clusters of Si atoms passivated with 
II, within the LDA [46]. They report a decrease in the energy gap E g with increasing 
cluster size d. Further, the oscillator strength is seen to exhibit a a non-monotonic de- 
pendence on d, similar to the results of Ohno and coworkers [42] which were described 
earlier. 

Some workers [47-49] employ a non-orthogonal tight binding framework to treat 
nearly spherical (closed shell) silicon clusters. They chart the gap dependence on size 
up to d = 4.3 nm, (2058 silicon atoms). The attempted fit does not follow a 1/d 2 
law, but 

E g (d) ~ l/d 7/5 

The calculated gap is compared with the experimental results on silicon nanocrys- 
tallites coated with hydrogen [50]. The PL emission energy is the difference of the 
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gap and the Coulomb energy (-3.752 e 1 fed) as estimated by Kayanuma [51]. Further, 
they evaluate the radiative recombination lifetime r r , and find it to be in the range 
10~ 4 - io- 6 sec. for E g > 1.8 eV and 1CT 1 - 10" 5 sec. for 1.6 < E g < 1.8 eV. 
The dependence of r r on the PL emission energy is found to compare favorably with 
experiment. 

A noteworthy attempt to model the disorder in PS is due to Sawada, Hamada 
and Ookubo [52] who employ a semi-empirical tight binding methodology to study 
a two dimensional randomly porous cluster. The number of atoms in the cluster is 
around 450. A fraction (say, 100 atoms) are removed at random. The dangling bonds 
are terminated by hydrogen. The gap for a finite regular cluster of 350 atoms is 1.51 
eV while that for 450 atoms is 1.39 eV. The difference of 0.12 eV can be ascribed to 
confinement effects. On the other hand, for the porous structure where 100 atoms are 
removed at random from the 450 atom cluster, the gap is 2.06 eV. The additional 0.55 
eV enhancement is a disorder induced effect. Thus porous silicon has an enhanced 
gap due to both confinement and disorder. 

Calculations have also been carried out in support of the chemical luminescence 
model. The electronic structure methods of quantum chemistry have been exten- 
sively employed for this purpose. Deak and coworkers [53] have considered siloxene 
structures which are planar [Si 2 H(OH)] 3n , linear [Si2H 2 0] 3n , and ring-like [Si6H 6 03]„ 
and substituted ring-like siloxenes [Si6H 6 _ x (0H) x 03] n . It is found that the planar 
structure is metastable and will transform to a linear or ring like structure on heat 
treatment. Thus the surface of porous silicon is more likely to be covered by ring like 
siloxene and its other derivatives. The band is indirect and has a value 2.84 eV for the 
unsubstituted x = 0 case and 1.82 eV for the fully substituted x = 6 case. An analy- 
sis of the HOMO charge density reveals that it is high along the rings. The authors 
thus speak of a “chemical” quantum confinement in contrast to the “geometrical” 
quantum confinement. 

In another quantum chemical study, Takcda and coworkers [54] consider a bridge 
model Sixo-Cbi5-Si 4 H 8 -SiioCbx 5 . Here Si 4 H 8 is the oligosilane bridge. Cb stands for 
a capped bond. The Six 0 silicon cluster is terminated not by H, but the capped 
bond which is a Is orbital with electron affinity identical to Si. The gap 3.3 eV is not 
identified with PL emission, but with PL excitation. On excitation a fresh calculation 
employing configuration - interaction is carried out. The excited structure is once 
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again allowed to geometrically relax. The gap is now in the range 1.2 - 2.1 eV with 
the LUMO 50% .s-like while the HOMO is 60% p-like. Given these symmetries, the 
optical transition matrix is non-negligible for the 1.2 - 2.1 eV de-excitation. This 
is identified with the PL emission. After de-excitation, it is found that the earlier 
structure is once again favored. The work thus attempts to explain why the PL 
excitation and emission energies generally differ by at least 0.5 eV. 

Allan, Dclcruc and Lannoo [55,56] have calculated the electronic structure of 
linear a bonded silicon atoms with H termination. They find that the radiative 
lifetime is independent of the size and falls in the nanosecond range. On the basis 
of the calculated energies and radiative lifetimes, they conclude that polysilanes may 
contribute to the blue band PL (1.8 - 2.8 eV) but certainly not to the dominant red 
band. From their earlier calculations [47,48] described on page 16, it is seen that 
the observed rod band energies and decay times are compatible with the predicted 
emission from silicon clusters of size ~ 10 3 atoms. It may then be possible to think 
of the red emission as originating from the core of the crystallite, whereas the blue 
band is due to the surface chemical complexes. 

Though almost all calculations mentioned above yield a direct, enhanced band 
gap and a non-negligible oscillator strength between the valence band maximum 
and the conduction band minimum, certain worrisome points need to be mentioned. 
Firstly, the band gap underestimation by LDA makes any absolute comparisons with 
PL suspect. Within the first principle calculations themselves, there has been some 
disagreement on key issues such as geometry optimization. Buda and coworkers claim 
relaxation of atoms [41] on the crystallite surface, whereas some other workers [39] 
assume negligible relaxation. The effective masses calculated by two first principle 
approaches [40,42] show the same trend (i.e., they are larger than the bulk effective 
mass) but are significantly different. As far as the radiative lifetimes, and PL energies 
are concerned, one expects at best a semi-quantitative agreement. Ihe presence of 
dangling bonds, excitonic binding and localization energies, non radiative relaxation 
and other radiative centers will defeat any attempt at absolute comparison. Finally, 
three elementary idealizations in most of the electronic structure calculations must 
be pointed out. The material considered is not disordered, perfectly passivated and 
not porous. 
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1.5 More on Luminescence from Silicon - Some 
Surprises. 

The wealth of experimental results in this area is a challenge to the theorist, with 
many workers reporting conflicting observations. For instance, the work of Kanc- 
mitsu, which reports a size independent PL emission [4], is at variance with the other 
experimental works in this field [17,57,58]. Some salient features of PL, about which 
a consensus seems to exist, have been described earlier in Sec. 1.3. The microscopic 
models of Sec. 1.4 have been able to explain - at least qualitatively - the observed 
emission energy range and the red and blue luminescence bands. However, several 
puzzles associated with the detailed nature of the light emitting process still remain. 
In this section, we list some such issues which need to be addressed within a common 
theoretical framework. 

The photoluminescence spectrum 

The microscopic theoretical models have shown that the observed emission ener- 
gies are justifiable within a quantum confinement scheme. The enhanced oscillator 
strength justifies the observed external quantum efficiency to some extent. However, 
the following points still need to be sorted out. 

• The spectra have a large FWHM of 300 - 400 meV. 

• The spectral line shape is asymmetric about the peak on the energy scale, with 
a shoulder on the high energy side. 

• The energies computed for realistic crystallite sizes in electronic structure cal- 
culations are much larger than the observed PL peak energies [39]. To explain 
this, they need to invoke unphysically large values of the exciton binding energy. 

It is evident that to explain the large linewidths, one needs to invoke a size 
distribution of crystallites. We find that such an assumption can also lead to the 

resolution of the other two problems listed here [59]. These calculations are discussed 
in chapter 4. 
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Intensity and decay dynamics of PL 

The intensity as well as the characteristic decay time of luminescence from porous 
silicon exhibit several features which cannot be accounted for within any one of the 
simple models outlined earlier in this chapter. Specifically, 

• The PL intensity exhibits a non-monotonic variation with temperature, with a 
maximum in the range 50 - 150 K [4]. 

• The PL intensity drops exponentially with increasing pressure. 

• The luminescence decay time falls by nearly two orders of magnitude as the 
temperature raised from 10 K to 300 K. 

• The luminescence decay time is seen to decrease exponentially with increasing 
emission energy. 

• The band gap of crystalline silicon is observed to red shift with increasing 
temperature. However, the observed temperature dependence of the PL peak 
is anomalous, with both blue and red shifts being reported by various groups 
[60,61]. 

• The crystalline silicon band gap redshifts with pressure. The PL peak of porous 
silicon may or may not show an initial blue shift depending on the pressure 
transmitting medium employed [62,63]. 

Several workers have attempted to model the luminescence decay dynamics by 
invoking a competition between the radiative and nonradiative process. We show 
that competing dynamics between an Arrhenius type (a exp[-T r /T]) radiative re- 
combination and a Berthelot type (oc exp[r/T s ]) carrier hopping mechanism can 
successfully explain the diverse observations listed above. T r and T B represent char- 
acteristic temperatures associated with these processes. This “unified 1 ’ model for the 
optical behavior of PS is discussed in chapter 5. 


Surface sensitivity of PL 

The PL from porous silicon is seen to be highly sensitive to surface treatments. The 
spectral shape and intensity are seen to vary with aging and surface oxidation. The 
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PL is quenched on annealing, but a brief dip in HF can restore light emission. The 
PL characteristics can be affected by so much as a dip in boiling water. Given these 
vagaries in optical behavior, it becomes necessary to achieve a proper characterization 
of the sample surface. 

The variation of PL with changing surface environments have also been substanti- 
ated by microscopic, theoretical calculations. In an LDA based approach to evaluate 
the energy levels of a planar polysilane, Van de Walle and Xorthrup [64] found the 
gap to be close to 2.75 eV and almost direct. When 50% of hydrogen was substituted 
with OH groups, the gap remained direct and was 1.7 eV. The lowered band gap in 
the OH substituted structure is attributed to the strong electronegativity of the O 
atom. 

Infrared vibrational spectroscopy is a powerful tool to characterize the surface 
properties of porous silicon. The vibrational modes of the surface species can be 
related to the electronegativity of the surrounding environment. In chapter 6 we 
propose a phenomenological model which employs purely theoretical atomic scales 
to predict the stretching mode vibrational frequencies of hydrogen in silicon. This 
model is seen to provide a valuable tool for characterizing the PS surface. 

Absorption chair acteristics of PS 

The excitation of light emission in PS also poses its share of problems. They are: 

• The absorption edge in PS is at least 0.5 eV higher than the observed peak of 
emission. 

• The PL spectral peak seems to blue shift linearly with increasing excitation 
intensity [65]. 

• The PL intensity initially increases linearly on increasing the excitation inten- 
sity, finally saturating to a constant value [66]. 

Besides these, there exists many other issues that stem from a lack of proper 
characterization of the material, leading to apparently contradictory observations. 
The dispute regarding the size dependence of the PL emission spectrum is a case 
in point. We also need to address the various unresolved questions regarding the 
formation process which were outlined in Sec. 1.2. In order :o resolve these issues. 
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one may need to look well beyond the confines of traditional semiconductor theory. 
We now proceed to examine some of these problems in greater detail. 


Key References 

The major references for the material contained in this chapter are listed below: 
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2. Y.Kanemitsu, Physics Reports 263, 1-91 (1995). 


Chapter 2 


Porous Silicon Formation: 

A Generalized Model for 
Aggregation and Etching 

2.1 Introduction 

The study of porous silicon formation spans four decades. Electrochemistry, semicon- 
ductor phenomenology, and more recently, computer simulations have been invoked 
to explain pore formation in silicon. The primary aim of porous silicon formation 
models is to explain, and reproduce, the diverse morphological patterns obtained un- 
der various experimental conditions, as described earlier in Sec. 1.2. We recall that 
a p-type silicon substrate gives rise to a uniform and highly interconnected network 
of random pores, resulting in a spongy structure. As the doping level is increased 
(i.e., p + substrate), the pores become shorter and wider. In n-type silicon, pores are 
formed only under illumination and are relatively linear compared to p-type silicon. 
Illumination is also seen to influence the anodization process in p-type silicon. For 
n + type substrates, the pore patterns formed resemble those grown on p + silicon. 

The various attempts to model porous silicon formation can broadly be catego- 
rized into two, namely (i) Phenomenological theories and (ii) Computer simulations. 
The phenomenological theories generally aim to provide a qualitative understand- 
ing of the pore formation process in terms of competing long-range and short-range 
phenomena. For instance, the Beale model [12] which will be described in greater 
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detail in section 2.5 as well as the next chapter, hypothesizes local modifications to 
the electric field as a result of microscopic irregularities at the acid - semiconductor 
interface. A host of similar models have been proposed, which have been reviewed in 
detail elsewhere [2]. 

There exists another class of models which attempt to obtain a quantitative for- 
malism by examining the stability conditions at the semiconductor surface [67.68]. 
This approach is patterned on the lines of the classical Mullins-Sekerka stability anal- 
ysis [69.70]. However, because of the inherent mathematical complexity of the system, 
such attempts have enjoyed only a limited degree of success. 

Another approach, which we wish to employ in tackling the problem of porous 
silicon formation takes recourse to computer simulations which have been extensively 
employed in the past to model aggregation and etching phenomena [71]. In the next 
section (Sec. 2. 2), we briefly survey existing studies of growth and aggregation phe- 
nomena. In section 2.3 we propose a generalized computational scheme termed the 
drift-diffusion model to simulate the classes of morphologies seen in porous silicon as 
well as in other processes such as electrochemical deposition (ECD), viscous fingering 
and amorphous annealing [72]. The results of these simulations are outlined in sec- 
tion 2.4. Section 2.5 constitutes a discussion where we examine the extension of the 
drift-diffusion model to the case of porous silicon formation. The subsequent chapter 
describes various phenomenological models proposed earlier to explain porous silicon 
formation and extends our computational model to simulate the specific instance of 

silicon anodization in HF. 


2.2 Existing Models of Aggregation and Etching 

The studies of various morphologies obtained in aggregation phenomena such as elec- 
trochemical deposition, viscous fingering, bacterial colony growth etc., have attracted 
a great deal of attention since the early eighties [73], In electrochemical deposition 
experiments, two parameters, the electrolyte concentration (C) and the applied po- 
tential (V) are tuned to obtain (1) dendritic structures, both thick and needle like (or 
stringv); (ii) Dense branching morpholog.es (DBM) (also known in literature as Eden 
trees [74)). which are homogeneous and (hi) randomly ramified self similar structures 
as in diffusion limited aggregation(DLA) [75,76). In a Hele-Shaw cell, similar transi- 



24 


A Generalized Model for Aggregation and Etching 


tions are observed when the pressure and surface tension are varied [ / ' ]. In silicon, 
anodic etching gives rise to differing pore morphologies depending on the anodization 
potential and the substrate doping level [11,12], as described earlier in Sec. 1.2. 

A first principles explanation of such phenomena will have to encompass (i) the 
diffusive field, (ii) the Laplace field, (iii) convective processes and (iv) surface ten- 
sion. curvature and underlying anisotropy effects. Hence, even numerical solutions 
may prove difficult. Simple formal approaches have been attempted for obtaining 
stability conditions at the interface [69,78.79]. It has been hypothesized that the 
emergence of different characterizing length scales is the result of the interplay be- 
tween the Laplace and diffusion fields governing the growth process in ECD [79;. 
Early phenomenological models for PS formation (to be described in detail in the 
next chapter), also invoke two scales in the form of a depletion layer width at the PS 
- HF interface and a barrier lowering at the pore tips depending on the curvature. 

It is in this context that the statistical models which can simulate physical pro- 
cesses taking place in a large system become relevant. The modeling of complex 
growth phenomena has become increasingly possible with the advent of fast comput- 
ers. Simple statistical algorithms have been able to reproduce the patterns obtained 
in varying phenomena like amorphous annealing and dielectric breakdown [77]. Some 
of the simplest "growth ” laws proposed in literature [71] are outlined below. In these 
models, the central seed on a lattice can grow by 

1. Random occupation of a surface site (i.e., nearest neighbor site) of the seed. 
This is known as the Eden model [80], which results in compact clusters. 

2. The ballistic transport (i.e., movement in a straight line) of a particle generated 
at a randomly chosen site on the lattice boundary to the seed surface. This is 
called ballistic deposition. 

3. The diffusive transport (i.e.. a Brownian random walk) of a particle generated 
at a randomly chosen site on the lattice boundary to the seed surface. This is 
labeled diffusion limited aggregation [81]. 

Ever since the diffusion limited aggregation (DLA) model was suggested in 19S1. 
a great deal of attenrion has been devoted to developing variations of the algorithm. 
Patterns generated by the DLA were found to resemble the growth geometries ob- 
tained in electrochemical deposition, dielectric breakdown and viscous fingering. The 
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present work constitutes a generalization of the DLA. We first outline the algorithm 
for the diffusion limited aggregation process: 

1. Start with a central seed on a lattice. 

2. A particle is launched on a random walk from a randomly chosen site on the 
lattice boundary (the step size of the walk should be much smaller than the 
lattice size). 

3. \\ hen the particle comes into contact with the seed (i.e., walks into a nearest 
neighbor site of the seed) the random walk is stopped and the particle is added 
to the seed. 

4. Steps 2 and 3 are repeated 

This algorithm employs no controlling parameters. The resulting pattern can be 
characterized by the fractal dimension df [71.32] defined as 

M = R d} 

where M is the mass of the aggregate (the total number of particles) and R the radius 
of the cluster. For the DLA pattern this was found to be approximately 1.7 [73]. This 
is in good agreement with the fractal dimensions calculated for zinc leaves obtained 
in electrochemical deposition experiments at very large dilutions. 

Generalizations of the DLA have been attempted by the introduction of addi- 
tional parameters into the model. Surface tension effects have been modeled by the 
incorporation of a sticking probability at the aggregate surface [83], In this model, 
a particle is allowed to stick to the cluster only after n contacts. Here, n is a scale 
characterizing surface tension effects. With increasing value of n, a transition was ob- 
served from the low density fractal clusters to regular patterns exhibiting a compact 
structure. 

Yet another scale to represent the particle concentration can be introduced by 
employing a multi-walker algorithm [84] wherein several walkers are released simulta- 
neously from the lattice boundary, instead of a single particle as in DLA. The number 
of particles at the lattice boundary is maintained constant throughout the simulation. 
This is indicative of a const ant particle concentration at the boundary ■ Recentk . a 
multi-walker DLA model has been proposed [85] wherein two parameters, namely 
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the particle concentration p and the width of a '‘migrational envelope l are tuned 
to obtain dense branching morphologies as well as DLA patterns. In this simulation, 
the particles undergo diffusive motion outside the migrational envelope. On reaching 
the envelope, they migrate to the nearest point on the aggregate surface. The migra- 
tional envelope is continuously modified to maintain a minimum distance l from the 
aggregate contour. By varying the values of the parameters l and p one can obtain a 
transition from DLA to DBM. 

It will be a desirable goal to arrive at a model which can reproduce the various 
morphologies obtained under different growth conditions as well as the transition 
from one morphological class to another. From the discussion in the preceding para- 
graphs. it is clear that the morphological variations result from the interplay of two 
or more parameters [S6.87]. For the sake of simplicity, a minimal set of parameters 
should be employed. In the next section, we propose a drift-diffusion model which 
employs the competition between two antithetical parameters to simulate the various 
morphological classes described in the beginning of this section. 


2.3 The Drift-Diffusion Model 

In the classical DLA model [81] described in Sec. 2.2, a particle starts its random 
walk at an infinitely large distance from a seed or cluster. The random diffusive mo- 
tion is terminated the moment it comes into contact with the aggregate cluster. It is 
conceivable that very close to the aggregate surface, particle movement is no longer 
controlled by the macroscopic diffusion field, but by microscopic, localized phenom- 
ena. To model such processes, additional parameters will have ro be introduced into 
the simulation. We attempt to develop an algorithm which models pattern formation 
as a result of the interplay of a macroscopic field and localized surface phenomena. 

The main control parameters in this model are: (i) a depletion layer width All' 
which controls the diffusion and (ii) a drift length l governing the held driven processes 
in the proximity of the aggregate. The design of the algorithm is outlined below: 

(i) A particle is launched from a randomly chosen lattice site beyond the depletion 
layer boundary which is at a distance All from the surface of the aggregate. This is 
illustrated in Fig. 2.1. We begin with a circular depletion zone or radius AIF around 
a central seed. As the aggregate grows, the depletion zone boundary is modified to 




Figure 2.1: The drift-diffusion model. A particle is released beyond the depletion 
layer boundary (the closed curve) which dynamically follows the aggregate (solid 
circle cluster) contour at a distance AIT. When the particle wanders to within a 
radius l (see dotted arrow) of the aggregate surface, it is “field driven” to become a 
part of the cluster. Note that the figure depicts the case l > AIT. 

follow the aggregate contour at a distance AIT as shown in the figure. 

(ii) The particle executes a random walk as in ordinary off-lattice DLA algorithms 
[88,89]. A circular contour (of small radius compared to the lattice size) is drawn 
around the the particle. The particle is then moved at random to any of the lattice 
sites which fall on this contour. 

(iii) The moment a part of the aggregate surface is encountered within a radius l 
of the particle location, the random walk is terminated. The particle is then moved 
to the surface (in a sense, “field driven") and becomes part of the aggregate. 

(iv) Steps (i) - (iii) are repeated. 
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The simulations were stopped when the aggregate had grown to a diameter of 
around 400 units, one unit being the diameter of a single particle. Approximately 
10 4 - 10 5 particles had to be launched for various cases. 

The two control parameters employed in this simulation have also been suggested, 
albeit separately, elsewhere in literature. To model porous silicon formation, a “finite 
diffusion length” similar to our AIT has been proposed [90]. The drift length l in our 
simulation is in a way similar to the width of the “migrational envelope” of Erlebacher 
et al. [85,91], Our algorithm interpolates between existing growth models in various 
limiting cases. In the limit / = 1 and AIT — ¥ oo, the algorithm is identical to the 
DLA. AW = 0 and l = 1 approximates the Eden limit described in Sec. 2.2. As I 
becomes larger than All' (the case depicted in Fig. 2.1), a fraction of the generated 
particles are directly transported to the surface without undergoing diffusion. This 
is akin to ballistic deposition (see Sec. 2.2). When l » AH', the ballistic process 
dominates the diffusive process. 


2.4 Results 

Fig. 2.2 depicts the variety of morphologies obtained with varying AIT and l. For 
small 1/ AW , we obtain DLA like patterns (Fig. 2.2a). A simple mass-radius scaling 
calculation for this pattern yielded a fractal dimension of 1.65 ± 0.05. Stringy struc- 
tures similar to those reported in electrochemical deposition at high voltages [75] are 
seen in the large AIT limit with l JAW ~ 1 (Fig. 2.2b). The mass-radius scaling 
exponent fluctuates around ~ 1.45 for this pattern. The stringy morphology observed 
in some experiments is essentially one dimensional in character [76]. This is repro- 
duced for very large AIT(= l ) > 25. As AW{= l ) is increased, the scaling exponenr 
approaches unity (see Fig. 2.3). On the other hand, for small AIT and l/AW ~ 1. 
we obtain short dense branches exhibiting a relatively smooth front which remains 
nearly circular throughout the growth period (Fig. 2.2c). This has been identified in 
experimental literature as homogeneous or tip splitting patterns or dense branching 
morphologies [75.76,78,92]. The mass radius scaling for these structures resulted in 
a fractal dimension ~ 2. For l/AW > 3, thick dendritic growths with side branching 
was obtained (Fig. 2.2d). For large AIT(> 5), the patterns become too inhomo- 
geneous to show any well defined scaling [75], However, for small AW, the model 




Figure 2.2: The simulated growth patterns obtained by varying the depletion width 
‘AH' ’ and drift length T. (a) Patterns resembling DLA (AW =10 ,1 = 2), (b) 
stringy structures (AW = 10, l = 10), (c) homogeneous, dense branching morpholo- 
gies with a nearly circular growth front (AW = 2, l = 2) and (d) thick dendritic 
growth with side branches (AW = 8 and l = 24). 

approaches the Eden limit and the resultant patterns are somewhat homogeneous. 
The occurrence of dendritic growth (see Fig. 2.2d) in such a simulation is sur- 
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Figure 2.3: The fractal dimension df plotted against AH’ (in lattice units) with / 
equal to AFT. df is seen to approach unity in the large AW(= l) limit. 

prising, since it is well understood that anisotropy is required in the interfacial dy- 
namics for side branching to occur [77]. The present algorithm admits an underlying 
anisotropy in the way a particular surface site is chosen, when more than one point 
on the aggregate surface falls within a radius l of the random walker. This selection 
was done in three ways, (i) The nearest site was chosen, (ii) A site was chosen at 
random, (iii) The site nearest to the radial line connecting the random walker and 
the central seed was chosen. The morphologies remained similar for the case l < AH’. 
The three methods yielded different results in the limit l » A W . Dendritic growth 
was observed only in case (iii), which corresponds to the imposition of a preferred 
radial field. In the case of electrodeposition, this could be interpreted as the applied 
radial electrostatic field. 

Figures 2.3 and 2.4 show the variation in the mass-radius scaling exponent dt 
with varying values of AIT and l. As AIT is kept equal to l and varied. d f decreases 
steadily with increasing AIT(= l) and approaches unity. This is depicted in Fig. 2.3 
and represents the emergence of string}' structures as shown in Fig. 2.2. The values 
obtained for various AIT when / is kept constant (/ = 2) are plotted in Fig. 2.4. Tire 
dimension df is seen to vary from 2.00 to 1.65 ± 0.05 as AIT goes from very -mail 
(~ 1) to large (> 15) values. This illustrates a smooth transition between Eden type 
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Figure 2.4: The fractal dimension df plotted against AW (in lattice units) with 1 = 2. 
For small AW, df ~ 2 (see Fig. 2.2c), which is indicative of compact clusters. In the 
large AW' limit, df stabilizes at 1.65 ± 0.05 (see Fig. 2.2a). 

growth and DLA. 

Fig. 2.5 codifies our observations into a “phase" plot l . The plot was constructed 
on the basis of approximately a hundred patterns, grown to aggregate diameters of 
around 400 units. The transitions from one phase to another being continuous, the 
boundaries are not rigidly defined. For example, in the / AW case, the dendritic 
patterns observed at large AH' become denser and more homogeneous as A IF ap- 
proaches zero. The phases depicted have been reported earlier in electrochemical 
deposition experiments [75,76] and related processes [77]. 

2.5 Discussion 

One can explore the phenomenological relevance of the parameters AW and Z. The 
existence of competing processes in electrochemical deposition and other phenomena 
is well established. The phase diagrams reported by Sawada et al. [75] and Grier et al. 
[76] represent various morphological classes resulting from varying the concentration 
and applied voltage. Ass i gnin g a direct correlation between our parameters and the 


x the term “ phas es” is taken to mean classes of aggregates as reported earlier in literature 
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Figure 2.5: “Phase” plot of patterns observed in the simulations. The phases depicted 
transform into one another when the parameters AW and l are varied. The bound- 
aries must not be treated as rigidly defining a transition and are hence depicted by 
shaded regions. Both AW and / are expressed in terms of lattice units. 


experimental parameters based on a comparison between the phase diagrams of the 
simulation and experiment is difficult. In all probability, AW’ and l are functions of 
concentration, applied voltage and the interface effects. We shall further explore the 
significance of these parameters in the anodization of silicon in the next chapter. 

In the diffusion field versus drift field dynamics governing growth in our simula- 
tion, AW is a control parameter for the macroscopic diffusion field, whereas / models 
microscopic surface effects. An example of a phenomenological model employing sim- 
ilar arguments is the Beale model [12] (this is discussed in greater detail in the next 
chapter) for porous silicon formation which hypothesizes the existence of a semicon- 
ductor depletion layer at the Si-HF interface. The width of this depletion layer varies 
as 


AH' oc 


n 0 


(2.1) 


where V& is the anodization potential, V B1 is a constant built in voltage and no is 
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the substrate doping level. On the other hand, the parameter / can be related to 
the barrier lowering (A <f> s ) due to microscopic surface irregularities. This lowering, 
Ao s \fE~s- E s being the enhanced local electric field. For planar interfaces, E s oc 
Os/All where o s is the overall barrier height. Due to surface irregularities, this 
height is enhanced locally as 


E s ca 


ds 

(Kw - /) 


(l < AH 7 ) 


( 2 . 2 ) 


As l —7 AM 7 , E s is very large at the tips, leading to stringy patterns. In a recent 
work [93], it has been shown that A W and l can be correlated to the depletion layer 
width and barrier lowering given in equations (2.1) and (2.2). 

Several studies of two dimensional growth on a one dimensional substrate have 
been reported in the past. Matsushita et al. have grown zinc “trees" on a linear car- 
bon cathode [94]. Such processes have also been modeled using a DLA like algorithm 
[95,96]. There exists a diffusion limited annihilation model [97] where a diffusing 
particle etches away a site on a flat surface. The etching process in porous silicon 
may also be viewed similarly. One approach [90] has been to introduce a “diffusion 
length” which is defined on similar lines as A W in the drift-diffusion model. The 
muti-walker model of Erlebacher and coworkers has also been extended to model sil- 
icon dissolution [91]. In the next chapter, we attempt to extend the drift-diffusion 
model to porous silicon formation by using a linear substrate which is represented by 
a side of the lattice. The deposition is no longer radial, but “quasi-one dimensional'’. 

In a sense, the two scales AH 7 and / are antithetical, one separating the particle 
from the aggregate, the other driving it towards the aggregate. Various growth models 
(DLA. ballistic deposition, Eden) can be obtained as limiting cases of this algorithm. 
We stress that the two scales do not contain any a priori bias towards any specific 
morphological structure. Nevertheless, a variety of distinct morphological structures 
observed in experimental growth and dissolution is obtained in our simulation. 


Key Reference 

The major reference for the material contained in this chapter is: 

• George C. John and Vijay A. Singh, Phys. Rev. E 53, 3920-3924 (1996). 



Chapter 3 


Porous Silicon Formation: 

The Diffusion Induced Nucleation 
Model 


3.1 Introduction 

As mentioned earlier, porous silicon is formed during the anodization of silicon in HF 
acid. The formation of diverse pore structures rather than a uniform ‘piling off of 
the silicon layer presents an intriguing problem. One needs to explain the preferential 
pore formation in comparison to a uniform etch, as well as the observed transition 
in morphologies with varying experimental parameters. The early phenomenological 
theories sought to accomplish this in a qualitative fashion, by means of arguments 
premised on a particular model of the semiconductor - electrolyte interface. 

The majority of the growth models proposed to explain pore growth in silicon do 
not concern themselves with the exact dissolution chemistry, but invoke physical argu- 
ments to explain the influence of anodizing conditions on porous silicon morphologies. 
This is because of the fact that the detailed dissolution chemistry of silicon is com- 
plex and not well understood. Several dissolution mechanisms have been proposed, 
an excellent review of which is found in the article by Smith and Collins [11], 

In the Si - HF system of interest to us, the etching process can be influenced 
by a variety of factors such as (i) the electrochemistry of dissolution, iii) the semi- 
conducting properties of silicon, (iii) transport of carriers in the semiconductor, (iv) 
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transport of ions in the electrolyte and (v) surface phenomena in the form of irregu- 
larities on the silicon surface or the surface tension of the liquid. The Si-HF system 
is forward biased when a p-type Si substrate is used as the anode. In n-type silicon, 
pore formation occurs only when some extraneous hole generation processes such as 
intense illumination or high electric fields are invoked. This seems to suggest that 
the presence of holes is essential for porous silicon formation. 

A substantial amount of work has been done on diffusion limited models in the last 
decade and a half [71.88.98]. The recent models of porous silicon formation have been 
based on detailed computer simulations of diffusive carrier transport within silicon 
and are briefly reviewed in section 3.3. These models are premised on similar lines 
as the diffusion limited aggregation (DLA) which has been described in the previous 
chapter. 

The drift-diffusion model outlined in chapter 2 is successful in reproducing the 
morphological classes observed in diverse growth phenomena. These patterns bear 
considerable resemblance to the pore morphologies in silicon. For instance, the dense 
branching morphologies seen in Fig. 2.2c are similar to the highly interconnected 
network of pores obtained with p-type substrates. The stringy structures (Fig. 2.2b) 
are characteristic of n-type silicon substrates, where pores grow in a relatively linear 
fashion. In degenerately doped silicon, thicker pores are obtained which resemble the 
AU’ < / case of the drift-diffusion model (Fig. 2.2d). However, one would like to go 
beyond a picture-to-picture comparison and evolve a model which should be capable 
of explaining other features of the experimental growth process as well. Some of these 
are: 


• Rate of growth: The porous layer is reported to exhibit a constant growth rate, 
with the velocity of the pore front remaining independent of time [99.100]. Fur- 
ther, one would also like to explore the dependence of the growth rate on the 
experimental parameters such as the anodization potential, electrolyte concen- 
tration etc. 

• Porosity profile: PS exhibits a constant porosity vs. depth profile with a sharp 
planar film front [101]. Considering the strong dependence of visible lumines- 
cence on porosity, it would be desirable to have a model which can predict the 
variation in porosity with experimental conditions. 
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• “Piping”: At high anodic potentials, the linear pores formed on n-type sub- 
strates widen and become pipe-like [11]. This phenomenon is known as “piping 1 
and is not well understood. 

• Electropolishing: Silicon, both p and n-type, are seen to be uniformly etched 
away at high voltages, instead of forming porous silicon. This phenomenon, 
termed “electropolishing” is also seen to occur at low electrolyte concentration 
(~ 10%). 

In this chapter, we extend the drift diffusion model outlined in chapter 2 to the Si- 
HF system, based on a careful examination of the earlier phenomenological theories. 
We attempt to pick out the key processes and incorporate them into the algorithm. 
Two simple processes are identified as contributing to the primary dissolution mech- 
anism, namely (i) the diffusion of holes from the bulk towards the Si-HF interface, 
and (ii) further propagation of a pore from the site of dissolution with a finite ter- 
mination probability p t (a nucleation process). The model incorporating these two 
processes, termed the diffusion induced nucleation (DIN) model is described in sec- 
tion 3.4. The various parameters and their numerical values are justified on the basis 
of physical arguments. Section 3.5 illustrates the different morphologies obtained 
using the algorithm and compares them with experimentally grown porous silicon 
structures. Section 3.6 constitutes a discussion, where experimental observations like 
the rate of growth of the film, variations in the pore diameters and electropolishing 
are compared with the results of the simulation. Some experiments for further study 
and elucidation are suggested. 


3.2 Phenomenological Theories 

Most of the early models of porous silicon formation are based on the I-V character- 
istic of the dissolution process. The general I-V behavior for the Silicon-HF system 
[11] is briefly summarized below (see Fig. 3.1): 

(i) Porous silicon formation occurs in zone A, when the current I is less than 
the value at the first peak l PS . Zone B is a transition region to the electropolishing 
regime depicted as zone C. 
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Potential (arb. units) 


Figure 3.1: Typical I-Y characteristic (schematic) for anodization of p-tvpe silicon in 
HF. The three zones represent the pore formation regime, a transition region and the 
electropolishing regime. 

(ii) Since Si is the anode, the system is forward biased for p-type silicon, whereas 
for n-type silicon it is reverse biased. Rectifying behavior is observed in both cases. 

(iii) Under illumination, n-type and p-type Si exhibit similar I-V characteristics. 
For the heavily doped case (i.e.. n + and p + ) the characteristics are similar even under 
dark conditions. 

(iv) In the porous silicon formation region (zone A), the curve is exponential. 

Several models for porous silicon formation proposed in the last decade try to 

correlate these facts to the complex pore morphologies observed in silicon. The most 
prominent of these is perhaps the Beale model [12]. This model assumes that a current 
flow is necessary for porous silicon formation, with Si as the anode. In p-t\ pe silicon 
this is easily possible as the Si-HF system will be forward biased. But for n-type 
silicon which will be under a reverse bias, current flow requires intense illumination 
or high voltages which can cause breakdown. The current flow in silicon is localized 
as a result of surface inhomogeneities. This leads to a non-uniform dissolution of the 

surface, initiating pores. 

The semiconductor - electrolyte system is modeled as a Schottky interface 
(Fig. 3.2) by Beale and coworkers [12], and some others [102]. The etched surface 
supports several surface states which can act as traps for carriers. The high density 
of surface states results in the pinning of the Fermi level in the mid-gap. This acts 
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Figure 3.2: (a) Representation of the semiconductor - electrolyte system according to 
the Beale model [12]. Energy band diagrams are drawn assuming a Schottky barrier 
for (b) p-type Si, (c) n-type Si. The solid lines represent the conduction (E c ) and 
valence (E v ) bands without any anodizing potential. The dashed lines show barrier 
raising (lowering) after a potential is applied. The Fermi level (dotted line) is assumed 

to be pinned to the mid-gap at the interface because of the high density of surface 
defects. 

as a barrier for current flow and can explain the observed rectifying behavior. 

Carrier transport across the barrier can then occur through thermionic emission 
or by tunneling. In the case of lightly doped Si. the tunneling probability is small 
and hence the current flow is determined mainly by thermionic emission, which is 
controlled by the barrier height. The barrier height depends on the applied bias as 
well as the local electric field distribution. Charges near the barrier produce image 
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HF 



SILICON 

Figure 3.3: The ‘'pinch off" effect due to the overlap of depletion layers (of width 
AIF) associated with adjacent pores. The pore diameter is d p , the pore tip radius is 
r. Note the reduction in the depletion layer width around the pore tip (AIF r < AFF). 

charges in the electrolyte, thereby enhancing the local electric field. The barrier 
lowering is hence, 


A Oi = q(E s /47re s e 0 q) 1/2 

where E„ is the electric field at the interface, e s is the relative permittivity of silicon, 
€ () is the permittivity of vacuum, and q is the carrier charge. For a planar barrier 
with a depletion layer width AH’, E s is given by 

E s = 2o s /AW 

where <p., is the barrier height. Around a surface irregularity which can be modeled 
as a hemispherical tip of radius r, (AIT > r) this becomes (Fig. 3.3) 

E s = <Ps/r 

Thus we see that there is a considerable field enhancement at regions of high curva- 
ture. This leads to enhanced silicon dissolution as a result of increased local current. 

For heavily doped silicon, the Si-HF system conducts in both directions. The 
behavior is non-ohmic. This is best explained by the theory that carrier transport is 
then governed by carrier tunneling and not thermionic emission over the barrier. In 
this case, the current is not controlled by the barrier height, but' by the width (be. 
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the width of the semiconductor depletion layer) which can be written as 


AW'=(2 


( 3 - 1 ) 


where N c is the carrier concentration. According to Eq. (3.1), the depletion layer 
width will be small for degenerately doped silicon. As in the case of lightly doped 
silicon, field enhancement around the pore tips can locally reduce the depletion layer 
width. Modeling the tip as a hemisphere of radius r, the width of the hemispherical 
barrier AW r is given by 


1 - 



2 /A W r \ 3 AU-' 
3VAirJ r 


( 3 . 2 ) 


The depletion layer width thus acts as a regulating parameter for the pore growth. 
When depletion layers associated with adjacent pores overlap (Fig. 3.3), there is a 
pinch off effect, leading to depleted pore walls. This results in increased resistance 
of the porous layer, and the current preferentially flows down the electrolyte. The 
porous film is thus stable against further dissolution. This pinch off can also lead 
to a columnar pore structure which is observed for n-type substrates. However, 
perturbations due to defects and grain boundaries can redirect the pore growth. 

If a pore outgrows its neighbors, tunneling can occur over an increased area of 
the tip, leading to a larger tip radius. This, in accordance with Eq. (3.2), results in 
a larger depletion layer, which serves to limit carrier tunneling resulting in a slowing 
down of the rate of growth. Such a regulating mechanism is considered responsible 
for the more or less planar porous film front. 

The Beale model, described so fax in this section, is thus constructed in terms of 
familiar semiconductor device physics, not accounting for the electrochemical reac- 
tions at the surface. In an attempt to explain pore growth in n-type silicon, Zhang 
[102] also takes recourse to a depletion layer as the controlling parameter. Carrier 
transport is assumed to take place through tunneling across the barrier. If a potential 
much higher than kT is applied, band bending occurs, which reduces the depletion 
layer width. The model is much the same as that of the Beale model and assumes 
that pore tips (irregularities) cause band bending leading to an enhanced current. A 
maximum limit is set on the pore wall thickness, which is twice the depletion layer 
width (Fig. 3.3) This is a consequence of the pinch off effect due to the depletion layer 
overlap. The main difference from the Beale model is that Zhang invokes carrier tun- 
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neling even in lightly doped silicon, whereas the Beale model assumes thermionic 
emission as the major carrier transport mechanism. 

These models are consistent with the following experimental observations: In- 
crease in dopant concentration leads to shorter pores with less spacing between them 
compared to the lightly doped case. The semiconductor depletion layer width A W 
varies inversely with \d\ c as indicated in Eq. (3.1). Besides, as the HF concentration 
is increased, the pore becomes narrower. This is expected, as the current at the pore 
tips increase. Also, p-tvpe silicon has a much smaller depletion layer compared to 
n-tvpe which can give rise to narrow and uniform pores. However, p-type silicon 
produces an interconnected network of pores which such models cannot satisfactorily 
explain. 

Both the above models share the idea of a depletion layer, which is substantiated 
further by others [103]. An alternate theory exists [104,105], which assumes the 
formation of a passivating layer on the Si-HF interface. In this case, the etching 
is limited by the local dissolution of the layer. The pore formation dynamics is 
controlled by the competition between the growth and the dissolution of the passive 
layer. Analytical and numerical calculations on this model have been successful 
in reproducing the observed morphological variations with changes in the applied 
voltage V and the HF concentration. However, the pore wall thickness d is found 
to 1.x* strongly dependent on the doping concentration. This cannot be satisfactorily 
explained by this model. 

An appealing model to explain the pore wall thickness is the quantum confinement 
model, proposed by Lehmann and Gosele [103.106-108]. The basic idea of this model 
is that as a result of carrier confinement in the nanostructures formed during silicon 
dissolution, there is a band gap enhancement. Only those carriers with energies high j 
enough to overcome this barrier can now be present in the crystallite. This results 
in hole depletion in the porous structure. As holes are essential for dissolution, no 
further etching occurs. 

We have seen that a simple depletion type model sets an upper limit on the wall j 

thickness of the porous silicon film. Similarly, the quantum confinement effects are j 

responsible for setting a lower bound on the crystallite size. According to Lehmann [ 
and Foil [99], the macropores formed in n-type silicon are covered with microcrys- 
talliies on the pore sides. The size of these microcrystallites could be limited by 
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quantum confinement. However, there exists other reports which indicate that this 
picture of n-type porous silicon is not necessarily true [11]. Reports of a possible 
fractal structure also exist [11], though this is contested by Lehmann and coworkers. 

It is evident that the quantum confinement model cannot be regarded as a primary 
mechanism for porous silicon formation. At best, it is a secondary effect putting a 
lower bound on the crystallite size. For the onset of pore formation, a depletion layer 
seems to be necessary, the width of which controls many large scale features of the 
final morphology. However, it could.be argued that in p-type silicon, the space charge 
region has a very small width, and hence the quantum confinement effects can play 
a major role. 

The semiconductor depletion layer is functionally equivalent to a passivating layer, 
which can prevent dissolution of the pore walls. Thus we see that the semiconductor 
depletion layer, the passivating oxide layer and the nanocrystallites with an enhanced 
gap serve the same function, namely create a carrier depletion at the surface. In short, 
one can view this depletion layer width as a control parameter which serves to “push" 
the pores apart. This is reminiscent of the role of the depletion parameter AW in the 
drift-diffusion model of chapter 2, which serves to keep the diffusing particle away 
from the aggregate. The idea of such a depletion layer has been gainfully employed 
by several computer simulation studies of porous silicon formation, which we shall 
now examine. 


3.3 Computer Simulations of Porous Silicon For- 
mation 

Diverse patterns are observed in natural processes such as dielectric breakdown, vis- 
cous fingering and electrochemical dissolution [88] which have been successfully mod- 
eled using the diffusion limited aggregation (DLA) model described in the previous 
chapter. Since the pore morphologies in porous silicon appear to resemble the DLA 
patterns, several diffusion limited mechanisms have been proposed to model porous 
silicon formation. 

In a computer simulation, Smith, Chuang and Collins [90] assume that the diffu- 
sion of an ‘electro-active species ; (holes) from the bulk to the Si-HF interface is the 
controlling factor in porous silicon growth. The bulk silicon is represented by a two 
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dimensional rectangular lattice of sides M x .V. Periodic boundary conditions are 
implemented in the lateral directions, which means that a particle 'walking off' one 
side of the lattice reappears at the opposite side. This algorithm, termed the "finite 
length diffusion model" is briefly described below: 

(i A hole is released from an isoconcentration surface within the lattice. This is 
similar to the depletion layer boundary of the drift diffusion model (see Fig. 2.1 on 
page 27). The isoconcentration surface is constituted by the locus of all points at a 
distance All’ from the Si-HF interface (in this case, a side of the square lattice). 

(ii : The hole executes an on-site random walk on the lattice, till it contacts a 
surface site. Then that surface site is dissolved with the annihilation of the hole. 

(iii . The isoconcentration surface is modified to follow the contour of the freshly 
dissolved surface. 

(ivi Another particle is released from the iso concentration surface. Steps (ii) and 
(iii) are repeated. 

The pore geometries resulting from this simulation bear some resemblance to 
TEM micrographs of porous silicon. The morphological features were controlled by 
the parameter AW, which Smith and coworkers assume to be a diffusion length. The 
pore spacing is found to depend on AW because of the following reason. When two 
pores are separated by a distance less than or equal to 2AW, the isoconcentration 
surfaces associated with each overlap. Now the holes released from this surface have a 
much lesser probability of reaching a site deep in the resulting wall like structure. This 
phenomenon is similar to the pinch-off effect caused by the depletion layer overlap in 
the Beale model (Fig. 3.3). Smith and Collins claim [11] that the diffusion length is 
in a sense related to the width of the semiconductor space charge region. 

One can also perceive an in-built mechanism wherein the hole diffusion is prefer- 
entiallv directed towards the pore tips. The sites at the tips sample a greater extent 
of the isoconcentration contour (depletion layer boundary) than a site on the pore 
sides, because of reduced interference from other surface sites. Therefore, the pore 
tips grow much faster, resulting in a uniform density of the porous film. 

This model does not yield the high porosity values observed in the case of light 
emitting porous silicon. The maximum reported values in the simulation are around 
40%. Besides, the pores are never interconnected as in the case of p-type porous 
silicon. The model does not produce wide pores which are seen in the ‘piping phe- 
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nomena in n-type silicon [11] described earlier in Sec. 3.1. 

This model is functionally similar to the Beale model, though it a\ oids man\ of 
its assumptions such as pinning of the Fermi level at mid-gap and the oversimplified 
model of a Schottky barrier at the semiconductor electrolyte interface. For a compar- 
ison between the two models, we refer the reader to the review article of Smith and 
Collins [11]. The microscopic features are better modeled by the diffusion limited 
scheme. This model also predicts a fractal to non-fractal transition at a depth AH 
from the surface of the porous film [109]. Experimentally it is known that porous 
silicon consists of a nanoporous layer at the top and a macroporous layer of uni- 
form porosity below [101]. However, the observed thickness of the nanoporous layer 

1 pm) is too large to be explained by a simple relation to the diffusion length which 
is of the order of the nanocrvstallite size as seen from the Beale model (Fig. 3.3). 

Further, the finite length diffusion model is extremely oversimplified, assuming 
only one control parameter AW for the entire process. Therefore, the claim that it 
is an exact equivalent of the Beale model [11] is invalid, as more parameters would 
be required to represent the effects of the dopant type, concentration and the electric 
field. Besides, a simple one to one correlation between parameters of both models as 
established by Smith and Collins does not seem convincing. 

A similar model employing a control parameter p, which represented the concen- 
tration was proposed in 1984 by Voss and Tomkiewicz [110] to simulate electrodepo- 
sition in metals. There the simultaneous diffusion of a large number of particles was 
employed, with the population of such walkers being controlled by the parameter p. 
The same model is applicable for porous silicon formation if we regard the aggregation 
of particles in an opposite sense, i.e., the removal of sites. The morphologies are sim- 
ilar to those obtained by Smith and coworkers, though a better variation is achieved 
by employing a second control parameter in the form of a ‘sticking probability' at the 
surface. A random walker, on contacting a surface site may stick to it with a finite 
probability or continue to walk. At low sticking probability, mossy deposits which 
resemble the transition region to the electropolishing regime in porous silicon was 
observed. When the sticking probability was maximum, dendritic growth resulted. 
Thus a transition between mossy deposits and fully developed dendritic growth was 
recorded. Note that these parameters are similar to those discussed earlier in Sec. 2.2 
for the case of radial growth simulations. 
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Erlebacher and coworkers [85.91] introduce a two scale model which employs a 
multi-walker algorithm as described above. They attempt to model the effect of 
the high electric field focused at the pore tips by a ‘migrational envelope' which is 
constructed as a locus of points equidistant from the pore surface. This follows the 
contour of the growth front at a fixed screening distance l. A particle, on contact- 
ing the migrational envelope during the course of its random walk is preferentially 
transported (field driven) to the nearest pore tip, as in the case of the drift-diffusion 
model. The migrational envelope is claimed to be analogous to the depletion layer in 
n-type silicon. The other control parameter is the concentration p of particles, which 
determines the number of holes in the system. The concentration is maintained con- 
stant at the edges of the lattice, far from the interface. Particles were not allowed to 
overlap during the course of their on-site random walk. 

Many of the characteristic morphologies of porous silicon are obtained in simu- 
lations with low concentration and large screening lengths. Highly directional pores 
are seen, which do not grow into each other. This is characteristic of porous silicon 
formed from n-type substrates. The simulation would correspond to illuminating the 
sample from the back where the carrier concentration due to hole generation will be 
maintained constant. As the concentration increases, the spacing between the pores 
decreases. The density profile with depth remains constant. When / is set to zero, 
the simulation mimics that of Voss and Tomkiewicz with sticking probability one. 
The porosity obtained in each simulation was approximately equal to p. If p is taken 
to represent tin 1 hole concentration, such a result is clearb unphvsical. 

We can also study two dimensional growth on a linear substrate using the drift- 
diffusion model described in chapter 2.This is achieved by using a linear substrate 
which is represented by a side of the lattice. The deposition is no longer radial, 
but “quasi-one dimensional.” Fig. 3.4 depicts the dependence of the mean aggregate 
density on the drift length l. Deposits were grown on a linear substrate of length 
300 units for a given AH' . A cross-over behavior is discernible, v ith a minima at 
AID = l. This can be understood in terms of our phase plot (Fig. 2.5 on page 32). 
The minima corresponds to the string}- region, where AW ~ l On either side of 
the stringy phase, there exists phases of higher density. A similar behavior is seen in 
porous silicon where transitions are observed from a uniform network of thm pores to 
shorter, wider pores on increasing the doping level in the substrate [12]. On the other 
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Figure 3.4: Plot of the percentage aggregate density vs. I (in lattice units) obtained 
in deposition on to a linear substrate of length 300 units. The two curves are for 
AW — 6 (solid line) and AW = 9 (broken line). In each case, a minimum is 
observed at AW = l. 


hand, increasing the applied potential leads to the formation of thick and relatively 
linear pipe-like pores [11]. These similarities suggest that the drift-diffusion model 
can be extended to model porous silicon formation. 

The computer simulations with an underlying diffusion limited process are thus 
seen to premise some kind of a depletion zone similar to the concept of the semicon- 
ductor depletion layer which is invoked extensively in the phenomenological theories 
of Sec. 3.2. However, to produce the complex pore geometries of porous silicon, it is 
quite evident that other competing mechanisms must also be at work. For instance, 
the formation of an interconnected network of pores as seen in p-type silicon is for- 
bidden by this very “depletion layer 1 ’ which tries to keep them separate. Hence we 
introduce a pore propagation probability which allows an etched site to propagate at 
random with a finite probability. Such a diffusion induced nucleation model (DIX i 
is outlined in the next section. 
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Figure 3.5: The diffusion induced nucleation model. Hole (filled circle) from the p- 
tvpe substrate diffuses towards the HF electrolyte and etches a surface site. The 
broken line represents the depletion layer edge at a distance AH' from the interface. 
The hatched area represents the pores. I (dotted line) is the drift-diffusion length. 
The substrate is kept at a positive potential with respect to the HF electrolyte. 

3.4 The Diffusion Induced Nucleation Model 

The experimental conditions for porous silicon formation outlined earlier (Secs. 1.2 
and 3.1) suggest that holes constitute the diffusion limited electro-active species [11]. 
A hole generated in the bulk undergoes diffusive motion until it comes within an 
interacting distance / of the surface. It is then field driven to the surface site and 
dissolution occurs. We designate l as the drift-diffusion length. 

Dissolution of an atom can create upto three dangling bonds in the nearest neigh- 
bor atoms. These dangling bonds can act as secondary electro-active species, promot- 
ing further dissolution with a finite probability p p . Thus the dissolution can either 
continue or terminate with a termination probability pt- 


Pt = l~ 4 p P 


(3.3) 
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The simulation is carried out on a two dimensional M x A* lattice with periodic 
boundary conditions (i.e., a wraparound configuration of modulo A in the lateral 
direction). Fig. 3.5 elucidates the geometry and the primary processes involved. 
Each unit of the lattice corresponds to one silicon atom with four nearest neighbors. 
The HF-Silicon interface is represented by a single side of the lattice which we term 
the ‘surface’. Close to the surface there is a depletion region where no free carriers 
are present. The algorithm used to simulate etching in p-type silicon is outlined as 
follows. 

Diffusion 

(i) Holes are generated in the lattice sites beyond the depletion region of width AH', 
with uniform probability. 

(ii) If a surface site is present within the drift-diffusion distance l of the hole, the 
hole moves to that site, causing dissolution. 

(iii) If no surface sites are located within a radius l , the hole moves to one of the 
nearest neighbor sites at random. This random walk is repeated till it moves to a 
site within a radius l of the surface. It then proceeds to dissolve the detected surface 
site. 

(iv) If the hole wanders beyond a distance 3 AW from the interface, that hole is 
aborted and a new one generated. 

Pore Propagation 

(v) Following the dissolution of a site, one of its nearest neighbors is chosen with a 
probability p p and removed. This process is repeated till termination occurs with 
probability p t where p v and p t are related by Eq. (3.3). 

Quantum Effects 

(vi) Due quantum confinement effects, holes have a smaller probability of being gen- 
erated in the porous silicon structure. We take this into account by employing a 
Boltzmann probability within a Metropolis like algorithm. 

For n-type silicon the algorithm stands modified as follows: 
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As mentioned earlier, it is observed that pores are formed in n-type silicon only 
under illumination. Therefore, we assume that optical generation of holes is the major 
controlling mechanism in this case. 

According to the Fermi thin slab approximation, intensity attenuation in silicon 
is given by 

I(z) = / 0 exp(- J u 2 ) (3.4) 

where / 0 is the incident intensity and /r is the attenuation factor associated with sili- 
con. The number of carriers generated within a region is assumed to be proportional 
to the optical intensity in that region. 

An optically generated carrier can recombine in several ways. The diffusion of 
these carriers is therefore limited by the carrier lifetime r. To take this into account, 
the hole is allowed to walk only for a proportionately smaller number of steps. 

As the surface region gets etched, confinement effects increase the band gap caus- 
ing these regions to become transparent to higher frequencies of illumination [111]- As 
a result, the carrier generation zone continually moves inwards. Taking into account 
these factors, the ‘diffusion’ part of the algorithm is modified as follows. 

(i) Holes are generated beyond the depletion region in the lattice, with an expo- 
nentially decaying probability from the surface in accordance with Eq. (3.4). 

lii) If no surface sites are found by the hole after a short number of steps (« 5), 
it is aborted. 

The rest of the algorithm proceeds exactly as described for p-type silicon. 

We also define a time scale for the simulation as follows: (i) each hole movement 
corresponds to one unit of time and (ii) each pore propagation chain also constitutes 
one time unit. This scale is used in later sections to make meaningful comparisons 
with experiment. 


3.5 Morphologies in Diffusion Induced Nucleation 

The diffusion induced nucleation model described in the previous section has been 
employed successfully to simulate structures similar to those resulting from the an- 
odization of p-type and n-type silicon. A large number of runs were made 1 with 
'The runs were carried out on a HP 900 workstation. The random number generator employed 

had a periodicity of 2 s2 . 
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varying parameters. The resulting pore geometries exhibited a strong dependence on 
the parameters AIT and 1. Typical results obtained are reported in this section. 



Figure 3.6: Typical morphologies obtained in the diffusion induced nucleation sim- 
ulations on 300 x 300 lattices, (a) Structure similar to the porous film formed on 
& P-type substrate, obtained using ATT = 6, l = 11. (b) shorter and wider pores, 
similar to those formed on p + substrate obtained using AIT = 3 and l = 4. (c) Simu- 
lation done for an n-type substrate. (AIT = 8,1 = 9). Note the relatively linear pore 
formation, compared to (a), (d) n + substrate, simulated using AIT = 3 and l = 4. 
The pore propagation probability is 0.15 in all cases. Note that (b) and (d) have 
almost identical structures. These morphologies are similar to TEM visualizations of 
porous silicon [11,12]. 

bsing the algorithm for p-type silicon, the simulations were carried out on a 
300 x 300 square lattice. With ATT = 6,1 = 11 and p p = 0.15, the structure shown in 
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Fig. 3.6a was obtained. This has a spongy nature with long and highly interconnected 
pores and resembles the TEM visualizations of porous silicon fabricated from p-type 
■ substrates [11.12]. 

To model the increase in doping concentration, one has to account for both the 
decrease in the depletion layer width and the shorter interaction distance due to 
screening. A screened Coulomb potential can be written as 

<^(r) = - exp (- k s r ) 

T 

where q is the carrier charge, r the distance from the carrier, and the inverse screening 
length k s oc 7i 0 , the number of carriers present. The depletion layer width AW 
decreases with increasing n 0 both in the Beale type depletion and in the quantum 
confinement case. For a semiconductor depletion layer, this relation is given by 
Eq. (3.1). In the quantum confinement theory, the probability of finding a carrier 
in the enhanced band gap zone increases with carrier concentration. This can also 
be modeled as a decrease in AW. To incorporate these dependences, smaller values 
were used for both AW(= 3) and l(= 4) with all the other parameters being the 
same as that in Fig. 3.6a. The resulting structure, depicted in Fig. 3.6b has wider 
and shorter pores with numerous side branches and compares well with the results 
obtained for p + material. 

In order to carry out the same procedure on n-type silicon, the algorithm is 


modified as explained in the previous section. Fig. 3.6c represents a typical structure 
that was obtained using the algorithm for n-type silicon, with ATI' = 8, / = 9 and 
p t> = 0.15. Holes were generated in the lattice according to the distribution given 
by Eq. (3.4). with the attenuation constant p = 1. The pores are no longer random 
as seen in p-type. but exhibit a strong tendency to form straight channels [11]. At 
higher potentials, this leads to “piping' 1 where we get straight, wide pores [11,99]. 
Our model can replicate this phenomenon if the depletion layer width is reduced 
(corresponding to a large applied potential) and the pore propagation probability is 
increased (a higher current would increase the etching probability). In all likelihood, 
the observed pipe formation has other processes contributing to an increased pore 
diameter, but the it is heartening to note that the basic structure can be obtained 


through our simulation. ' ‘ 

The pores in Fig. 3.6c are found to propagate in a linear fashion because of the I 
following reason. Since the carrier lifetimes are very short, the6SM:Il^ it tA^R^RY 

Z a rruut 
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through to an interior site are very small. If the holes are unable to locate a tip site 
as soon as they are formed, they recombine and hence do not contribute to silicon 
dissolution. This virtually assures pore propagation in the tip direction, leading 
to straight channels. However, in the case of n + material, the inter-pore spacing is 
smaller and the pores get interconnected. This is shown in Fig. 3.6d, which shows the 
result of a simulation carried out with AIT = 3 and l = 4, all the other parameters 
remaining the same as in Fig. 3.6c. Note that Fig. 3.6b and Fig. 3.6d are almost 
identical, in conformity with TEM visualizations [11]. 

To sum up. our diffusion induced nucleation model successfully reproduces the 
morphologies of porous silicon films grown on both p-type and n-tvpe silicon sub- 
strates. 


3.6 Discussion 


In Sec. 3.4. we have outlined the diffusion induced nucleation model used to simulate 
porous silicon formation on both p-type and n-type structures. Sec. 3.5 illustrates the 
differing morphologies obtained by varying the parameters and compares them with 
experimentally obtained structures. The values chosen for the parameters in each 
simulation were based on the experimental conditions. In this section, we attempt 
to place our efforts on a firmer footing by highlighting the agreement between the 
simulation and experimental growth literature. The results reported in this section 
are averaged over typically 10 simulations. Runs were made for lattice sizes ranging 
from 300 to 1000. Finite size effects were seen to be negligible. 

We first investigate the rate of growth of the porous silicon film as a function of 
the controlling parameters. To monitor the growth of the film, we pause to define a 
root mean square height 


1 



where hi is the pore height on the i th row of the lattice. 

Fig. 3.7 depicts the variation of the height h of the film as a function of the time 
t. Bej ond a certain initial period, the plot is linear, which is indicative of a constant 
growth rate. This behavior is in conformity with several experimental observations 
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Figure 3.7: Plot of the film thickness h vs. time t. The plot is linear after beyond a 
small initial period. This figure is for a 300 x 300 lattice with p p — 0.15. Similar be- 
havior was obtained for sizes upto 1000. Such linear rate growths have been reported 

by several workers [12.99,100]. 

which report film growth rates which are linear in time. Beale and coworkers [12] 
report that the film thickness is proportional to the anodization time. Searson [100] 
has found that this linearity is sustained over a large range of values for the film 
thickness (up to ~ aOCfian). These observations are confirmed by Lehmann and Foil 
[99], who report that the rate of growth is a function of parameters other than time. 

Further, we note from Fig. 3.7 that the rate of growth ( h ) depends on the width 
of the depletion region (AW) and the interaction distance (/) associated with the 
hole 

h = h(AW, l ) 

The growth rate is seen to increase with increasing values of l , as illustrated in Fig. 3.7 
for / = 4 and / = 6, all the other parameters being held constant. Also, h increases 
with decreasing AW, as portrayed by the two plots for AW = 5 and AW = 6. j 

The dependence of h on AW is investigated in Fig. 3.8 which depicts a log-log j 
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In AW 


Figure 3.8: Log-log plot of the rate of growth h (^) of the film vs. the depletion 
layer width. AW (l = 5). For large A IV, the graph follows the solid line with a slope 
~ —2, but exhibits a steeper slope for smaller values of AW (broken line). p p was 
kept constant at 0.15. A semi-quantitative discussion for this behavior is given in 
Sec. 3.6. 


plot of h vs. AIV. In the large AW range, h is seen to vary as 

h a AIV 2 


(3.5) 


whereas for smaller values of AW, the plot becomes steeper. This behavior maybe 
explained by invoking the dependence of the depletion region width on the applied 
potential across a semiconductor electrolyte junction. 

AW oc (V BI - V a )2 (3.61 


where AW is the depletion layer width, V BI is the contact potential across the junc- 
tion and V a is the applied potential. From Eq. (3.6) we obtain, for V 0 < V B1 


and for higher values of V a 


Va_ 

Vbi 


~ -2 In AW 



V a 




Vbi 


(3.7) 


—(2 + e) In AW e>0 


(3.8) 
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Figure 3.9: Porosity vs. depth plot for a typical simulation. The surface layer has 
high porosity (70-90%) and exhibits a fractal character [109]. This is followed by a 
constant porosity region (~ 60%). In the growth region AZ G) there is a sharp fall in 
the porosity. This simulation was done for a p-type substrate with l = 5, AW = 3 
and a 300 x 300 lattice. The porosities are comparable to those of samples exhibiting 
visible photoluminescence. The broken line shows the porosity profile when a reduced 
barrier height is employed. 

Searson [100] has shown h to be linearly dependent on the current density. In 
the pore formation region, a typical silicon I —V a characteristic [11] is exponential. 
Therefore, h should vary as A W~ K , where K is a constant. Thus we see that Fig. 3.8 
is in good agreement with Eqs. (3.7) and (3.8). 

In the context of the quantum confinement model, we expect AW to decrease 
as the voltage is increased. AW is inversely proportional to exp [-(£,, - A V)/kT). 
The barrier height is effectively lowered because of the higher carrier energy. This 
will also lead to a similar monotonic dependence of h on AW. 

This model yields much higher values for the porosity (~ 80%) unlike the finite 
length diffusion model proposed by Smith and Collins [90], which typically reported 
porosities in the range 20-40%. In Fig. 3.9 we show a typical porosity vs. depth 
plot. The surface region has high porosity (70-90%). Porous silicon samples with 
porosities in this range are capable of causing photoluminescence in the visible range 
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[1,112]. Beyond the surface region, the porosity stabilizes around 60% and remains 
depth independent. It has been reported recently that further passive etching (i.e., 
HF can dissolve Si very slowly in the absence of any anodizing current) can result in 
depth inhomogeneities [113], but that regime is not of interest to the present work. 
The growth zone, denoted by A Z G in Fig. 3.9 shows a sharp fall in porosity. The 
dependence of porosity on quantum confinement is illustrated by the dotted line. We 
assume a smaller barrier height in our simulation, which results in a higher probability 
of hole occurrence in the porous silicon layer. This leads to higher porosity values as 
argued by Lehmann [108]. 

The high porosity surface region exhibits a fractal character, as has been reported 
earlier by other workers [109]. This region is probably mainly responsible for the 
high frequency component of the photoluminescence spectrum of porous silicon. Ex- 
periments have reported the suppression of the high energy part of the spectrum if 
the surface layers are made nonradiative [114]. It has also been shown that such 
inhomogeneities can lead to a ‘Lifshitz tailing’ of the photoluminescence spectrum 
[59] (See discussion in chapter 4). 

It would be desirable to define a statistical measure which can differentiate be- 
tween the morphologies depicted in Fig. 3.6 in a quantitative fashion. As seen in 
Fig. 3.9, the porosity is nearly depth independent beyond the initial high porosity 
region. Hence a mass-depth scaling exponent will have value ~ l. This points to 
a possible non-fractal nature of the pore network. Hence the use of a conventional 
fractal dimension may not prove fruitful. We examine the degree of branching in 
terms of a dimensionless quantity B which is defined to be the ratio of the distance 
between two adjacent nodes (branching points) on a pore to the mean thickness of the 
pore. Preliminary calculations yield approximate values of B as follows: for Fig. 3.6a 
(p-type silicon) B 9, for Fig. 3.6b (p + ) B ~ 3, for Fig. 3.6c (n-type) B « 17 and 
for Fig. 3.6d (n + ) B ss 5. 

Unlike other models, our simulation can successfully reproduce the electropol- 
ishing regime in the anodization of porous silicon. At high anodic overpotentials, 
it has been observed that silicon uniformly etches away instead of forming porous 
structures. In our simulation, this corresponds to small values of AW, which leads 
to structures with porosities greater than 90%. The simulation carried out using 
AW = 2, l = 5 and p p = 0.15 exhibited a uniform removal of the surface. The etched 
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region in this case resembled the aggregates obtained using the ballistic deposition 
model [71], 

Eailier models attributed electropolishing to the formation of an insulating oxide 
at the pore tips, which leads to preferential dissolution in the lateral directions. In 
the finite length diffusion model [90,115] the pore density is a function of the diffusion 
length AW. Smaller values of AW result in high porosity values. The limit AW -> oo 
mimics DLA algorithms-, whereas the AW -> 0 limit yields the Eden [80] cluster. The 
compact Eden pattern may be invoked to explain electropolishing, though this is not 
the contention of Smith and coworkers. However, no convincing physical correlation 
has been established between varying the value of AW and changing the experimental 
conditions. While this may be one of the processes that is taking place, we emphasize 
that in our model, electropolishing occurs as a natural consequence of decreasing AW, 
without having to invoke any additional parameters. The dependence of porosity on 
l is illustrated in Fig. 3.4. 

Besides the diffusion limited approach, percolation models have also been em- 
ployed for simulating pore growth in silicon [110,117]. Given the nature of the Si-HF 
system where the acid eats into the crystal lattice, a percolation like approach seems 
appealing. However, such simulations have not been as successful as their diffusion 
limited counterparts in reproducing realistic morphologies. Further, no serious at- 
tempt was made in any of these simulations to establish a meaningful correspondence 
with experimentally relevant parameters. The present work (a) successfully repro- 
duces porous silicon morphologies similar to those observed in TEM, including the 
well known ‘piping’ phenomenon in n-type silicon, (b) obtains structures of porosity 
60-80% which in porous silicon are responsible for visible photoluminescence, (c) ac- 
counts for constant growth rates and its dependence on the anodizing potential, and 
(d) explains electropolishing in a natural fashion without taking recourse to additional 
processes or parameters. Further, the role of quantum confinement in determining 
the porosity is illustrated. 

We have presented a simple, pragmatic diffusion induced nucleation model which 
explores the generic features in porous silicon formation. Our model successfully 
captures the essential physics of the dissolution process by employing mainly three 
parameters, namely the depletion layer width AW, the drift-diffusion length l and the 
pore propagation probability p p . It must be noted that the simulation is of a classical 
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nature. The quantum confinement effects are included by invoking a depletion region 
wherein hole generation occurs in a probabilistic fashion. These effects are seen to 
play a role in determining the lower limit of the crystallite size (Fig. 3.9). Further 
investigations need to be carried out on the model to work out better correspondences 
with the experimental conditions. For example, it may be worthwhile to look into the 
pore morphologies when inert masks are placed on parts of the silicon surface [97]. 
This will enable us to determine the feasibility of fabricating optical interconnects 
using porous silicon, a procedure which requires sharp edges. At a theoretical level, 
the fractal character of the surface layer may also be explored. Det ailed experiments 
could be conducted to verify the dependence on the depletion region width, doping 
and temperature. 
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Chapter 4 


Shape of the Luminescence 
Spectrum in Porous Silicon 


4.1 Introduction 

Primarily, porous silicon is a disordered system. It consists of an intricate network of 
crystallites with varying sizes and shapes. The large surface supports hydrogen and 
oxygen complexes. At the microscopic level there exists dangling bonds and voids. 
Constructing a theory for the photoluminescence (PL) spectrum which encompasses 
all levels of disorder is a difficult task. Further, it may not be a desirable goal 
as it would entail the use of a large set of parameters whose numerical values arc 
not accessible either by microscopic calculations or by experimental observations. 
Hence, a modest approach to explain the overall features of the PL spectrum, which 
operates within a circumscribed and limited set of plausible parameters is a desirable 

alternative [59]. 

The salient features of the PL spectra of porous silicon have been outlined in 
chapter L We briefly recapitulate them below (also see Fig. 4.1): 

. The dominant red band in the PL spectra falls in the energy range 1.2 - 2.3 
e V. This constitutes about 97% of the observed luminescence intensity. 

. The PL spectra are broad, with FWHM around 300 - 400 meV. • 

. The spectral line shape is asymmetric on the energy scale, with a shoulder on 
the high energy side. 
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Figure 4.1: Salient features of the PL spectra in porous silicon. Note the asymmetry 
in shape, with a shoulder on the high energy side. The FWHM is approximately 300 
meV. Some spectra show a low energy tail, which is represented by the broken line. 

• A small pumber of experimentally reported spectra show a low energy tail, 
which extends into to the infrared. 

Our work is based on the quantum confinement model [1,19,103]. Most workers 
employ a simple quantum confinement model (Sec. 1.3.1) where the PL peak alone 
is sought to be explained. We believe that disorder plays a key role and model it by 
a distribution of crystallite sizes. This work, described in Sec. 4.2, is similar in spirit 
to the mean field approach of Kane and the probabilistic arguments of Lifshitz to 
explain the Urbach tail in optical absorption [118,119]. A large body of theoretical 
work in the Urbach tail problem assume that the absorption edge a(E) is simply 
prbportional to the electronic density of states and ignore the energy dependence 
of the transition matrix element and the electron-electron interaction [120]. There 
are several features which we choose to ignore in our formalism. Some of these are, 
relaxation of carriers, gap states due to voids and defects, thermal disorder and the 
distinction between hole and electron contributions to oxeiton energies. Wo wish to 
employ a small set of plausible parameters, and further, observe that these explain 
several reported PL spectra (Sec. 4.3). Our formalism can be extended to include the 
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above mentioned factors as well as important experimental parameters such as the 
frequency and intensity of the incident excitation. These are indicated later (Sec. 4.4). 

In Sec. 4.2, we describe our theoretical formalism. We derive expressions for PL 
spectra originating from both column and dot distributions. The presence of local 
inhomogeneities where either columns or dots might dominate requires a Lifshitz like 
argument for modeling the edges of the PL spectrum. The PL spectra, which may 
appear Gaussian when plotted against the wavelength has a distinct asymmetry on 
the energy scale. This is naturally obtained in our formalism. A further noteworthy 
feature is a shift in the PL peak due to the distribution of crystallite sizes (Fig. 4.2). 
This enables us to employ a small and physically reasonable exciton binding energy. 

In Sec. 4.3 we outline how the numerical values of the parameters employed are 
obtained from experimental observations and microscopic calculations. Our calcula- 
tions compare favorably with the PL spectra obtained by several workers. 

In Sec. 4.4, we briefly discuss earlier works based on the quantum confinement 
model. We outline extensions of our formalism to include some of the features which 
have been neglected in Sec. 4.2. We also indicate how some insight into related ex- 
perimental observations such as photoluminescence excitation (PLE) can be attained 

based on our model. 


4.2 Theoretical Framework 

Our aim is to explain the photoluminescence (PL) spectra employing a minimal set of 
broad and plausible assumptions. There exists in literature a variety of approaches to 
the Urbach tail in optical absorption. We propose to work along the lines of two sister 
approaches : the one due to Kane and the well known Lifshitz argument [118,119]. 
We base our arguments on the quantum confinement model (see Sec. 1.3.1) but avoid 
the oversimplifying assumptions made by some of its proponents. 

Columns of mean diameter d 0 in the nanometer range have been reported by 
several independent groups [103,121,122]. The growth of these columns is a stochastic 
process as discussed in chapter 3 and it appears reasonable to assume columns of 
silicon with a Gaussian distribution of diameter d centered around a moan d 0 , with 
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variance a. 

{d - d 0 ) 2 ' 

2a 2 

The number of electrons in a column of diameter d participating in the PL process 
is proportional to d 2 . The heights of the columns depend only on the growth time 
and are approximately the same. Hence 

N e = N e {d)=ad 2 (4.1) 



where a is a constant. 

For a PS sample consisting of varying column diameters the probability distribu- 
tion of electrons participating in the PL process is given by a product of the above 
two expressions: 


Ped = 



(d - d 0 ) 2 ' 

2a 2 


(4.2) 


where b is a suitable normalization constant. 

In the quantum confinement model, the PL process is attributed to the energy 
upshift of the electrons and is proportional to 1/d 2 . 

The PL energy hu> is given by 


ft,U) — Eg — Ef, + — 
or 


where E g is the bulk silicon gap (1.17eV), Eb the exciton binding energy and c, an 
appropriately dimensioned constant. The energy upshift due to confinement AE U is 


A B. = huj - (E, - E t ) = £ 


(4.3) 

(4.4) 


where we have also paused to define a mean upshift A E 0 related to the mean column 
diameter do- 

The PL line shape is then determined by transforming Eq. (4.2) to the energy 
axis as is commonly done. 


P(AE U ) = 



(d-d q ) 2 

2a 2 


c\(d) 



4.2 Theoretical Framework 


63 



Figure 4.2: Theoretical PL spectra. The solid line is due to columns (peak at AE P = 
0.46 eV) and the dashed line is due to dots (peak at A E p = 0.67 eV). The mean 
diameter is d 0 = 30 A and variance a = 4 A for both spectra. The quantum 
confinement constant c (in c/d 2 ) = 485.816 eV.A 2 (the same as calculated by Read 
ct al. [39]) and the binding energy Eb ~ 0.14 eV. The horizontal bar at the peak of 
the column spectrum indicates a downshift (~ 0.1 eV) from hoj 0 = c/dl on account 
of the statistical distribution of the column sizes. (See text in Sec. 4.2 for further 
discussion). Both spectra are normalized to unity. 


The Dirac delta function facilitates a straight forward integration using the rela- 
tion 



5{f{x))g(x)dx = Y,9(xi) 

Xi 


>/ 


dx 

X 


) 


where .r; are the roots of f(x). This yields, 


P(AE U ) *. 


1 b ( c \ 3/2 c:: 
yJ/kKcr 2AE U \AE U ) 



(4.5) 


The PL line shape is approximately Gaussian if a is small. For finite a the 
\J\JaEu factor in the exponential outweighs the polynomial dependence in the pre- 
factor, resulting in an asymmetric curve with the shoulder on the high energy side 

(see Fig. 4.2). 




64 


Shape of the Luminescence Spectrum in Porous Silicon 


Another aspect we need to note is that the mean energy of the upshift AEq 
(Eq. 4.4) and the location of the PL peak A E p {= hu> p -{E g - E b )) are not identical. 
To see this, we differentiate and set to zero the above expression (Eq. 4.5). This 
yields 

2 


A E p = AEq 


1 

10 


do 
, a 



1/21 


■ 4 * 20 


T 2 


(4.6) 


For a /d 0 -> 0, A E p = AEq as expected. However, for reasonable a the above 
expression can be Taylor expanded to yield 


A£ * = AB »(i" 10 © 2 ) 

From the values reported by Read et al. [39]: d 0 = 30 A, a a 3 A, a/do ~ 0.1, 
A E p = 0.9AE 0 . 

Thus, there is a downshift. This is depicted by a horizontal bar in Fig. 4.2. It 
is important to realize the physical significance of this. As will be explained in the 
next section, the downshift dispenses with the need to invoke large and physically 
unreasonable exciton binding energies. 

The peak in PL intensity is P(AE p ). We can obtain an approximate expression 
for the FWHM (A E FW hm) of the PL spectrum. We can write the expression for the 
PL spectrum (Eq. (4.5)) as 


P{AE U ) = 


K 


AEu 


, 5/2 


exp 


AEq (do 

2 U 


1 


.y/AK VOTE 


where K is an appropriate normalization constant. If the pre-factor energy depen- 
dence AEj is ignored, this expression represents a Gaussian distribution of the 
variable £ = 1/ y/AE u , with a variance a = a/(d 0 y/AE 0 ). The FWHM on the £ scale 
is then given by 2a. Transforming this to the energy axis, we find, 

ab fwkm = ae 0 [{i-(£)} '-{i +(!;)} 

For small (a /do), 

AEpwHM — 




(4.7) 
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A similar result can be obtained using a simple quantum confinement approach. 
Noting that A E u = c/d 2 (Eq. (4.3)), 


AE fwjim = 2 | 6{AE U ) \\ d=:do 
4 c6d 
~ . dl 

4AEq<j 


The above expressions are identical for large ( d 0 /a ). Using c = 485.816 eV A 2 , 
dd = ct — 3 A, and rfo = 30 A(see the next section, Sec. 4.3 for a discussion of the 
numerical values), we obtain from Eq. (4.8), AE FWHM ~ 200 meV. A larger and 
experimentally reported AE FWHM is obtained if the exact expression Eq. (4.5) is 

employed. 

Some workers have hypothesized the existence of dots instead of columns of silicon 
[58,123]. An analysis similar to the above can be carried out for dots. Employing 
N r a d 3 in Eq. (4.1) and carrying out the integration, the line shape for the dot is 


^dots^-^") ~ *^£3 6X P 





(4.9) 


where K is a normalization constant. 


4,2.1 The Low Energy Tail in PL 

It is conceivable that both columns and dots are present in porous silicon. The exper- 
imental PL spectrum is a weighted sum of the column and dot spectra. Statistically 
speaking, there would be inhomogeneities resulting in a region (say of volume V) 
having an overwhelming column concentration c c (> c c0 , the mean column concen- 
tration), The dot concentration in this region Cd (•< Cao, the mean dot concentration). 
One can then invoke a methodology along the lines of the Lifshitz argument [119]. 

One attempts to evaluate the probability of occurrence of n column crystallites 
in a volume V, which contains a total of N crystallites. This is given by 

^ - 55 ( 

= cxp[ln A! - In n! - In (N -»)!■+ h In c c0 + (N - n) In c rf0 ] 
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Making use of the Stirling approximation IniV! ~ iVlniV - N, we can write, 


P(n) = exp 


N (c c In — + Cd In — \ 

\ c c o Qo / J 


Note that we have substituted n = Nc c and (N — n) — Nc-d ■ Taking the limit Cd — * 0, 
we get 


P(n) ~ exp 


-Nc c In 


^c0 


If v c is the mean volume of the column crystallite in this region, then 

V 

N a - 


(4.10) 


(4.11) 


cx d 2 for a column crystallite, hence the confinement energy varies as (Eq. 4.3), 


A E u a l/v c 


(4.12) 


Using Eq. (4.11) and (4.12) in Eq. (4.10) 

P(n) ~ e.xp [— A'i AE U ] (4.13) 

where K\ is a constant. The low energy part of the PL spectrum is modulated 
by this factor, and hence exhibits a tail. The presence of this tail accounts for 
some discrepancy between the theoretical PL spectrum based on Eq. (4.5) and the 
experimentally reported ones. 


4.3 Results 

We shall now present calculations based on our expressions in the previous section 
and compare them with experimentally reported photoluminescence spectra. Almost 
all reports on porous silicon report at least, one PL spectrum. We have selected a 
representative set which highlights both the promises and problems associated with 
the natural theoretical framework outlined in the previous section. Variations in the 
proposed model and alternative explanations are outlined in the next section. 

We note that the majority of the PL spectra reported have been recorded against 
the wavelength A on the x— axis. In our comparison, we have faithfully transformed 
A to the energy hu. Note that hu (in eV) = 1.24/A where A is in fxm. Recall that 
hu = AE U - ( Eg — Eb) from Eq. (4.3). 
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Figure 4.3: Comparison of experimental (broken line) and theoretical (solid line) PL 
spectra. The experimental curve, from Cullis and Canhain [124] (their Fig. 1, sample 
3), has also been discussed by Read and coworkers [39]. The theoretical spectrum is 
obtained using a mean diameter d 0 = 30 A and variance a = 4 A for both dots and 
columns. The dot and column concentrations are 10% and 90% respectively. 


The numerical values of physically important variables in our calculations are as 
follows: The band gap of silicon E g ranges from 1.17 eV at 0 K to around 1.14 K at 
300 K The existence of a distribution of column diameters results in a downshift from 
the mean peak(see the horizontal bar in Fig. 4.2). This allows us to use a smaller 
exciton binding energy ( E b * 0.06 eV) which is physically more reasonable (than, 
for example, 0.16 eV proposed by Read and coworkers [39]). Thus E - E b - 1.10 
eV. The constant c associated with the confinement energy (- c/d is taken rom 
the calculations of Read et ai [39]. Its value is 485.816 eV-A 2 . 


Electron microscopy measurements suggest that the crystallite sizes are in the 
range 20 - 40 A [39,124]. The distribution can then be modeled as a Gaussian, with 
variance a » 10/3 a 3 - 4 A. Thus 99% of all columns or dots lie within 3<r. We take 
* ~ 30 A, a number reported by severai workers [39,103,124], The speerhe values of 
R,cr} for a given spectrum are cited in the figure captions and these happen to he 



68 


Shape of the Luminescence Spectrum in Porous Silicon 


close to the above quoted values. In all our calculations we take the mean diameter 
(d 0 ) of the column and/or dot to be the same. A priori, there is no reason for them 
to be different. 

In Fig. 4.3 we compare the theoretical spectrum (solid line) with one experimen- 
tally obtained by Cullis and Canham (broken line) [124]. The experimental spectrum 
with a peak at hu p = 1.48 eV, a full width at half maximum (FWHM) of 325 meV 
has also been theoretically studied by Read and coworkers [39]. The theoretical spec- 
trum is obtained with mean diameter do = 30 Aand a = 4 A. These values are the 
same as the ones suggested by Read et al. and so no “fitting” or “adjustment” on 
our part has been carried out. The agreement is excellent ( Tiuj v = 1.45eV, FWHM = 
300 meV) except at the low energy end. The Lifshitz-like argument outlined at the 
end of the previous section provides a substantial improvement at the lower end with 
constant K\ = 0.13 eV -1 . We do not wish to overemphasize its importance at this 
juncture except to indicate that it suggests inhomogeneities in porous silicon. 

In Fig. 4.4 the experimental PL spectra of Vial and coworkers (broken line) is 
compared with our theoretical calculations (solid line) [58]. The experimental PL 
spectrum has a peak at hu p = 1.54 eV with a FWHM = 270 meV. It was obtained 
for an oxidation rate of Q 0 /4 where Q 0 is a threshold exchanged charge and is pro- 
portional to the layer thickness. The theoretical PL spectrum was obtained for a 
mean diameter of do = 28 A whereas the experimental work quotes a theoretically 
calculated value of 30 A. The variance <7 = 4 A(same as for Fig. 4.3). Vial et al. pro- 
pose a purely quantum dot model. We need to invoke dot and column concentrations 
of 15% and 85% respectively with both dots and columns having the same {do, < 7 } as 
mentioned above. Our calculations are in excellent agreement with the experimental 
curve. 

In Fig. 4.5 an experimental PL spectrum (broken line), due to Zhang et al. is 
compared with our theoretical calculations (solid line) [125], The theoretical spectrum 
is obtained on employing mean diameter d 0 = 28 A, variance a- 3 A, assuming that 
only dots are present in the sample. The peak position is consequently much higher 
(- !- 9 eV ) compared to the previous figures. The agreement with the experimental 
curve is good. 

The above results have been demonstrated for representative experimental PL 
spectra. We have also performed a number of calculations for experimental spectra 
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Figure 4.4: Comparison of experimental (broken line) and theoretical (solid line) 

PL spectra. The experimental curve, is due to Vial and coworkers (their Fig. 2, 
with oxidation level Qq/A) [58]. The theoretical curve was generated employing the 
presence of both dots (as conjectured by Vial and coworkers) and columns, each with 
mean diameter do = 28 A and variance cr = 4 A. Dots and columns are present, the 
concentrations being 15% and 85% respectively. 

reported by several other workers in the field. These include the ones by Banerjee et 
al, (their Fig. 2, curve c) [126], Hummel and Chang after spark erosion (their Fig.3) , 
[127], Tsai et al on Si:H passivation and immersion in HF (their Fig.l) [128], and 
Jung et al. (their Fig. 3) [ 129 ]. The numerical values of the parameters used have ; 
not been adjusted and are the same as those outlined in the beginning of this section. , 
The agreement in all cases is good. 


4.4 Discussion 

The simplified assumption of free standing quantum wires due to Canham [1] cannot 
explain the broad photoluminescence (PL) spectra in porous silicon, which typically 
exhibits a FWHM of 300 - 400 meV. Though there has been suggestions in the 
literature that a distribution of crystallite sizes may be responsible for this [45], 
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Figure 4.5: Comparison of experimental (broken line) and theoretical (solid line) PL 
spectra. The experimental curve is from Zhang et al. [125] (their Fig. 1, sample g). 
The theoretical curve (Eq. (4.5)) assumes a distribution of quantum dots alone with 
mean diameter d 0 = 28 A and variance a = 3 A. 

formal calculations have been attempted by only a few groups. In this section, we 
briefly review these attempts and examine our work critically in comparison with 
them. 

In an earlier model, Fishman et al. [130] model the PL spectrum on lines simi- 
lar to the work described in Sec. 4.2, but assume equal population of states in the 
crystallites. The energy upshift A E u is taken to be 


where n= 2 if a simple infinite potential well is assumed, and K is a constant calcu- 
lated from an effective mass approximation. They further account for a variation in 
the nonradiative rate with the photon energy E , 

Rn(E) oc exp ( E/p NR ), 


where p^ii is a constant (=0.29 eV). The PL spectrum is then given by 


m = 


P(E) 

ME)' 


( 4 . 14 ) 
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where 


P(E)dE = P(d)d(d), 


P(d) being the crystallite size distribution, which is assumed to be Gaussian. They 
plot (4.14) for two different n (n = 1.64 and n = 2). The FWHM in both cases 
is seen to compare well with the experimentally obtained values. However, the low 
energy tail perceived in the PL spectra is not reproduced accurately by (4.14). 

Xie and coworkers [31] account for the nonradiative recombination by modeling 
the probability that a given crystallite can luminesce according to Poisson statistics 


P r {d) = exp 



(4.16) 


where (d) 3 is the mean volume occupied per dangling bond. Using this modulating 
factor along with a size distribution, the emission spectrum is calculated and is found 
to yield a PL efficiency of 0.5%, which is slightly lower than the generally observed 
external absolute quantum efficiency. The expression (4.15) may not be appropriate, 
since the number of non-radiative centers are more likely to depend on the surface 
area than the volume. In this case the exponent of d could lie between 2 and 3. In 
addition, if the oscillator strength dependence is included in this calculation, the PL 
efficiency can fall further. 

Our work, described in Sec. 4.2, has attempted to obtain analytical expressions 
for the lineshape by assuming a Gaussian distribution of sizes and a simple quantum 
confinement model. The calculations have been performed for both columnar and 
spherical crystallites. One could also consider the case where A E oc d~ n , 1 < n < 2, 
as suggested by several workers [47,130] or consider polynomial dependences 


A E u 


Cl C2 

d I d? 


as reported by others [46]. In the extreme case of a d 1 dependence, Eq. (4.5) is 

modified to 


P(AE U ) 


N 

A El 


exp 



/ AE 0 

LU E u 


- 1 


which has a comparatively narrower FWHM. ( N is a normalization constant). 
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The oscillator strength f vc of the radiative transition may also vary with the 
crystallite size. Effective mass theory based approaches claim a size dependence of 
the form f vc oc d~ 6 for the oscillator strength [131,132], whereas some other workers 
claim a dependence of the form f vc oc d~ 5 [45]. LDA based calculations [42, 4G] ascribe 
a strong, but non-monotonic size dependence to the oscillator strength. We ignore 
the variation of oscillator strength with crystallite size in our formalism so that a 
minimal set of parameters can be employed. Despite this simplified approach, the 
success of Eq. (4.5) and (4.9) in modeling the observed PL spectra is heartening. 

There are other attempts to model the effect of the size distribution on the PL 
spectrum. Behrensmeier and coworkers attempt to model the PL spectra as a su- 
perposition of distinct peaks due to discrete wire diameters [133]. Wang et al . , in 
a similar calculation, obtain pinning of peaks at certain discrete energies [49]. The 
growth process is however, stochastic, and wire diameters have a range as observed in 
transmission electron microscopy and scanning electron microscopy [39,103]. Hence 
our approach has, as its basis, a continuous distribution of sizes. 

, The irregular shape of the nanocrystallites are also seen to contribute to the gap 
enhancement. Electronic structure calculations employing a randomly porous cluster 
[52] report a disorder induced enhancement of the gap in addition to the upshift due 
to confinement. This can be understood as follows. In order to minimize the ground 
state kinetic energy, it is necessary for the wave function to have near zero amplitudes 
in the constricted regions such as “necks” or protrusions in an irregular structure. In 
a recent study [134] Sapoval, Russ and Chazalviel have shown that this can lead to 
localization of the wave function in the broader regions of the crystallite. Thus, the 
presence of disorder can reduce the ‘effective size’ of the crystallite as perceived by 
the carrier, and lead to a further enhancement in the energy gap. Thus, the disorder 
in shapes may also contribute to a broad PL spectrum. 

One approach is to fit the PL spectra to a single or series of peaks with a Gaussian 
convolution. Employing this, Narasimhan and coworkers [61] have reported three 
peaks. A sum over exponentials representation is non-unique [135]. Further, the 
physical origin of these peaks needs to be elucidated for such an endeavor to be 
meaningful. In the absence of this, one merely has a curve fitting exorcise. In contrast, 
our approach assumes a physical basis wherein not the spectra, but the underlying 
structure of crystallites has a Gaussian distribution. 
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In an earlier work, Murayama and coworkers [66] interpret the broad luminescence 
band by a model in which the electron hole pair created in the silicon nanostructure 
recombines ladiatively through an excited luminescence center with strong phonon 
coupling. 1 he excited lattice center has an energy minimum away from the ground 
state minimum, because of the lattice distortion induced by strong coupling with 
phonons. I he visible transitions occur to the ground state according to the Franck- 
Condon principle. This can lead to an enlarged spectral broadening because of the 
vibration of the atoms around the luminescence center. It must be noted that lattice 
vibrations yield a typical broadening ~ kT (< 25 meV) which is too small to explain 
such a large FWHM. An explanation of this kind would involve unphysical values 
of the Huang-Rhys factor and is improbable. Further, the majority of PL spectra 
observed in porous silicon show an inhomogeneous broadening as opposed to the 
Gaussian lineshape assumed by Murayama et al. 

We observe that in seeking an agreement with the reported PL spectra, we need 
to invoke both columns and dots. Further, as mentioned in Sec. 4.2.1, a Lifshitz like 
argument could ensure that the tail states are properly accounted for. Thus, local 
inhomogeneities where columns pre-dominate over dots (or vice versa) exist. 

There has been speculations about a possible fractal character for porous silicon 
[11], Lehmann [136] points out that the experimentally observed size distributions of 
crystallites cannot account for a fractal character. It must be emphasized that the 
mere fact that the structure is porous and random does not imply a fractal geometry. 
The porosity is reported to be almost depth independent beyond the nanoporous 
layer [101]. This has also been supported by our simulation studies discussed in 
chapter 3. This results in a mass-depth scaling exponent of one, which implies a; 
non-fractal character. Any possible fractal geometry will be present only in the top 
nanoporous layer. Fishman, Mihalescu and Romcstain [130] in their work on PL 
broadening assume a fractal character and a resultant crystallite size distribution of 
the form d~ A > where d is the crystallite diameter and d f the fractal dimension. The: 
computed PL lineshapes show an unphysically large FWIIM. Because of this they; 
rule out a fractal nature. ; 

Our model can be extended further to incorporate several phenomena associated 
with the luminescence process, but at the cost of analytical transparency. Recom-j 
bination to defect states is a distinct possibility. Defect states, primarily those due 
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to dangling bonds and voids, lie in the band gap (see Sec. 1.3.3). These states may 
also be broadened due to the existing disorder in PS. Relaxation processes where the 
excited electron decays to an appropriate energy state and subsequently recombines 
is perhaps also present. The distribution of sizes in PS would exercise its influence 
over physical phenomena other than the PL spectrum. Murayama and coworkers 
[66] have reported an exponential Urbach tail in the PLE spectrum and this may be 
attributed to the above mentioned distribution and inhomogeneities. The effect of 
the surface contribution to the PL may be accounted for in our model by modifying 
the pre-factor in the distribution in Eq. (4.2) from bd 2 to (bid 2 — bid 0 ) where bi and 62 
are volume and surface dependent constants respectively. If the surface has a fractal 
character, the exponent a may be different from unity. Further, Eq. (4.3) for the PL 
energy may be modified to 

hu> = Eg — Eb + — + £(I,u>j,' • ■) 

Here, £ accounts for non-linear processes arising out of the influence of experimental 
parameters such as the intensity I and the frequency c oj of the incident radiation. 
These extensions which are not included in our work may be relevant in explaining 
the PLE spectrum and an associated exponential Urbach tail. We hope to extend 
our work to explain these and other experimental observations. 

In conclusion, we emphasize that the attractive aspect of our model is the choice 
of a minimal set of parameters whose numerical values are dictated not by the exi- 
gencies of an individual PL spectrum, but by independent microscopic calculations 
and experiments. 
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Chapter 5 


Model for the Photoluminescence 
Behavior of Porous Silicon 

5.1 Introduction 

A great deal of attention has been devoted to an understanding of the microscopic 
mechanisms responsible for the visible PL from porous silicon [2]. A considerable 
amount, of data has accumulated from investigations of the dependence of PL on the 
temperature, pressure, excitation energy, luminescence energy etc. The present work 
is an attempt to explain and systematize this large body of data. 

The properties of porous silicon have been known to depend on a range of factors. 
[2,57]. Briefly, some of these are: preparation conditions, size distribution of the 
crystallites and the surface chemistry of the nanostructures. Nevertheless, there 
is near unanimity about the photoluininescence behavior of high porosity samples 
(> 60% porosity). Further, a careful study of the literature reveals systematic trends 
in the optical properties under various conditions for the dominant red band (hu> € 
[1.2 - 2.2] eV). These trends have been discussed earlier in sections 1.3 and 1.5. We 
recapitulate some of them below: 

1. The PL intensity shows a maximum with temperature, the peak being in the 

range 50K - 150K 

2. The luminescence decay time falls by an order of magnitude or more as the 

temperature is increased from ~ 10K to room temperature. I 
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3. The PL intensity falls steeply with pressure. 

4. The luminescence decay time decreases systematically with emission energy in 
the range 1.4 - 2.5 eV 

5. The PL peak position exhibits an anomalous behavior (both blue and red shifts) 
with temperature. 

6. The PL peak position exhibits an initial blue shift and a subsequent red shift 
with pressure, when the pressure transmitting medium is alcohol. 

7. The electroluminescence from porous silicon is weak, the quantum efficiency 
being at least 3 orders of magnitude lower than that of PL. 

8. Radiative efficiency of porous silicon is higher than that of GaAs, but the lumi- 
nescence lifetime is in the micro to millisecond range unlike GaAs (nanosecond 
lifetime). 

Studies on carrier transport in porous silicon have asserted that the current pre- 
dominantly follows a surface mechanism [137-139], The carriers localized on the 
surface can (i) recombine radiatively, (ii) recombine nonradiatively, or (iii) escape to 
another localized state from where they may recombine either radiatively or nonra- 
diatively. Since the electroluminescence is found to be weak despite a strong surface 
current, one can assume that the majority of charge transfer occurs through either 
process (ii) or (iii). 

The luminescence decay time is normally expressed in terms of the competing 
radiative decay dynamics and nonradiative decay dynamics [58]. In recent years, 
some workers [140] have hypothesized the existence of a hopping term which also 
plays a role in determining the luminescence decay time. The resultant luminescence 
lifetime r can be written as 

— = Rr + R n 4- Ri lop 

T 

where Rr represents the radiative recombination rate, Rn the nonradiative recombina- 
tion rate and R hop the hopping escape rate. Earlier studies [141,142] of the radiative 
process posit a temperature dependence of the Arrhenius type 

Rr = u r exp (-jr) 
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where T r is a characteristic activation temperature and v r a characteristic frequency. 
Since the photoluminescence efficiency of porous silicon is high, we propose that the 
nonradiative recombination term R n is small compared to the hopping term R hop . 
Sevei al workers [143,144] report that the conductivity of porous silicon exhibits a 
Beithelot type [14o] temperature dependence. Therefore, we now hypothesize that 
the hopping term varies with temperature as 

Rhop = v B exp 

where 1 B is the characteristic Berthelot temperature associated with the escape pro- 
cess and v B a characteristic frequency. The varied optical behavior can then be 
analyzed in terms of the competition between the radiative and the hopping process. 
This approach is similar to the recent suggestions of a competition between radiative 
and non-radiative processes [140,146-148]. 

Section 5.2 describes the proposed model. In section 5.2.1, we describe an ex- 
pression for the time integrated PL intensity in terms of {T t ,T b } and the reduced 
frequency uq = u B /v r . Following Hurd [149], we then relate this expression to the 
properties of the barrier namely its height, width and frequency of vibration. This is 
described in section 5.2.2. In this section, we also discuss and circumscribe the range 
of numerical values that the parameters {v r ,v B ,T r ,T B } may assume for porous sil- 
icon. This a priori assignment of numerical values enables our model to possess a 
firm predictive status. 

In section 5.3, we analyze the temperature and pressure dependence of the PL 
intensity. Section 5.3.1 explores the temperature dependence of the time integrated 
luminescence intensity based on our model. It is found that the PL intensity exhibits 
a maximum at T m ~ 50 - 150K. In section 5.3.2, we investigate the temperature 
dependence of the luminescence decay time and show that it falls by an order of ; 
magnitude or more as the temperature is raised from 10K to 300K. In section 5.3.3, - 

we demonstrate that the PL intensity falls exponentially with increasing pressure. | 
In each of these investigations, we make a detailed comparison with experiments i 
conducted independently by a number of workers across the world. j 

A number of microscopic mechanisms have been proposed to explain the photolu- 
minescence in porous silicon. A popular microscopic model is the quantum confine- 
ment model which assumes an enhanced, direct band gap due to carrier confinement I 
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in the nanocrystallites. In our scheme as outlined in sections 5.2 and 5.3, we do not 
commit ourselves to any specific microscopic model. In section 5.4, we interpret our 
scheme to include an energy upshift due to quantum confinement. This enables us 
to explore the dependence of lifetime on luminescence energy (section 5.4.1) and the 
dependence of the peak position on temperature (section 5.4.2) and pressure (section 
5.4.3). In each case we again carry out wide ranging comparisons with experiment. 

Section 5.5 constitutes the discussion. We examine the validity of the approxima- 
tions made by our model and suggest further experimental verification. We demon- 
strate that our model is analytically transparent and can successfully explain a wide 
variety of experimental measurements on the dominant red band in nanocrystalline 
and porous silicon. A preliminary account of our work along these lines has been 
reported earlier [150]. 

5.2 The Model 

5.2.1 The Photoluminescence Intensity 

The intensity of the PL line is expressed as 

/(f) = N(t)R r 

where R r is the, radiative recombination rate and N(t) is the population of the excited 
carriers at time t, 

N(t) = N „ exp [-(!)'] (5.1) 

Here r is the characteristic lifetime associated with the Kohlrauseh (stretched ex- 
ponential) decay. Porous silicon is a disordered system and the time dependence of 
the PL intensity has been reported to be of the Kohlrauseh type with 0 < p < 1 
[58,140,141,151], Thus 

/(f,T) = JV„fl r exp[-( ; t) P ' 

The time-integrated luminescence intensity is then given by 

roo 

/ ( r ) = J 0 I(t,r)dt 
iVor^Rr 
P 
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Note that the above expression reduces to N 0 i^r for the pure exponential case (p = 
1). Defining I 0 = N 0 T(^)/p we have 


J(r) = I 0 R rT 


( 5 . 2 ) 


I he characteristic lifetime is determined by the radiative recombination rate (Rr), 
the nonradiative recombination rate (/l n ) and the hopping escape rate ( Rhop )• As 
mentioned in section 5.1, R n is taken to be small compared to Rhop • We can then 

express the luminescence decay time as 


The radiative rate has a weak Arrhenius dependence [142] 


R r = u r exp 



( 5 . 4 ) 


Suemoto and coworkers obtain a value of 3.5 meV for the activation energy (kT r ) 
based on the time dependence of PL [141]. They also obtain a value of 5000 s -1 
for R r - Porous silicon is a highly disordered sample dependent system. The singlet 
triplet splitting energy of the exciton is in the range 3-10 meV. We therefore assume 
that the activation energy is in the range 3-10 meV and this translates to the range 
25-135 K for the characteristic temperature T r . This is in consonance with the 3-17 
meV range reported recently by Kanemitsu [142]. 

As mentioned in Sec. 5.1 we hypothesize that the escape rate has a Berthelot type 
behavior 


Rhop = vb exp 


T I 
-Tb. 


( 5 . 5 ) 


Employing Eqs. (S.3-5.5) in Eq. (5.2) we obtain 


I(T) = 


Jo 

1 + 1/0 exp [(£ + ?*•)] 



where = vnh'r is the reduced frequency. The temperature dependence of the PL 
intensity is explicitly displayed in Eq. (5.6). I his expression also embodies informa- 
tion on pressure, emission energy and related dependencies which we shall examine 

in the present work. 
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Figure 5.1: Possible model of tunneling through a vibrating barrier of width S and 
height Vo. 

5.2.2 Dynamics of Carrier Hopping 

The carriers can recombine radiatively or nonradiatively, or escape recombination via 
a hopping mechanism. Our hypothesis is that the hopping escape occurs as a result 
of tunneling across a vibrating barrier (see Fig. 5.1). 

The transfer probability of thermally assisted tunneling has been obtained earlier 
by Tredgold [152] and Hurd [149]. We briefly recapitulate the salient features of their 
derivation. The tunneling probability T of a carrier of mass m across a stationary 
barrier of width S and height V and in the quantum state having energy AE U (see 
Fig. 5.1) is 


T = ^oexp(-2a;S , ) 


(5.7) 


where a 1 is the extent of the carrier wave function and is given by 


a 


2 _ 2m(V - A E u ) 
h 2 


(5.8) 


For a harmonically vibrating barrier the thermally averaged tunneling rate <r) is 
determined by two competing scales, the thermal energy kT and the harmonic energy 
MH 2 /2ol 2 where M is the mass of the vibrating ion and il its frequency. Thus 

'2a 2 fcTl 


(T) = r 0 exp(-2o;5) exp 


Mil 2 


(5.9) 


Here, T 0 is a characteristic tunneling rate. A more rigorous derivation by Hurd [149] 
yields 


(r) = r 0 exp(— 2aS) exp 


ha 2 


LMfitanh(^)J 


(5.10) 
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Note that in the classical limit Ml < kT, Eq. (5.10) reduces to Eq. (5.9) 

We identify the escape rate R ho? with (r). Thus, 

Rhop = exp (5.11) 

with 


T b 


r 0 exp(-2o;5) 
Mil 2 
2 a 2 k 


(5.12) 

(5.13) 


where u q and Tq define the characteristic Berthelot frequency and temperature of 
the sample. 

A number of assumptions are made in the above derivation. The Born- 
Oppenheimer approximation is employed. A single frequency fl is assumed for the 
vibration which implies that the escape process is probably linked to a defect related 
phenomena. The microscopic nature of the phenomena, whether it is vacancy or void 
or surface related, can be speculated on, but we shall at present refrain from doing so. 
Earlier studies of carrier tunneling across barriers exist [58,141,153]. These studies 
do not consider vibrations as we have done above. 

Before we proceed further we examine the possible range of numerical values of 
the above-mentioned quantities. The mass M maybe associated with the silicon mass 
[149] and is 4.7 x 10“ 26 kg. The wave function extent a~ l is taken to be 3-7 A as 
argued by Deresmes and coworkers [144] in the context of electrical behavior of porous 
silicon. The earlier work of Hurd [149] in Si0 2 assigns a -1 the value of 3.4 A. We 
assume a to be in the range 0.2-0. 8 A -1 . The frequency of the vibrating barrier f2 
is assumed [144,149] to be 10 u - I0 12 rad/s . These considerations yield the force 
constant Mil 2 to be in the range 3 x 10~ 5 to 3 x 10 -3 eV/A 2 and hence the Berthelot 
temperature Tb (= Mil 2 /2a 2 k ) to be in the approximate range 10-200K. This is 
similar to the range suggested by conductivity experiments [143,144], The latter is 
also the quantity most directly accessible from experiments. Later, we will modify 
these assignments based on comparisons with actual experiments. For example, the 
results of Sec. 5.3.1 suggest MQ, 2 to be in the higher range [3.5 x 10~ 4 — 1.65 x 10 -2 ] 
eV/A 2 . 
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5.3 Photoluminescence Intensity 

5.3.1 Temperature Dependence 

The explicit temperature dependence for the PL intensity is given in Eq. (5.6), namely, 

I(y ^ l + ^o exp [(^ + ^)] 

The asymptotic limits of the above expression are 

7(T)r->o — * 0 (as exp[-T r /T]) (5.14) 

J(T)t- oo — ► o (as exp [-T/T B ]) (5.15) 

It is clear from the above expressions that the PL intensity will have a maximum at 
some intermediate temperature T m . This can easily be obtained to be the geometric 
mean of the characteristic temperatures 

T m = sJtVTb (5.16) 

The possible range of the characteristic temperatures have been discussed earlier. 
Based on these a priori assignment of ranges the maxima for a PL emission line 
T m ~100K. 

A normalized PL spectrum based on equation (5.6) is depicted in Fig. 5.2. For the 
solid line we had assumed T r = 75 K and T a — 100K. The maximum temperature 
T m = 87K. Notice that the intensity drops gradually for T > T m and somewhat 
steeply for T < T m . The latter behavior is on account of the strong exponential 
character at low temperatures (~ exp[— T r /T\). Also depicted in the figure is the 
experimental plot of Kanemitsu and coworkers for the PL intensity variation with 
temperature at 750 nra and for crystallite size of roughly 20 A in diameter. The 
agreement is good. 

The variation of PL intensity with temperature for crystallite of sizes 35 A and 90 
A also shows a similar behavior, with maxima around 80K. Intensity maxima have 
been reported at lower T by other workers. For example, Banerjee has reported that 
the PL intensity for a PS sample at 810 mn shows a maximum at 31K [154]. The 
dashed line in Fig. 5.2 depicts the theoretical curve with T r - 20 K , T B = 50 K- 
The agreement with Banerjee’s result is encouraging. Banerjee also reported studies 



5.3 Photoluminescence Intensity 


83 



100 

80 

60 

40 

20 


Figure 5.2: The variation of the luminescence intensity with temperature. The boxes 
denote data from Ref. [14] and crosses, data from Ref. [154]. The solid (dashed) 
line shows a theoretical plot based on Eq. (5.C) with v T — 9000 (6000) s ~ l , = 
5000(1000)s- 1 ,7; = 75(20)K and T B = 100(50)K. 

on emission lines of lower wavelengths which exhibited well defined maxima around 
125K. 

We note that similar temperature dependence of the PL line consistent with our 
model has been reported by a large number of workers. Xu, Gal and Gross [155] 
report that the PL from high porosity samples exhibited a maximum at 150-200K, 
whereas for low porosity samples, the PL intensity seemed to be constant below 
80K. Mauckner and coworkers [156] also report a maxima at 150K for the 684 nm 
emission line from a p-type sample. Sucmoto, Tanaka, and Nakajima [141] report an 
intensity maximum for the 850 nm line around 80K. Their PL intensity falls slowly 
below 80 K and by a factor of 40 from 80K to 300K. A similar behavior with T m near 
100K has also been reported for n-type samples [157]. We also note that there are 
reports of a monotonic increase in PL intensity [158] from 2K to 200K. In studies 
on partially oxidized porous silicon samples, Tsybeskov and coworkers [159] found a 
steady increase in the PL intensity for the orange band from 15 - 300 K indicating that 
T m could be over 300 K in this case. Some other groups [31,160,161] have observed 
that the PL intensity remains nearly constant for low values of T, and decreases 
above T m . This behavior is similar to that observed by Xu, Gal and Gross [155] for 
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their low porosity samples. We examine the expression for the intensity in the limit 
of low T to see that this behavior can be attributed to a low activation energy (fcX r ) 
associated with the radiative process. 

I(T) |r->o ~ — exp(“) 

Uq 1 

The T m is thus low in these samples. 

Some workers have chosen to plot the PL peak intensity with temperature [162] 
which yields a maximum at 210K and a sharp fall for T > 210 K. Zheng, Wang 
and Chen [60] report maxima in the PL peak intensity between 100-200K for various 
p-type samples. A majority of the samples show a sharper dip in intensity below T m 
than at T > T m . The PL peak position is itself a function of temperature and we 
shall address this issue later in this chapter. 

We also note that in most cases T r < Tb ■ For T > Tb the expression for the PL 
intensity (Eq. (5.6)) reduces to 



k_ 

1 4- v o exp 



This type of behavior, e.g. 


j o 

m 


— 1 oc exp 


\ T ] 

■T b . 


(5.17) 


(5.18) 


has been reported earlier in chalcogenide glasses, amorphous silicon [163], siloxene 
[164], and porous silicon [21,165]. From several fits to experimental data we also 
obtain a revised range of Mti 2 to be 3.5 xlO -4 - 1.65 x 10~ 2 eV/A 2 . 


5.3.2 The Luminescence Decay Time 

The luminescence decay time is given by equations (5.3-5. 5) to be 

1 

T = — 

-ft<r "4” R hop 



v r exp(-Zk) + V B exp(^) 

In the low temperature limit, the lifetime approaches a constant value, 

r -* — (T -f 0) 

1 


(5.19) 


(5,20) 
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Figure 5.3: The temperature dependence of the luminescence decay time. Shown are 
the experimental values obtained from Ref. [156] (crosses) and Ref. [157] (boxes). 
The solid (dashed) line represents a theoretical fit based on Eq. (5.19) with v r = 
1 5000 (35000) s -1 . v B = 100(700)s~\ T r = 125(125)K and T B = 150(150)K. 

For the small but finite temperature limit where T < T rt we can approximate 

Eq. (5.19) as 

1 


^exp(^) 



The lifetime is thus seen to decrease as T is increased from T - 0. A theoretical 
curve based on Eq. (5.19) is plotted in Fig. 5.3. This depicts a log(n) vs. log(T) 
plot which compares well with the experimentally reported values (boxes) for the < 00 ; 

nm line due to Finkbeiner and Weber [157] (their figure 3). The lifetime r ~ 10 
s at T = 10K. As the temperature is increased to ~ 200K, it falls by two orders of 
magnitude, and then levels off near room temperature. The low temperature value : 
of the luminescence decay time (Eq. (5.20)) serves to fix the value of the hopping ; 
escape rate v B - lO 2 *" 1 - The values of the other parameters employed in computing ; 
the theoretical curve are: v T — 15000s l , T r = 125K and T B — 150K. j 

The detailed behavior of the plot in Fig. 5.3 can be understood from the expression j 
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for the slope of the log r vs. log T curve 


d(logr) 

d(logT) 


T r 


T 


[l + u 0 exp + ft)] T [l + 1 ex P ( r + t b )] Tb 


(5.21) 


where v 0 = ^ ~ 10~ 2 - 10"" 3 . In the intermediate range y ft « 1, therefore the 
second term'on the right hand side of Eq. (5.21) is negligible. We can approximate 
the slope in this range to be 


d(logr) Tr 

d(log T) ~ [l + i/ 0 exp + ft)] T 


Comparing the slope at T = T r 


d[ logr) 
d(log T) 


T-T r l 1 + ^ eX P i 1 + ft) 


and at T — Tb 


d(logr) 


d(logT) 


T r 


t=t b [l 4- vo exp (ft + l)] T b 


We find that the absolute value of the slope decreases as we increase the temperature. 
This is clearly discernible in Fig. 5.3, where we observe a flattening of the curve around 
T ~ 150K. Note that the slope always remains negative. At very large temperatures 
(T > T b ), the second term of Eq. (5.21) is again dominant, i.e.. 


djlogr) ^T 

dQogT) T b 

Thus the lifetime falls rapidly to zero for temperatures over 400K. This range is not 
of much experimental interest. 

A similar dependence of the luminescence decay time on temperature has been 
reported by other workers. Pavesi and Ceschini [140] observe lifetimes to be indepen- 
dent of the excitation energy. Their measurements of the luminescence decay time 
vary over an order of magnitude of time from 10K to 300K. Mauckner and coworkers 
[156] measure the lifetime at 689.3 nm. Eq. (5.19) yields an excellent fit to their data. 
This is depicted in Fig. 5.3. They find r to vary from 10 -3 s at 25K to ~ 5 x 10"“*s 
at 200K and seek to explain this as being caused by the excitonic exchange splitting 
proposed earlier by Calcott and coworkers [166]. Amato. Di Francia and Menna [160] 
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have carried out luminescence decay studies on samples of varying porosity. They 
observe a larger decrease (~ two orders of magnitude) for 80% porosity samples as 
compared to 60-70% porosity samples (~ 1.5 orders of magnitude) in the tempera- 
ture range 10 - 200 K. A similar decrease of one and a half orders of magnitude in 
the lifetime is reported by t’Hooft and coworkers [167]. Tsuboi, Laiho and Pavlov 
[16S.169] postulate the observed red band luminescence decay time to be due to two 
states with lifetimes r 2 and r 3 , with r 2 ~ 0.1r 3 . In agreement with our results, r 3 
increases slowly with decreasing temperature from 300 Iv upto 100 K. and then in- 
creases rapidly. Our formalism can successfully reproduce these experimental results, 
working within a numerical range of parameters that has bee’n a priori fixed. 


5.3.3 Pressure Dependence 


Our formalism enables us to gauge the pressure dependence of the PL intensity. We 
once again write down our expression for the PL intensity (Eq. (5.6)) 

h 


I(T) 


1 + Vq exp {y + jr-} 


we recall that 


Mb 


Vr 

£0 

is r 


exp(— 2 aS) 


(5.22) ! 


where we have used Eq. 5.12. The application of pressure will modify the barrier : 
width S ; 


S -4 5(1 - A'P) 


(5.23) | 


where K is the isothermal compressibility. Hence from equations (5.6), (5.22) and 
(5.23), the pressure dependence of intensity can be explicitly stated. 

I{F) = 1 + [Sjejcp {F + exp{-2aS(l - KP )} J ( °' 24) 

It is clear that the intensity will drop with pressure. At room temperature, one 
can ignore the first term (unity) in the denominator of Eq. (5.24) and write 

I(P) ~ I 0 C(T)e~ 2aKSP 
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Figure 5.4: The variation of luminescence intensity with pressure, computed from 
Eq. ( 5 . 24 ) with = 1000 ,T r = 75K,T B = 100/C a = G.25A“ ! .S = 20Aand I\ = 1 
Kbar -1 . The inset shows a fit to the data reported by Zhou and coworkers (Ref. [170]). 

where C(T) contains the temperature dependent factors. Thus we see that 

logI(P)~-P 

The intensity falls exponentially with increasing pressure. In Fig. 5.4 we have 
plotted the pressure dependence of the intensity based on Eq. (5.24). Because of 
the uncertainties in the value of the compressibility for porous silicon (note that 
for silicon, K is 1.012x10 3 /Kbar) and the arbitrary scale for intensity, a detailed 
quantitative comparison with experiments cannot be attempted. However, the trend 
predicted by our formalism is experimentally observed. 

Zhou and coworkers [1/0] report that the PL intensity falls exponentially upto 20 
Kbar for two different samples fabricated from p-type silicon. This is indicated in 
the inset of Fig. 5.4 This is in agreement with the experimental results of Ookubo. 
Matsuda and ICuroda [171] who note that the intensity is strongly reduced by the 
application of pressure and is almost quenched for pressures over 40 Kbar. This 
reduction is not entirely reversible. On restoring the pressure to its original value, 
the intensity returns only to a quarter of its initial magnitude. Interestingly, thev 
note that the decrease in PL intensity due to increased pressure can originate partly 
from the increased non-radiative recombination due to an enhancement in tunneling. 




5.4 Model Dependent Results 


89 


Ryan. Wamslev and Bray [172] report a rapid quenching of luminescence upto 90 Kbar 
where the signal becomes very difficult to detect. In a recent work, however, Zeman 
and coworkers [173.174] have been able to observe PL upto 150 Kbar of pressure. 
These experimental results are in consonance with the predictions of our theory. 

5.4 Model Dependent Results 

Several microscopic models to explain the PL in porous silicon have been proposed 
(see Sec. 1.3). A popular model is the quantum confinement model which posits an 
enhanced direct gap on the basis of carrier confinement in silicon nanocrystallites. 
The results presented in sections 5.2 and 5.3 are not necessarily partial to any mi- 
croscopic model. They can be derived solely in terms of the competition between 
an Arrhenius type radiative recombination term and a Berthelot type escape term. 
It would be interesting to explore our proposed scheme further in the light of the 
existing microscopic models. In this section we undertake this task in the light of the 
quantum confinement model. 

In the quantum confinement model, the luminescence energy hu is related to the 
bulk gap of silicon E g , the confinement energy of the carrier A E u , and a subtractive 
term (~E S ) coming from the exciton binding energy, the transverse optical phonon 
energy, etc. Thus 

fiuj — Eg + Aj — E s (5.25) 

Recall from Fig. 5.1 that the escape rate is given in terms of the barrier height 
(V-AE U ). The carrier is in the level marked by A E u and we can ascribe to A E u the 
upshift due to the quantum confinement of the carrier in the nanocrystallite. From 
Sec. 5.2, the detailed form of the escape rate is (Eq. (5.5)) 

Rhop = v B exp (5.26) 

(5.27) 

2 2 m(V - A Eu) 


= To exp(— 2aS) exp 


2 a 2 kT 


Mfl 2 


where from Eq. (5.8) 
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Taking the logarithm and making only the energy related term explicit 

log(J4o P ) = 7 o - 7 i (V - A E u ) 1/2 + 72 O ' - A E u ) (5.2S) 

where 70 = logT 0 , 71 = 2y/2mS/h, and 72 = (j^) and the constants 7 are all 

positive. This establishes the framework within which we can examine the lumines- 
cence decay time dependence on emission energy, and the variation of the PL peak 
with temperature and pressure. Note that the subtractive term of Eu. (5.2,1 1 is small 
(~ 0.1 eV [59]) and serves as a correction term to the barrier height. Thus we take 

1 •’ - A E u -* V + E g - ft-c « 5.29 1 

5.4.1 Dependence of the Lifetime on Luminescence Energy 

The dependence of the lifetime on the luminescence energy has been studied by several 
workers. In our formalism this dependence occurs explicitly in the non- radiative term 
and perhaps implicitly in the prefactor I 0 term. We recall Eq. (5.191 for the lifetime 

1 

r = — — 

Rr ~f~ R- hop 

The radiative recombination rate Rr is small compared to the hopping escape rate 
Rkop at T ~ 300 K. Hence 

log(r) = -log (R r + R hop ) (5.30) 

= - log {Rho?) + O (5.31) 

Substituting in the above equation, the expression for log R hvp given in Eq. (5. 28) 

log(r) = -[70 - 7 i(P - A £ u ) 1/2 + 72 (V - A E u )) ( 5 . 32 ) 

where the dependence of the activated radiative term R r on the confinement energy, 
if any, will be ignored for the present. The behavior of log; r ! on the confinement 
energy is depicted in Fig. 5 . 5 . Recall that the hopping term is dependent on the 
confinement energy as indicated above. The plot in Fig. 5,5 is almost linear and the 
lifetime falls over an order of magnitude in the confinement energy range 0.2 to 1.0 eV 
(corresponding luminescence energy range is 1.4 to 2.2 eV). The parameters used are 
within the circumscribed ranges outlined in the previous sections. For comparison we 
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Figure 5.5: The variation of the luminescence decay time r with confinement energy 
A E u . The experimental data (boxes) is from Ref. [146]. The solid line depicts 
a theoretical plot based on Eq. (5.32) with R. = 4000s -1 , T 0 = 5 x 10 8 s -1 ,7 1 = 
8(eV) -1/2 . 72 = 0.4(eV) -1 and V = 2.5eV. 

plot the experimentally reported lifetime data due to Pavesi and coworkers [146] for a 
S7% porosity sample at room temperature. They extracted the lifetime from a fit of 
the time dependent PL intensity to a stretched exponential form as given in Eq. (5.1). 
The uncertainties inherent in extracting the lifetime inhibit a detailed quantitative 
comparison. Nevertheless the agreement between theory and experiment is excellent. 

Past studies of the dependence of the lifetime on the energy reveal uncertainties in 
assigning a functional relationship. t’Hooft and coworkers [167] postulated a power 
law dependence r ~ AE~ 3 ^ 2 . We have been able to replot their data along lines 
similar to Fig. 5.5. On the basis of effective mass theory, Hybertsen has proposed 
[132] that r ~ A E~ 3 . He has also suggested that for phonon dominated processes, 
r ~ A /7-3/2. Q n t j ie 0 th er hand, several experimental groups have proposed that 
the inverse of the lifetime depends exponentially on the luminescence energy (lnr = 
-£huj). This has been the observation of Vial and colleagues [58] who obtained £ = 
3.57 eV -1 , Ookubo, Hamada and Sawada [148] (their Fig. 2) and Kanemitsu [4] who 
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observed a linear variation of ln(T _1 ) with the emission energy at high temperatures 
(f = 3.4 eV -1 ). The data of Pavesi and coworkers depicted in Fig. (5.5) has also 
been modeled by a similar expression with £ =1.1 eV 1 in a later work [1/5]. Since 
hu ~ E g + A E u , the reports from these three distinct experimental groups are in 
conformity with the our model. 

The influence of the confinement energy on the lifetime can be understood in a 
simple way. Since A E u < V we Taylor expand Eq. (5.32) to obtain 


log(r) = -[tq - 7 iV 1/2 + j 2 V] - A E u 


2j 2 W 


7i 


2\V 


ChAE 2 ) (5.33) 


From the numerical values quoted earlier 72 -C 7i, and the second term in square 
brackets on the r.h.s. is positive. Hence, using Eq. (5.29), 

log(r) ~ —hui — O(AEl) 


This is the behavior observed in Fig. 5.5 and also in the experiments discussed above. 


5.4.2 Temperature Dependence of the PL Peak 


In Sec. 5.3.1 we described the temperature dependence of the PL intensity. Several 
workers have also mapped the temperature dependence of the peak PL emission 
energy hui p . There are reports of a blueshift with perhaps a subsequent redshift with 
temperature. Our model provides some insight into this anomalous behavior of the 
peak emission energy. 

Recall our expression for the PL intensity in terms of the tunneling parameters a 
and 5 (equations (5.4), (5.6), (5.12) and (5.13)), 


7(T) = 


Io 


1 + S e xP (ft ) exp(-2aS) exp 
The emission energy is contained in the parameter a. Recall that 

2m (V - A E u ) 


(5.34) 


a 2 = 


_ 2m(V 4- Eg - hu p ) 

r? 

where we have absorbed the small subtractive term due to exciton binding energy, 
phonons etc. into V. The peak energy is located by differentiating Eq. (5.34) with 
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r.espect to a, 


SMQ? 

a P ~ 2kT (5.35) 

~ f (5.36) 

Combining the above expressions we obtain the PL peak emission energy 

hu p = V + E g -^ (5.37) 

where Co is a constant 

i 

_ h 2 ( SMQ 2 \ 2 
° 2 “ 2m { 2k ) 

The above expression for hu> p implies a blue shift of the peak energy with tem- 
perature. It is known that the band gap of silicon decreases with temperature 
Eg = E° — CiT where C\ = 2.3 x 10 -4 eV/K. Using this as a guide and working 
within the quantum confinement model as noted earlier, we have 


hu, p = V + iJ-cT-l I 

= Co — CiT — ~ (5.38) 

Clearly this expression implies an initial blue shift followed by a red shift. The 
crossover is determined by T = ( 2 c 2 /c{) 5. Depending on the value of this cross over 
temperature, the PL peak may exhibit an initial blue shift or a red shift. In other! 
words the anomalous dependence of the peak PL emission energy is possible. 

This discussion must be viewed as a qualitative way of understanding the 1 
anomalous temperature dependence of the peak energy. The barrier frequency 
(f> ~ 10 11 — 10 12 sec -1 ) occurs as the fourth power in the coefficient C 2 . There is 
a large uncertainty in the parameter S. Besides, all parameters may be affected by 
post-anodization treatments. The temperature coefficient of the band gap in a silicon 
quantum wire may not necessarily be the same as of bulk silicon. Some uncertainty ini 
the quantitative determination of the crossover temperature is thus inevitable. The 
crossover temperature may range between 20 to 1000K depending on the value of cfci 
Eq. (5.38) may be taken as encapsulating the correct physical behavior. The term 
C 2 /C 1 can be obtained from experimental fit to data and can be used as a character- 
izing parameter. ! 
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An anomalous temperature dependence of the peak position, similar to that pre- 
dicted by Eq. (5.38) has been reported in earlier experimental studies. Zheng. Wang 
and Chen [60] perform temperature dependent studies of the PL peak from four 
different samples of porous silicon, prepared with different current densities and an- 
odization times. They observe an initial blue shift in the PL peak in three of their 
four samples as the temperature is increased from 20 K. One of the samples show a 
consistent red shift with increasing temperature. Xu. Gal and Gross [I55j also observe 
both positive and negative shifts of the PL peak with temperature. Xarasimhan and 
coworkers have reported an increase in the PL peak energy till 310 K. The PL spec- 
tra reported by Ivesan and coworkers [158] also exhibit an anomalous behavior in the 
temperature dependence of the peak position. Shimuzu and eoworkers have reported 
a steady blue shift for one sample and a steady red shift for two ot hersamples [161]. 
Based on experiments carried out on silicon nanocrystallites Kanemitsu [142] reports 
that the PL peak energy is almost temperature independent and exhibits only a small 
red shift from 100K to 300K. 

5.4.3 Pressure Dependence of the PL Peak 

A number of studies of the pressure dependence of the peak PL emission energy hu> p 
have been conducted. Some workers have found that hu> p exhibits a small blueshift 
with applied pressure and then redshifts significantly [63,170.172,176]. Others have 
found only a red shift. These anomalous reports can be understood on the basis of 
our model. 

In Sec. 5.4.2 we have obtained an expression for the peak energy (Eq. (5.38)) 

>^ P = V + E,-^ 

Mote that c 2 contains information on the barrier width S (e.g. c 2 oc S'-). The 
application of pressure will modify S as discussed in Sec. 5.3.3 (Eq. (5,23)). It will 
also modify the gap E g . The indirect gap in crystalline silicon redshifts with pressure 
and with a pressure coefficient -1.41 to -1.5 meV/Kbar [177,178]. In amorphous silicon 
a larger redshift has been reported [179] with pressure coefficient ~2±0.5 meV/Kbar. 
Theoretical calculations by Yeh. Zhang, and Zunger [180] on silicon quantum wires 
also suggest a redshift but with a pressure coefficient smaller than in crvstalline 
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silicon. Thus 

E, = E° - r,P 


where r\ may range from 0.1 to 2 meY/Kbar. The peak energy will then display a 
pressure dependence of 



(V+E°-^) + P 

Go *f" CL\P — CLiP' 



P 2 K 2 c 2 

J' 2 


(5.39) 


Working within the numerical range of values chosen earlier, one can see that d, 
the co-efficient of the term linear in P can assume both positive and negative values. 
This in turn depends on the specific surface conditions. If ai is positive, the r.h.s of 
Eq. (5.39) has a maximum and an initial blueshift can thus be explained. 

Experimental studies of the pressure dependence of the PL peak are strongly 
dependent on the hydrostatic pressure transmitting medium employed. The ex- 
periments conducted in alcohol [63.170.172.176,181] show an initial blueshift upto 
pressures of ~ 20 Kbar. beyond which there is a significant redshift. In another 
work employing a fluorocarbon medium. Ookubo, Matsuda and Kuroda [171] report 
redshifts above 20 Kbar. but are uncertain about the PL peak behavior at lower pres- 
sures. Sood, Jayaraman and Muthu [62] report a redshift at all pressures when the 
medium is paraffin oil. Their peak energy vs. pressure plot appears to be parabolic, 
which is again consistent with Eq. (5.39). 

Cheong and coworkers [63] argue that due to these inconsistencies, pressure exper- 
iments could not be used as critical tests of the quantum confinement model. They 
conducted pressure studies using both alcohol and helium as the transmitting media. 
The PL peak shift in alcohol was consistent with the earlier results, where a blue 
shift was observed upto ~ 25 Kbar. When helium was used instead of alcohol, only 
a red shift was seen. Based on this medium (and therefore, surface) dependence of 
photoluminescence, they asserted that the results are inconsistent with the standard 
quantum confinement model. In our framework, these results can be explained within 
a quantum confinement scheme. 
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5.5 Discussion 

We have proposed a simple theoretical framework to understand and systematize 
the behavior of the dominant red band in the photoluminescence of porous silicon 
and silicon nanocrystallites. Our basic model as outlined in Sec. 5.2 and 5.3 uses 
a limited set par am eters {v r , i/ B ,T r ,T B }, each of which can be ascertained in an 
independent, a priori fashion. These values can further be confirmed and refined by 
comparing our predictions with a variety of experiments. We do not take recourse to 
computer modeling or simulations. Our framework, which employs only analytical 
manipulations, is physically transparent. 

The photoluminescence in porous silicon has been seen to be extremely sensitive 
to post-anodization surface treatments. This has led many workers to hypothesize 
a radiative recombination process involving intermediate surface states. Our model 
explains the surface dependence by hypothesizing an escape mechanism involving 
the surrounding environment. The parameters M , Q. and V are surface dependent. 
The fact that the observed variations in photoluminescence with temperature and 
pressure can be reproduced by our model supports this assumption. 

Porous silicon is a disordered system and its properties are dependent on sample 
preparation conditions and on surface chemistry. This is reflected in our model. 
We have assumed in Sec. 5.2.1 that the time dependence of the PL intensity is of 
the Kohlrausch type exp[-(t/r) p ] with 0 < p < 1. The parameters {v r , v B ,T r ,T B } 
assume a range of values reflecting the disordered nature of the sample. In Sec. 5.4.3 
we have noted that the PL peak energy has an initial blue shift when the hydrostatic 
pressure transmitting medium is alcohol. If - as noted recently [63] - the medium is 
helium, the peak is uniformly redshifted. In our formalism this is accounted for by 
the sign of the linear coefficient of the peak (Eq. (5.39)). 

A key hypothesis in our study is the presence of a Berthelot component in the PL 
intensity. As noted earlier, some previous studies on a-Si, siloxene. and porous silicon 
have reported an exp[-T/T 0 ] dependence of the PL intensity I(T). These studies do 
not ascribe this dependence to tunneling of carriers across a thermally vibrating 
barrier. We have shown in Sec. 5.3.1 that their expression for If T) maybe viewed 
as a limiting case of our formalism. Some studies on the conductivity of porous 
silicon [137,144] do explicitly mention a Berthelot process. An issue of relevance 
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then, is the nature of the barrier. The phenomenological models fpr PL behavior 
by Vial et al. [5Sj. Sacillotti ei al. [182. 1S3] and Qin and Jia [184] do describe non- 
radiative tunneling of carriers to a SiO- coating surrounding the crystallite. In a 
recent work Schuppler et al. [57] have claimed that PL occurs in nanocrystallites of 
sizes ~ 15 A. Our barrier width S has a range 5-20.4. It is tempting then to think of 
the carriers in Fig. 5.1 to be hopping from the core to the surface of these luminescing 
nanocrystallites. At this stage we can only speculate on the nature of the barrier. Our 
results in Sec. 5.3.1 and Sec. 5.3.2 hold if we assume an Arrhenius and a competitive 
Berthelot term irrespective of the microscopic origin of the latter. We refrain from 
ascribing a definite microscopic material character to the barrier at this stage. 

We have obtained the tunneling probability (T) in terms of the barrier parameters. 
We had noted that a more rigorous derivation by Hurd [149] yielded 


(P) = r 0 exp(— 2aS) exp 


ha 2 


|_AfO tanh(44F) 


Even for a barrier frequencies as large as 0 ~ 10 12 sec -1 , as suggested by Hurd, the 
approximation t&nh(hQ./2kT) ~ hCljlkT is excellent, upto temperatures as low as 
1QK. Hence, this approximation is valid in our case. 

The barrier width S has been assumed to be temperature independent. The 
linear thermal expansion coefficient 3 for silicon has been reported to be 1 - 5 x 
10 -6 K -1 . Thus the thermal expansion S3 AT of the barrier is less than 0.1 c 7c when 
a temperature variation of 100K is assumed. Hence the thermal expansion of the 
barrier can be safely ignored. 

The radiative recombination has been assumed to follow an activated process 
R- = u T exp(— T t /T). This is in agreement with the studies of Suemoto and coworkers ' 
[141] who experimentally obtain a fit to Rr, 

-E a - 


R~ = .4 exp 


kT J 


+ B 


where .4 = 5000s -1 and B = ISOs -1 . Since B -C -4, we have left out this temperature 
independent term, for the sake of mathematical simplicity. If B is included in the 
model, the low temperature limit of the luminescence decay time (Eq. (5.20)) will 
become 
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PL studies in amorphous silicon (a-Si) have been conducted for the past two 
decades [185-188]. Peak PL emission energies have been reported between 1.0 - 
1.5 eV. The mechanisms contributing to PL have been debated ranging from donor- 
acceptor recombination to self-trapped excitonic models. Disorder plays an important 
role in these models. We have not undertaken a detailed survey of a-Si results yet. 
However, we point out two features of these studies. Street [187] has reported a 
maximum in the PL efficiency at 50K. He has attributed the fall in the PL efficiency 
below 50K to Auger recombination of geminate electron-hole pairs. There has been 
a recent suggestion by Delerue and coworkers [189] that the non-radiative process in 
silicon nanocrystallites is dominated by fast, Auger recombination. Our work however 
does not explicitly invoke an Auger mechanism. Another feature of" a-Si studies 
[185.186] is to ascribe the competitive non-radiative process to tunneling to defects. 
This is in consonance with our model. 

Since porous silicon is a highly disordered system some studies have invoked Mott's 
variable range hopping model [4,190]. It has been customary to analyze the temper- 
ature dependence of the lifetime in terms of an activated or hopping process. 


1 

- oc exp 
r 



al 


Here a = 1 is a pure Arrhenius dependence while a fractional value of a indicates a 
Mott type variable range hopping process. Kanemitsu [4] has suggested a value of 
a = l/3 corresponding to two-dimensional exciton hopping in the interface region of 
the silicon nanocrvstallites. Our model attributes the temperature dependence of the 
lifetime to a mathematically well-defined competition between radiative and escape 
mechanisms where the latter is determined by tunneling across thermally vibrating 
barriers. An attempt by Mares [137] to fit the conductivity data to a Mott type 
process was unsuccessful in as much as it resulted in hopping distance of a thousand 
meters. 

Based on our model we can suggest a range of experiments. The pressure depen- 
dence of the lifetime can be examined. Eq. (5.3) and Eqs. (5.10-5.12) yield 

i = T 0 exp(— 2aS] exp[T/T s ] + R r 


From Eq. (5.23) S t S(1-A7P ) with pressure P. This suggests that the lifetime falls 
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almost linearly with increasing pressure. 

log(r) ~ —(2aKS)P 

The conductivity of the sample has been observed to follow a Berthelot type 
law as noted above. The pressure dependence of the conductivity term can also be 
examined. 


a = exp[— 2aS'(l — KP))a\ 

where temperature and related dependencies are present in o\. Thus the conductivity 
should increase with pressure. The role of the pressure transmitting medium in 
governing the extent to which the lifetime and conductivity maybe significant as 
noted in the pressure studies in Sec. 5.4.3. 

In summary, we have presented a transparent analytical framework to explain a 
diverse range of experimental studies on the PL behavior in porous silicon. This, we 
hope will provide a unified perspective to workers in the field. Systematic studies to 
verify, sharpen, and critique our model are possible. We hope to extend this frame- 
work to other aspects of porous silicon such as electroluminescence and to related 
nanocrvstalline and disordered semiconductors. 


Key Reference 

The major reference for the material contained in this chapter is: 

• George C. John and Vijay A. Singh. Phys. Rev. B 54, 4416-4419 (1996). 
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a proper characterization of the porous silicon surface. In this section, we discuss 
various efforts employing infrared vibrational spectroscopy to validate these models. 

The formation process of porous silicon involves Si, HF and an organic medium 
such as C 0 H 5 OH. Consequently, the surface composition is largely unpredictable. 
Even the hydrogen passivated porous silicon surface is unstable and susceptible to 
contamination by other elements. Visible PL from porous silicon has been found to be 
extremely sensitive to surface treatments. This has prompted the school of thought 
which attributes the luminescence from porous silicon to the presence of luminescing 
chemical complexes on the large surface. Several molecules have been proposed as 
possible sources of PL. 

The contention that certain siloxene derivatives constitute the source of the vis- 
ible room temperature luminescence in porous silicon [193] was made largely on the 
basis of a comparison between the infrared transmission spectra of porous silicon and 
siloxene, which appeared to be remarkably similar. Further, as mentioned earlier in 
chapter 5, the temperature dependence of PL intensity in both the materials is well 
described by the relation 

* — & 

The large width of the PL spectra of porous silicon (see discussion in chapter 4) 
also supported the siloxene case, as all reported spectra fell well within the siloxene 
luminescence energy range. 

Tischler and Collins showed that the infrared transmission spectra of freshly pre- 
pared porous silicon samples do not show any oxygen related modes [22]. At the same 
time, these samples also exhibit efficient visible luminescence. The pronounced lack 
of oxygen indicates the absence of siloxene or other oxygen containing complexes on 
the surface. On oxidation, porous silicon spectra was seen to develop oxygen related 
modes and become similar to the siloxene infrared spectrum. They concluded that 
the presence of siloxene is not necessary for obtaining efficient PL. Thus the refutation 
of the siloxene model was also premised on careful infrared studies. 

In a later study [194], they observed that during UV illumination of the PS sample, 
hydrogen desorption occurs (as illustrated by a drop in intensity of the Si-H vibration 
band) along with photoluminescence degradation. It was also observed that thermal 
flnnPflKng causes degradation of PL due to hydrogen desorption. This led to the belief 
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that surface Si-H n molecules (silanes) were the luminescing agents in porous silicon. 
However, in another infrared study [195] carried out with various etching times of 
silicon in HF, it was seen that the intensity of the PL spectra bore no correlation with 
the Si-H vibration mode. The observed degradation of PL on hydrogen desorption 
could be attributed to an increase in the number of dangling bonds which can act as 
centers for non-radiative recombination. 

To account for the sensitivity of PL to post-anodization surface treatments within 
a quantum confinement model, a hybrid model invoking the recombination of carriers 
confined in the nanocrystallites via the surface states has been suggested. It has 
been proposed that various chemical ligands such as H. OH. O and other radicals 
can saturate the Si cluster. This results in a wide spectrum of disorder induced 
trap states which act as mediators in the light emitting process. A convenient way 
to characterize these surface states is provided by infrared vibrational spectroscopy. 
This can yield information on the electron states as well as on the vibrational modes 
of the surface spectra. 

The above discussion underscores the importance of impurity related studies in 
discriminating between the surface related models for PL in porous silicon. It would 
be useful to develop a phenomenological scheme for the assignment of the vibrational 
frequencies of impurities in porous silicon. In this chapter, we propose a simple model 
for the stretching mode vibrational frequencies of hydrogen and chlorine in varying 
environments found on the silicon crystallite surface. The model is calibrated by 
appealing to well established vibrational spectra assignments in crystalline silicon. 

It has been known for sometime that the stretching mode vibrational frequencies 
of Si-H (i/si-//) in substituted silane molecules SiHR^Rs [196]. amorphous solids 
such as a-Si and a-Si0 2 [197], and in crystalline silicon (c-Si) [198] correlate with the 
electronegativities of {RiRoRs}. The R* here may denote an element or a radical. The 
following explanation is usually proffered for this correlation. As the electronegativity 
of the substituting species {R1R0R3} increases, the s-character of the Si-H bond 
increases. A calculation of the s-component of the bond-order matrix substantiates 
this [199]. Because of the enhanced s-character. the Si-H distance dsi~n decreases 
and the effective force-constant increases. Thus vsi-H increases. In this chapter, we 
develop a phenomenological approach to predicting the vibrational frequencies of the 
Si-H bonds depending on the immediate molecular environment. Section 6.2 outlines 
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this- model. In section 6.3. we predict the frequencies of possible surface environments 
in porous silicon and critique the existing assignments in literature. 


6.2 The Renormalized Electronegativity (REEL) 
Model 

6.2.1 Electronegativity and Orbital Radius 


To motivate the application of the orbital-radii scales to the systematization and 
prediction of the vibrational frequencies of H and other impurities in semiconductors, 
we shall first relate the electronegativity to the orbital radius of the element Rj. The 
three orbital radii listed by Zunger [200] are r a , r p . and v d which are respectively the 
crossing points of the first principles pseudopotentials for l = 0, 1, and 2. The role 
of the a character of a vsi-H bond in determining the Si-H vibrational frequencies 
has been pointed out by other workers [198,199]. Further, vsi-H has been related 
to Sanderson electronegativity SR of element R 3 [197], This motivated us to seek a 
relation between the electronegativity and the Zunger orbital radii r s (Rj ) since the 
s-character of the Si-H bond appears to play a key role [199]. 

Pauling defined electronegativity as "the power of an atom in a molecule to attract 
electrons to itself" [201]. Several numerical scales have been devised for a quantitative 
representation of the electronegativity scale. Allred and Rochow chose to represent 
electronegativity by the electrostatic force exerted on the valence electron by the 
nucleus of the atom [202]. The Allred - Rochow electronegativity \ AR is expressible 
on the Pauling scale (i.e., with values close to the Pauling electronegativity xp values 
in magnitude) as 


Xar = 


0.3592Z 2 

Rl 


+ 0.744 


( 6 . 1 ) 


where Z is the effective nuclear charge, and R c is the single bond covalent radius. 

Another similar approach is due to R.T.Sanderson who defines a stability ratio 
(SR) [ 203 , 204 ]. This is the ratio of the average electron density (p e ) of the atom to 
that of an inert atom having the same number of electrons. 

3 Z 
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Figure 6.1: Row-wise plots of the Sanderson electronegativity SR(Rj) and l/r s (i? 7 ), 
where r s (Rj) is the l = 0 orbital radius of the element. The principal quantum 
number of the valence shell is denoted by n. 

whereas the electron density of the corresponding inert atom is calculated by inter- 
polating the noble gas electron densities. SR is reported to scale linearly with ^/\j> 
[197], 

^Xp = 0.21SR + 0.77 (6.2) 

We find a similar linear relationship, but with the slope 0.19 (instead of 0.21) and 
intercept 0.79 (instead of 0.77). 

Eqs. (6.1) and (6.2) suggest a linear relation between l/r s and SR of the form 

SR(Rj) a -tVt 

r s{Rj) 

can be established for each row of the periodic table [205]. This is depicted in fig- 
ure 6.1. But to establish a universal linear relationship incorporating several rows of 
the periodic table, we need to postulate another quantum scale. 

According to Zunger [200] the orbital radii are characteristic of the atomic cores 
whose defining quantum numbers are 1, 2. .... n — 1. In this spirit we attempt to 
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Figure 6.2: Linear relationship between the Sanderson electronegativity SR(Rj) and 
the REEL quantum scale V(Rj) = ^Jn(Rj) — l/r s (Rj ) where n(Rj) and r s (Rj) denote 
the principal quantum number of the valence shell and the l = 0 orbital radius of the 

element respectively. 


correlate the electronegativity SR(Rj) with the functional form f((n(Rj) - 1 ),r 5 ). 
An inspection of figure 6.1 reveals that there is a systematic shift in the slope of the 
row wise plots. This insight, and some trial and error, led us to the discovery that 
SR scales linearly with the quantity [205] 


V(Rj) = 


\/n(Rj) ~ 1 
r*( Rj) 


(6.3) 


which we call the valence shell renormalized electronegativity (REEL). 

Fig. 6.2 depicts a linear scaling behavior: 

SR(R j )=aV(R j ) + S 0 (6.4) 


Here the slope a and the intercept S 0 are 2.07 (± 0.02) a.u. and -1.66 (± 0.006) 
respectively. The uncertainty in the calculated electronegativity is ± 0.01 and is of 
the same order as the uncertainty in the original electronegativity data [206]. For 
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Table 6.1: Stretching mode Si-H frequencies ( v$i-n ) for well established environments 
in crystalline (c-Si) and amorphous (a-Si) silicon. Data for c-Si are taken from Shi et 
al. [198] and for a-Si from Ivniffler et al. [207]. The symbol 'd‘ represents a dangling 
bond. 


No. 

System 

i 'si-H (cm -1 ) 

Remarks 

1 

(SiSiSi)SiH 

1980-1990 

c-Si 

2 

(CSiSi)SiH 

2028-2030 

c-Si 

3 

(SiSi)SiH 2 

2055-2066 

c-Si 

4 

(CCSi)SiH 

2083 

c-Si 

5 

(CSi)SiH 2 

2105-2107 

c-Si 

6 

(OSi)SiH 2 

2160-2162 

c-Si 

7 

(OC)SiH 2 

2210-221S 

c-Si 

8 

(0 2 Si)SiH 

2193 

a-Si 

9 

(0 3 )SiH 

2247 

a-Si 

10 

(0 2 )SiH 2 

2219 

a-Si 

11 

(OSi 2 )SiH 

2100 

a-Si 

12 

(Sidd)SiH 

1830-1840 

c-Si 

13 

(CSid)SiH 

1925-1931 

c-Si 

14 

(Sid)SiH 2 

1957-1965 

c-Si 


larger values of V, the stability ratio saturates. Our attempts to correlate SR{Rj) 
with the orbital radius r p or (r s + r p ) were not as successful. 

6.2.2 Vibrational Spectra 

We now establish the relationship between the Si-H stretching frequencies vsi-H 
observed in silicon and the quantum scale defined by us. Table 6.1 lists some of the 
well known Si-H stretching bands observed in crystalline silicon (c-Si) and amorphous 
silicon (a-Si). The frequencies listed in the first seven rows are observed for c-Si and 
the next four for a-Si. The last three entries involve dangling bonds and are observed 
in c-Si. The values are taken from the data culled by Shi et al, [195; and Ivniffler 
et al. [207]. As demonstrated by Fig. 6.3. one obtains a linear relationship between 
V si-H and our quantum scale. 
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Figure 6.3: Linear relationship between the Si-H band stretching frequency vsi-H and 

4 

the quantum scale r )^y(R 1 ) where the summation runs over the four bonds of the 

■ 7=1 

central silicon atom. 


"Si-H = m Y, V(Rj ) + i/ 0 (6.5) 

7=1 

Here the summation is over all four nearest neighbors of the central Si atom and 
includes H. For H we have taken the relevant REEL value V to be 2.86 (a.u.) -1 . For 
the dangling bond neither a n(Rj ) nor a r s (Rj ) can be defined. To include it into 
our scheme we take V to be 0.54 (a.u.) -1 . These values have been obtained from 
empirical fit to data. In Eq. (6.5), the slope m = 61.03 cm -1 -(a.u.) and the intercept 
i/ 0 = 1456.45 cm -1 . The uncertainties in the slope and the intercept are Am = ±0.07 
cm - 1 -(a.u.) and A z/ 0 = ±7.3 cm -1 respectively. The relationship given by Eq. (6.5) 
has predictive value with the attendant uncertainty Az/ = ±(11 - 16) cm -1 . We 
shall employ it to obtain Si-H stretching frequencies in environments likely to occur 
m porous silicon. 
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6.3 Discussion 

In Sec. 6.2, we have developed a phenomenological model to compute Si-H stretching 
mode vibrational frequencies for various environments that ma\ occur within poious 
silicon. In this section, we discuss the existing assignments in literature lor the 
infrared spectrum of porous silicon as reported by different groups. V e locus on the 
Si-H stretching mode vibrations and suggest revisions of the existing assignments. 

The Si-H n vibrational modes in a purely silicon environment reported in literature 
fall in three frequency ranges (i) 600 - 835 cm -1 (bending modes), (ii) $< 0 - 91 C cm 1 
(scissors mode) and (iii) 2050 - 2150 cm -1 (stretching modes). As seen from Sec. 6.2. 
the frequencies depend largely on the immediate environment. In consonance witli 
this, the Si-H stretching modes in a Si-0 environment is shifted away from 2050 - 
2150 cm -1 and is seen to lie in the range 2150 - 2300 cm -1 . A summary of the 
assignments reported for various vibrational modes is depicted in Fig. 6.4. Besides 
Si-H vibrational frequencies, Si-0 modes are also seen in the figure. 

There exists a consensus in literature regarding the assignments of the Si-H scissors 
and bending modes. The Si-H 2 scissors mode is reported by almost all workers to 
be in the range 870 - 916 cm -1 (see Fig. 6.4). Ogata and coworkers also report the 
presence of an Si-H 3 bending mode in this range. All other workers report Si-H n (n 
= 1,2,3) bending modes in the range 600 - 835 cm -1 as depicted in Fig. 6.4. 

While there is a general agreement in the scientific community regarding the 
scissors and bending modes, there exists some disputes in literature over the detailed 
nature of the stretching mode vibrations. In particular, three peaks are observed 
in the range 2050 - 2150 cm -1 , which were initially assigned to Si-H (2087 cm -1 ), 
Si-H 2 (2110 cm -1 ) and Si-H 3 (2140 cm -1 ) [210,211], but this has been contested by 
several works on amorphous silicon [212-214] as well as on porous silicon [191]. In 
this context, the values we have cited in Table 6.1 have been chosen from the data 
that has been obtained from crystalline silicon [198]. A linear relation of the form 
given by Eq. (6.5) is seen to hold extremely well for the given set of data. 

The theoretical approaches to the evaluation of vibrational modes in silicon have 
relied on phenomenological schemes as well as ab initio electronic structure calcula- 
tions. The first attempt to systematize and predict the vibrational modes in silicon 
clusters was based on a semi-empirical electronic structure calculation due to Singh 
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Figure 6.4: The vibrational spectra of Si-H and Si-0 bonds in porous silicon reported 
by various groups: (a) Gupta et al. [191] (b) Brandt et al. [193] (c) Tischier et al. [22] 
(d) Tischier et al. [194] (e) Robinson et al. [195] (f) Schuppler et al. [57] (g) Dubin 
et al. [208] (h) Ogata et al. [209]. The figure divides the assignments made by the 
respective groups into three broad categories as shown. The detailed assignments are 
discussed in Sec. 6.3. 
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and coworkers [215]. Larer. first principles calculations were carried out. along with 
experimental studies on the Si(100) surface by Chabal and Raghavachari to study the 
stretching modes of monohydride and dihydride silicon [216.217]. The monohydride 
symmetric mode was assigned to 2099 cm -1 and the dihydride symmetric mode to 
2091 cm -1 . The calculations were limited to very small clusters, the largest duster 
being Si 9 H 14 . The calculated frequencies were also found to be dependent on the level 
of the theoretical approximation. In another study. Deak et al. [53] have considered 
siloxene structures which are planar [Si 2 H(OHi] 3n . linear [SLHA V, n . and ring-like 
[Si 6 H 6 0 3 ]„ as well as substituted ring-like siloxenes [SifiH^^OH^O A judidous 
mix of quantum chemical techniques were employed to evaluate bond -trengths. force 
constants and make infrared frequency predictions. A stretching frequency of 2105 
cm- 1 was predicted for Si-H in the planar structure. This is in consonance with 
the dominant Si-H stretching band in porous silicon. A peak in this range is also 
observed m m-situ siloxene. In a recent work. Ogata and coworkers [209] perform 
detailed calculations similar to those of Chabal and Raghavachari [216] to predict 
the stretching as well as various deformation mode frequencies for the Si-H n species 
in porous silicon. They find that the vibrational frequencies increase with increasing 
hydrogen content and assign the 2087 car 1 line to the Si-H stretching mode, 2108 
cm 1 line to Si-H 2 and the 2142 cm -1 line to Si-H 3 . 

The earlier phenomenological approaches to the prediction of infrared modes in 

silicon are due to Lucovsky [197] and Sahu and coworkers [199]. These are similar in 

.spirit to the present work. Lucovsky employs three independent linear relations akin 
to Eq. (6.5) 

vsi- Hn = ™ n 52 SR(R_ 7 ) + u 0n n = 1, 2, 3 (6.6) 

where R, are the various chemical ligands attached to silicon (other than hydrogen), 

,-) the respective Sanderson electronegativity and m„ and the characteristic 
constants of the Unear relation. In a later work. Sahu and coworkers [199] replace 
Eqs. (6.6) by a single linear relation: 

"Si-H = 39.7 52 SR(Rj) - 1512.7 
i= i 

Note that the sum is over all nearest neighbors including hydrogen. Further, the 
erS0 ” electr0ne ? at ™j- Of hydrogen is reassigned to 4.25 from 3.55. However, this 
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fit does not work well at high frequencies as there is a pronounced saturation effect 
in high electronegativity environments. In comparison, the linear relation established 
between REEL and the Si-H vibrational modes (Eq. 6.5) is seen to hold good for a 
wide range of frequencies (see Fig. 6.3). 

The present study shows that the vibrational frequency increases with increasing 
hydrogen content (Eq. (6.5)). This is also supported by most of the theoretical 
calculations described earlier. A notable deviation from this trend is the assignment 
due to Gupta, Colvin and George [191] who notice that during thermal annealing, the 
peaks at 2087 cm -1 and 910 cm -1 exhibit an identical temperature dependence, where 
their magnitudes start decreasing from 550 K. Both peaks disappear completely by 
720 K. In isothermal annealing, the time evolution of the integrated absorbance of 
these two peaks is nearly identical. Therefore, they assign the 2087 cm -1 line to the 
Si-H 2 stretch mode, as it is well known that the 910 cm -1 peak corresponds to the 
Si-H 2 scissors mode. It was also found that the 2110 cm -1 peak initially increased 
in magnitude and eventually shifted down to 2102 cm -1 . The variations in the peak 
intensities were shown to be consistent with a hypothesized conversion of Si-H 2 to 
Si-H and H 2 . Hence, the 2110 cm -1 peak was assigned to Si-H. The 2140 cm" 1 
peak was assigned to the Si-H stretching in an O-Si environment. This assignment 
compares favorably with the predictions of Chabal and Raghavachari [216,217], if we 
assume that only the symmetric modes of Si-H and Si-H 2 stretching are observed in 
the infrared spectrum of porous silicon. 

The Si-H stretching mode frequencies for possible environments in porous silicon 
as predicted by Eq. (6.5) is reported in table 6.2. Based on these values, we attempt to 
reconcile the various systematics reported in different studies. Gupta and coworkers 
establish a strong correlation between the Si-H 2 scissors mode and the 2087 cm" 1 
line. However, the predicted increase of the stretching frequency with increasing 
hydrogen content also needs to be accounted for. We note from table 6.1 that the 
Si-H stretching frequency in crystalline silicon is reported [198] to be ~ 2000 cm" 1 . 
The frequency for Si-H 3 stretching mode calculated from Eq. (6.5) is ~ 2100 cm" 1 
(see Table 6.2). This suggests the assignment of 2087 cm -1 to Si-H 2 and 2110 cm" 1 to 
Si-H 3 stretching modes, which is seen to be consistent with the experimental results 
of Gupta, Colvin and George and the theoretical trend predicted by the present work 
as well as Ogata and coworkers [209]. The 2140 cm" 1 line can be assigned to an 
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Table 6.2: Stretching mode Si-H frequencies ( vsi-H ) f° r possible environments in 
porous silicon. The u Si - H have been calculated on the basis of Eq. (6.5). The symbol 


sents < 

No. 

l aangimg uuna. 

System 

usi-H (cm -1 ) 

1 

(Si)SiH 3 

2106 

2 

(OSi)SiH 2 

2146 

1 

3 

(OC)SiHo 

2175 ! 

4 

(SiOo)SiH 

2186 

5 

(F)SiH 3 

2223 j 

6 

(0 3 )SiH 

2273 

7 

(F 2 )SiH 2 

2293 

8 

(Od)SiH 2 

2052 

9 

(Fd)SiH 2 

20S2 

10 

(0 2 d)SiH 

2092 

11 

(F 2 d)SiH 

2152 


oxygen complex (OSi)Si-H 2 as calculated from Eq. (6.5) (see Table 6.2). Recall from 
our earlier discussion that a similar suggestion has also been made by Gupta and 
coworkers [191]. 

While the location of the hydrogen related peaks in the vibrational spectra of 
porous silicon samples agree with that of H passivated Si(100) surfaces, the peaks 
are substantially broader with typical FWHM of ~ 20 cm -1 as opposed to 1 cm -1 in 
crystalline Si. This may be due to the high degree of disorder present on the surface, 
as a result of surface contamination as well as geometrical relaxation of the surface 
Si atoms. Oxygen, fluorine and organic radicals can get introduced into the silicon 
system during the anodization of silicon in HF. Substituted oxygen and fluorine may 
give rise to i/si-H in the range 2050-2150 cm. -1 . Table 6.2 enumerates some of the 
possible environments with frequencies calculated using Eq. (6.5). The first four 
entries of vsi-h fall in the range 2050-2150 cm -1 or are close to it. It has been 
reported that the number of dangling bonds in porous silicon increases upon UV 
illumination [194]. Hence we have also considered dangling bond environments. The 
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l<ist foui entries in table 6.2 suggest four such environments which mav contribute 
to the observed v s ,-n hne (2050 - 2150 cm' 1 ) in porous silicon. Note that the 
substitution of an oxygen atom in (0 3 )SiH by a dangling bond reduces v Si - H from 
2273 cm -1 to 2092 cm. -1 Broad lines are also observed at ss 850 cm -1 and 1100 cm -1 . 
These lines have been assigned to Si-H 2 scissors mode and Si-O-Si stretching modes 
respectively in the past. However, Si-F related frequencies [209,218] also lie in the 
range 800 cm -1 to 1030 cm -1 . In view of the possible presence of fluorine in porous 
silicon, caution must be exercised in assigning frequencies to the observed infrared 
lines. 

The changes in the surface chemical environment can also lead to geometrical 
relaxation of the cluster. In the ab initio calculations of Chabal and Raghavachari 
[216,217], it was found that the Si-H distance decreases slightly because of the strain 
present in the enhanced Si-Si distances (Si-Si in the cluster was found to be 2.51 
A whereas the bulk value is 2.35 A). In the first principles calculations of Buda and 
coworkers [41], the surface Si-Si distances were found to be relaxed by 0.15 A, and 
hence d (Si-Si) = 2.5 A. It is therefore reasonable to assume that in porous silicon, 
the Si-H distances are somewhat lower. Using the phenomenological formula [197] 

i = 7074Acm~ 1 (6.7) 


we obtain, the variation in the vibrational frequency Au to be 



21222 V J 

— i Adsi-H 

USi-H 


For a value of Ad Si - H = 0.3 A, Av works out to be 1257 cm 1 . Thus a small variation 
in the Si-H distance can lead to a large variation in the vibrational frequency. This 
could lead to an upshift in the frequencies as well as a broadening of the peaks if there 
is a distribution of Si-H distances. Further, based on the observed range of vsi-H in 
porous silicon, Eq. (6.7) suggests that d Si -H is confined to {1.49-1.51}A. This should 
be a useful guide in electronic structure calculations where the dangling bonds are 
passivated by hydrogen. Values such as d$i-H — 1-637 A or dsi-H = 1-17A employed 
in these calculations [39,45] are clearly out of range. 
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6.4 Conclusion 

The present work postulates a renormalized electronegativity scale REEL based on 
the principal quantum number n and the orbital radius r s . Consequently, this new 
scale is purely non-empirical in character. A universal relationship between the 
Sanderson electronegativity and our quantum scale V defined by Eq. (6.4) has been 
established. This relation is seen to be valid for almost the entire periodic table. A 
detailed discussion of the renormalized electronegativity scale in the context of ex- 
tant phenomenological scales is to be found elsewhere [205]. In this chapter, we have 
highlighted the utility of our scale in assigning vibrational modes to Si-H stretching 
in porous silicon. This approach can be gainfully extended to study the infrared 
spectra of ligands other than hydrogen as well as other semiconductors. 

One may extend this work to v Si - H in silane molecules substituted by organic 
radicals such as CH 3 , C2H5, C 6 H 5 etc. The SR values for these radicals can be 
calculated using Lucovsky's prescription [197], 

SR= (f[SR(Rj)) 1/n 

j = 1 

where SR(Rj) represents the Sanderson electronegativity values of the component 
elements. The linear fit in Fig. 6.2 enables us to define an effective V(Rj) using these 
SR values. One may then look for a correlation between vsi-n and these values of 
V(Rj). We have carried out such an exercise for SiH(jRCl„) (n = 0.1.2) where R 
stands for the radicals and found an approximate linear relationship akin to Fig. 6.3. 

Chlorine is known to be a good dangling bond passivator in silicon. We have also 
discovered a linear relationship between the Si-Cl stretching mode frequency and our 
quantum scale. We use the data for v S i-ci in a-Si cited by Wu et al. [219]. Fig. 6.5 
depicts this correlation. Except for the highly electronegative environment (ClaJSiCl. 
the data is linear to a good approximation. 

Our model has focused mainly on the identification of Si-H modes in porous 
silicon. Hydrogen plays a key role in silicon. In c-Si, it is known to passivate both 
deep and shallow impurities [220]. In porous silicon, the presence of hydrogen is 
linked to surface passivation which facilitates efficient luminescence. We also show 
that this model can be used to study the vibrational modes of other impurities. 

The past assignments of the vibrational spectra of PS to specific transitions exhibit 
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Figure 6.5; Linear relationship between the Si-Cl bond stretching frequency usi-ci 

4 

and the quantum scale ]T]V(,Rj). 

j = i 

some inconsistencies. Our model has limited value in making exact identifications, 
but suggests systematic variations in the vibrational spectrum with the addition 
(or removal) of other elements to the Si-H environment. For example, Shanks and 
coworkers (212] state that the shift of the 2000 cm -1 line on annealing to 2100 cm -1 is 
due to the introduction of two dangling bonds. However, our work seems to indicate 
that the introduction of dangling bonds can only lower the vibrational frequency as 
mentioned earlier. Thus the trends suggested by the present work can be used to 
improve upon existing results. Again, our model contradicts some of the assignments 
due to Gupta, Colvin and George [191] who posit a decrease in the vibrational fre- 
quency as the number of neighboring hydrogen atoms are increased. The main utility 
of this model thus lies in correlating changes in the vibrational spectrum with specific 
changes in the surface environment, which are known to play a major role in porous 
silicon photoluminescence. The PL spectra of siloxene, a complex considered respon- ( 
sible for the luminescence in porous silicon, can also be varied by the introduction j 
of halides or organic radicals [193]. In porous silicon, post anodization surface treat- 
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ment modifies the PL spectra [221]. This is perhaps due to the formation of molecular 
complexes on the surface of porous silicon. Some insight into this phenomena can be 
obtained by studying the vibrational spectra. Our approach provides a valuable tool 
to correlate the PL spectra of various samples with their vibrational spectra. This 
can serve to characterize the porous silicon surface. 


Key Reference 

The major reference for the material contained in this chapter is: 

• George C. John and Vijav A. Singh, Phvs. Rev. B 53, 9831-9837 1996). 


Chapter 7 


Conclusion 


Traditionally, the physicists have adopted a reductionist approach to materials 
physics, focusing on the extrapolation of individual atomic properties to macroscopic 
collections of atoms [222]. The quantum chemists carry out electronic structure cal- 
culations on molecules consisting of < 10 2 atoms. Solid state physicists compute 
the band structure of infinite solids by applying periodic boundary conditions on a 
unit cell consisting of < 10 atoms. Over the last decade, the increasing technolog- 
ical importance of nanocrystalline materials have called for a detailed, microscopic 
analysis of atomic clusters of the size 10 3 — 10 4 atoms [223]. Such clusters are found 
to exhibit properties surprisingly distinct from the molecular and bulk limits. This 
regime, bridging the microscopic (< 10 2 ) atoms and the macroscopic (> 10 6 ) atoms, 
has been termed mesoscopic. (See Fig. 7.1). 

Porous silicon, the material that we have examined in some detail in this thesis, 
belongs to this category of nanocrystalline materials. However, unlike the III-V 
quantum wells or dots which are normally fabricated under controlled conditions using 
sophisticated techniques such as molecular beam epitaxy (MBE) or metal organic 
chemical vapor deposition (MOCVD), porous silicon can be fabricated easily and ; 
economically - in fact in most undergraduate laboratories. In this sense, porous ; 
silicon could be described as the “poor man’s quantum dot.” j 

Porous silicon is easily fabricated. Its PL is present for all to see. But the agony j 
of the situation is not merely confined to the lack of success in presenting a unified 
model for the optical behavior. The holy grail of porous silicon research, namely j 
a commercially viable electroluminescent device continues to remain elusive. The j 
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Figure 7.1: Transition from the molecular regime to the bulk solid. Nanocrystallites 
and quantum dots fall in the intermediate mesoscopic region. The shaded area de- 
picts the size range of luminescing crystallites in porous silicon. Also shown are the 
schematic variation in properties. The mass of the free electron is represented by m e . 
The diameter d of a cluster of N atoms is d = (3.V/47r) 1 / 3 ao where ao = 0.543 nm. 
the silicon lattice constant. Uncertain values are indicated by question marks. 
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l&rgt? luminescence deca\ times (10 s) may render it unsuitable for optoelectronic 
applications requiring smaller response times. Efforts to obtain efficient electrolu- 
minescence have been disappointing and we face a technological bottleneck. At the 
time of writing this thesis, there has been reports of the successful fabrication of 
an integrated optoelectronic circuit incorporating porous silicon [224]. Though the 
speed and efficiency of this circuit needs to improve further, this achievement seems 
to hold out some promise for the future of porous silicon. 

The lack of control over the formation process leads to a high degree of disorder in 
the system as mentioned earlier in the thesis. To obtain a reasonable theory of porous 
silicon, one needs to look beyond the properties of individual silicon nanocrvstallites. 
The hierarchical disorder (see Sec. 1.4) has to be accounted for, and one should at- 
tempt to study the "collective microstructure” of the system consisting of crystallites 
of various sizes and shapes. 

Phenomenological theories for porous silicon formation employ experimentally fa- 
miliar parameters such as the anodization potential, electrolyte concentration etc., 
to qualitatively argue the case for pore formation in silicon. On the other hand, com- 
puter simulation models of chapters 2 and 3 employ rules of growth which are not- 
explicitly related to the experimental growth parameters. Thus, these theories are 
in a sense "indirect." The drift-diffusion model described in chapter 2 successfully 
reproduces the various morphological transitions observed in porous silicon with the 
tuning of experimental parameters. The PS formation problem is thus identified as 
one belonging to the larger class of non-equilibrium growth phenomena, where the 
competition between two antithetical parameters give rise to various morphological 
classes. It would be worthwhile to explore this further so that the competing phenom- 
ena mav be identified and their relevance to the simulation parameters established. 

The diffusion induced nucleation model for porous silicon formation discussed in 
chapter 3 successfully explains (i) a constant rate of growth, (ii) the dependence of 
the rate of growth on the anodization potential, (iii) high porosity structures similar 
to samples exhibiting visible photoluminescence and (iv) electropolishing of silicon in 
the high potential limi t. It would be desirable to carry out further studies using this 
model to establish a proper correspondence between the parameters of the simulation 
and the experiment. One needs to explain the occurrence of the constant porosity 
profile as well as the nearly planar film front. 
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A desirable goal is the formulation of a computational model which can simulate 
the fabrication process on a computer. Simulation studies of molecular beam epitaxy 
(MBE) growth of AlGaAs/GaAs semiconductor wells is a good example [225]. The 
DIN model discussed in chapter 3 can be used to study the effect of substrate rough- 
ness on porosity. The effect of placing inert masks on the substrate surface during 
the dissolution process can also be explored. Thus, the model can be extended to 
design experiments to fabricate porous silicon layers of varying geometry. 

The porous silicon community continues to remain divided as to the exact nature 
of the radiative process. Of the three categories of models described in Sec. 1.3, 
namely (i) the quantum confinement model, (ii) the chemical luminescence model, 
and (iii) the hybrid model, the quantum confinement hypothesis seems to be the 
simplest and the most appealing. Recently, the quantum confinement model has 
been reinforced by two experiments. Schuppler and coworkers [57] have obtained 
visible photoluminescence from Si nanocrystallites of about 1.5 run diameter. They 
contend that the luminescence in porous silicon is due to dots /columns of size <1.5 
nm, as opposed to the earlier belief that the sizes are typically 3 nm. In Si/SiC >2 
superlattices grown using molecular beam epitaxy, an unambiguous blue shift in 
the PL peak position with decreasing layer thickness in the range d ~ 1-5-3 nm 
has been observed by Lockwood and coworkers [17], Further, employing soft X- 
ray edge absorption spectroscopy and X-ray photo emission spectroscopy, they have 
measured the conduction band upshift and the valence band downshift respectively. 
For d= 1.5 nm, the valence band downshift is 0.1 eV and the conduction band upshift 
is 0.3 eV. Note that these measurements are in conflict with the effective mass based 

calculations described in Sec. 1.3.1, which predict a larger shift for the valence band 
[19]. 

The influence of the surface in the recombination process cannot be ruled out. 
Consider the QCM where an electron is radiatively excited into the conduction band. 
The experimental luminescence time constant r > 1 /j.sec. A simple particle in a box 
picture indicates that within this time the electron can sample the nanocrystallite sur- 
face roughly ten million times. Alternatively, we know that the electron mobility //„ 
in bulk silicon is 500 cm 2 /volt-sec. From the Einstein relation the diffusion constant 
D n = UnkT/q = 12.5 cm 2 1 sec. Given the crystallite size of 3 nm, the transit time 
to the surface is in the femtosecond range and much smaller than the luminescence 
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Figure 7.2: Role of surface states in the luminescence process. Fig. (a) depicts the 
schematic structure of a crystallite. The electron hole pair (ellipses) can get trapped 
in the surface states depicted by the shaded regions (solid arrow). Possible recom- 
bination pathways are indicated by the broken arrow. Fig. (b) shows the density of 
states n(E) in a realistic nanocrystallite, with a U shaped distribution of surface states 
S near the band edges CB (conduction band) and VB (valence band). The peak near 
the mid-gap region is due to dangling bonds (D) (also see Sec. 1.3). According to 
the Berthelot hopping model discussed in chapter 5, the surface states could act as 
pathways for carrier hopping. 

time. 

Surface related chemical mechanisms would be a cause for concern where the ap- 
plications of PL are involved, and perhaps, the scientific . community wants a more 
"physical theory/’ The presence of quantum confinement is clearly established by the j 
shift in the absorption edge from 1.1 eV for bulk silicon to over 1.6 eV for a lumi- 
nescing porous silicon layer. Hence, the variety of hybrid models which invoke both ; 
quantum confinement and the surface chemistry, seems an attractive compromise. 
These models assign roles to the core, surface, and dangling bond related states, but j 
differ in the identification of radiative transitions responsible for the visible PL (see j 
Fig. 7.2). The possible transitions are many - band to band (CB — > VB), bound j 
state to band (S/DB -4 VB, CB -4 S/DB) or bound state to bound state (S/DB -4 
S/DB). Hence these models do not come up with specific predictions, are vague, and 
flexible enough to accommodate experimental contingencies. In short, these theories j 



122 


Conclusion 


are non-falsifiable and hence not -‘scientific” in Karl Popper’s sense of the term [226]. 

From the course of the discussions on photoluminescence, it is clear that the 
source of luminescence consists of nanocrystalline structures which could be as small 
as 1 nm. This translates to around 10 2 - 10 3 atoms per crystallite (see Fig. 7.1). 
This is the regime where the molecular calculations of quantum chemistry meet the 
traditional solid state calculations. The question now arises whether one should view 
the problem of porous silicon in a reductionist fashion, i.e.. as a simple extrapola- 
tion of the properties of a constituent nanocrystallite. This approach has enjoyed a 
degree of success as elucidated in chapter 4. In this chapter, we have shown that 
the assumption of a Gaussian distribution of crystallite sizes can successfully explain 
the observed shape of the emission spectrum. This calculation also explains the dis- 
parities in the excitonic binding energies which need to be assumed by electronic 
structure calculations reviewed in Sec. 1.4. Thus, the necessity of incorporating the 
disorder inherent to the porous silicon system is highlighted. 

The spectral linewidth is an important parameter in the fabrication of light emit- 
ting devices. Chapter 4 explains the observed line width as a simple sum of the 
contributions due to individual crystallites. The carrier hopping length scale S of 
the model in chapter 5 is of the order of the crystallite size (1 — 5 nm.). Workers 
have also tried to model carrier transport using a variable range hopping model (see 
Sec. 5.5). Hence the calculations of chapter 4 can be improved upon by accounting 
for the collective behavior of the crystallites. Recently, there have been reports of a 
reduction in the PL line width due to electronic coupling in a stack of indium arsenide 
quantum dots separated by gallium arsenide layers [227,228]. One needs to look for 
similar phenomena in porous silicon. These considerations assume importance if one 
attempts to fabricate semiconductor lasers from porous silicon. 

A consequence of the ongoing debate about the photoluminescence mechanism has 
been the accumulation of a mountainous volume of experimental data on the optical 
behavior of silicon. The results from various experimental groups, which have often 
been advanced in support of a particular hypothesis, need to be explained under a 
common “unified” model. Chapter 5 presents a phenomenological scheme which lays 
the foundations for such a unification. Our model invokes the competition between (i) 
an Arrhenius type radiative, process and (ii) a Berthelot type carrier hopping process 
to explain wide ranging phenomena such as. 
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• The variation of luminescence intensity with temperature, with a maximum in 
the range 50 — 150 K. 

• The exponential drop in intensity with increasing pressure. 

• The faster luminescence decay with increasing temperature. 

Assuming a quantum confinement model (see Fig. 5.1), one can extend the above 
phenomenological scheme to explain 

• The variation of the luminescence decay time with emission energy. 

• The anomalous variation in the PL peak position with temperature, which 
shows both blue and redshifts. 

• The dependence on the PL peak energy on applied pressure as well as the 
pressure transmitting medium. 

Thus our approach is seen to support the quantum confinement hypothesis. This 
model can also explain, to some extent, the high surface sensitivity of PL as the 
Berthelot process is most likely to be a surface related mechanism. For instance, the 
carrier hopping could take place between the surface states in the gap (see Fig. 7.2), 
whereas the light emission occurs as a consequence of re-excitation into the bands 
and subsequent band to band recombination. 

There does exist several observations that fall outside the ambit of the model 
proposed in chapter 5. The disparities in the absorption and emission edges cannot 
be addressed within our scheme. The electroluminescence from PS is found to be ex- ; 
tremelv weak. This can be attributed to a surface dominated transport process which 
results in poor carrier injection into the crystallites. The model should be extended : 
to address transport related phenomena such as persistent photo conductivity [229]. ; 

The exact nature of the vibrating barrier assumed in chapter 5 is still unclear, j; 
However, it is well known that several vibrating species are present on the PS surface, j 
the frequencies of which have been linked to a theoretical electronegativity scale in ! 
chapter 6. Using the renormalized electronegativity scale defined in chapter 6, one j 

<• j 

may attempt to predict the possible environments that can give rise to vibrating 
barriers in the desired frequency range. However, such an approach needs to be | 
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augmented by appropriate experiments to vary the surface composition of the silicon 
crystallites in a controlled manner. 

The main utility of the phenomenological model outlined in chapter 6 lies in cor- 
relating changes in the vibrational spectrum to changes in the surface chemistry. 
This scheme can be used to predict the surface composition by observing the infrared 
vibrational spectra after post-anodization treatments. The variations in the lumines- 
cence behavior can then be attributed to the presence (absence of specific chemical 
species on the surface. 

Summarizing, we have, in this thesis 

• Proposed a generalized computational scheme to model aggregation and etching 
phenomena. 

• Successfully extended this model to explain porous silicon formation during the 
anodization of silicon in HF. 

• Proposed a theoretical framework to explain the observed spectral shape in 
porous silicon photoluminescence. This framework provides two parameters, 
the mean diameter d 0 and variance a of the nanocrystallite size distribution 
which can perhaps be used to characterize the porous silicon sample. 

• Proposed a novel phenomenological model to explain a large body of data on 
the optical behavior of porous silicon under a “unified theory.” This model 
suggests a valuable characterizing parameter. Tg, the Berthelot temperature, 
which is also experimentally easily accessible. 

• Defined a novel quantum scale which is then successfully employed to predict 
the vibrational spectra of defects in porous silicon. 

The computational models of porous silicon formation can be gainfully extended 
to aid device design and fabrication. The phenomenological models of chapters 4. 5 
and 6 suggest novel characterizing scales and provide a firm mathematical basis for 
predicting the optical behavior as well as the vibrational spectra of porous silicon. 
This can constitute a basis for further microscopic calculations which need to take 
into account the mesoscopic, disordered nature of the system. 
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